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FOREWORD 


This  rsport  was  prapared  by  th©  Douglas  Aircraft  Coinpanyj  Aircraft 
Division,  Long  Beach,  California,  for  the  Analysis  and  Optimization  Branch, 
Structural  Mechanics  Division,  Air  Force  Flight  Dynamics  Laboratory,  Air 
Force  Wright  Aeronautical  Laboratories,  Wright-Patterson  Air  Force  Base,  Ohio 
under  contract  F33615-77-C-3101 .  This  research  was  conducted  under  Project 
2401,  "Structural  Mechanics",  and  Task  240102,  "Design  Analysis  Methods  for 
Aerospace  Vehicle  Structures".  Mr.  L.  J.  Huttsell  of  the  Structural  Mechanics 
Division  was  the  Project  Engineer. 

The  work  reported  herein  was  conducted  during  the  period  of  August  1,  1977 
to  March  1,  1978.  _ 

The  Principal  Investigators  were  J.  A.  McGrew,  J.  P.  Giesing  and  R.  M. 
Pearson  assisted  by  K.  Zuhuruddin,  M.  E.  Schmidt  and  T.  P.  Kalman. 

The  preliminary  theoretical  developments  necessary  to  incorporate  part 
of  the  improved  aerodynamic  time  lag  representations  were  developed  under 
IRAD  funding  prior  to  this  contract  award.  Sections  2. 5. 1.1  through  2. 5. 1.4 
and  Appendices  A  and  B  document  these  efforts  for  the  purpose  of  completeness. 
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SECTION  I  INTRODUCTION 


Modern  aircraft  wing  designs  are  increasingly  taking  advantage  of 
recently  developed  supercritical  airfoil  sections  which  are  more  aerodynami- 
cally  efficient  than  conventional  airfoil  sections.  This  increased  efficiency 
is  achieved  by  delaying  the  drag  rise,  thus  allowing  the  aircraft  to  cruise 
at  higher  Mach  numbers  with  less  sweep.  Also,  higher  lift  coefficients  are 
obtainable  with  wings  using  these  advanced  airfoil  sections.  The  YC-15  and 
YC-15II  prototype  aircraft  are  recent  examples  using  supercritical  wings. 

There  has  been  concern  that  the  flutter  characteristics  of  supercritical 
wings  may  be  significantly  different  from  wings  with  conventional  airfoils. 
Recent  NASA  test  data  on  a  TF-8A  flutter  me  del  showed  a  substantially  larger 
decrease  of  flutter  speed  near  the  Mach  number  of  maximum  lift  curve  slope 
than  predicted  by  theoretical  subsonic  oscillatory  aerodynamic  theory. 

Further  studies  were  needed  to  determine  the  aerodynamic  mechanisms  which 
contribute  to  this  large  decrease  in  flutter  speed. 

The  objective  of  this  study  was  to  develop  and  evaluate  transonic 
flutter  prediction  methods  for  supercritical  wings.  The  study  was  divided 
into  five  phases: 

a.  Phase  I  -  TF-8A  Flutter  Analysis  (Steady  Corrections)  - 

Conduct  a  flutter  analysis  of  the  TF-8A  flutter  model  using  correct¬ 
ion  factors  to  the  aerodynamic  influence  coefficients  (AIC's) 
based  upon  steady  aerodynamic  data  from  the  pressure  model  and 
correlate  with  test  results. 

b.  Phase  II  -  TF-8A  Flutter  Analysis  (Unsteady  Correction)  - 

Apply  an  improved  unsteady  correction  factor  technique  (see  Phase 
III)  to  the  flutter  analysis  of  the  TF-8A  flutter  model  and  corre¬ 
late  with  test  results. 

c.  Phase  III  -  Transonic  Method  Improvement 

Improve  the  time  lag  representation  and  include  shock  wave  motion 
in  the  transonic  method. 
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d.  Phase  IV  -  YC-15II  Flutter  Analysis  - 

Use  the  final,  proven  methods  of  Phases  I,  II  and  III  to  perform 
a  flutter  analyses  of  the  YC-15I1. 

e.  Phase  V  -  Flutter  Test  Plan  - 

Use  the  results  of  the  analytical  effort  to  define  an  economical 

flutter  model  test  plan. 
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SECTION  II  TECHNICAL  DISCUSSION 
2.1  BASIC  AERODYNAMIC  DATA 

The  basic  theory  used  for  the  flutter  and  static  aeroelastic  analyses 
reported  herein  was  the  Doublet  Lattice  Method  (DLM)  described  in  Reference 
1.  This  method  is  applicable  to  general  configurations  and  is  well  proven 
and  correlated.  The  DLM  requires  that  the  lifting  surfaces  (wings,  tails, 
etc.)  be  broken  up  into  boxes  and  the  bodies  (fuselage,  nacelles,  etc.)  be 
idealized  with  both  slender  body  and  interference  elements.  Figures  2.1-1 
and  2.1-2  present  the  idealizations  for  the  two  configurations  considered 
in  this  report,  i.e.,  the  TF-8A  wing-body  and  the  YC-15II. 


Experimental  wind-tunnel  data  are  available  for  both  of  these  configura¬ 
tions.  The  YC-15II  data  were  taken  from  unpublished  McDonnell  Douglas 
sources  while  the  TF-8A  data  were  taken  from  Reference  2.  The  TF-8A  wind  tunnel 
model  was  found  to  be  quite  elastic  considering  that  it  was  to  be  a  nominally 
'rigid'  model.  Aeroelastic  corrections  were  applied  to  the  data  in  an 
attempt  to  remove  the  effects  of  flexibility.  These  corrections  were  essential 
to  obtain  the  correct  section  lift  curve  slopes  across  the  span  and  consisted 
of  determining  the  actual  o**  'effective'  angle-of-attack,  of  each 

section  and  using  these  in  the  distributions  of  c^and  c^^.  The  aeroelastic 
correction  was  based  on  a  bending  curve  found  in  Reference  2  for  one 
particular  loading.  This  curve  was  generalized  for  all  cases  by  normaliza¬ 
tion  with  respect  to  the  angle-of-attack  and  dynamic  pressure.  The  twist  for 
the  reference  condition  was  obtained  from  the  bending  curve  by  taking  the 
slope  and  scaling  it  with  the  sine  of  the  sweep  angle. 


e  =  (dz  /  dyjsin  A, 
e 

r 

% 

where  z  is  the  local  bending  deflection  and  y  is  distance  along  the 


elastic  axis. 


The  resulting  normalized  elastic  twist, 1^.  is  given  as: 


T 

a 


r  r 


and  is  plotted  in  Figure  2.4-4.  Subsection  2.4  presents  additional  details 
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1 

I 

! 

of  the  aeroelastic  correction  procedure  including  the  corrected  or  effective 
angle-of-attack.  at  each  spanwise  section  which  is  given  by: 

Ct  rr  ~  Cl(l“T  g) 

eff  a 

Subsection  3.1.3  presents  plots  of  section  c^  and  c^^  data  versus 
for  various  conditions. 

The  data  of  Reference  2  were  obtained  for  a  low  Reynolds  Number, 
specifically  1.96  x  10^,  based  on  the  chord  at  the  three-quarter  span  loca¬ 
tion.  Thus,  the  data  is  not  considered  reliable  for  predicting  the  stall  or 
near  stall  conditions.  Because  of  this, only  lift  curve  and  moment  curve 
slopes  in  the  low  angle-of-attack  range  were  used,  i.e.,  where  c^  was  a 
near-linear  function  of  The  moment  curves  were  not  linear  anywhere  but 

the  non-linearity  in  the  small  angle-of-attack  range  was  due  to  transonic 
effects  and  thus  were  considered  more  reliable.  The  points  on  the  lift  and 
moment  curves  where  the  slopes  were  taken  are  given  in  the  results  section 
(Section  3). 
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2.2  VIBRATION  MOOAL  DATA 

Normal  modes  of  vibration  of  the  TF-8A  wing  and  the  YC-15II  wing  were 
used  as  the  basis  for  flutter  and  static  aeroelastic  calculations.  For  the 
YC-15II  modes  were  available  from  the  design  analyses  process  and  were  used 
directly.  These  modal  data  had  been  correlated  with  past  aircraft  vibration 
test  data  and  good  agreement  shown. 

NASA  (Langley)  provided  measured  mode  shapes  and  generalized  mass  for 
six  elastic  modes  of  the  TF-8A  model  wing  In  the  form  of  normal  modal  defl¬ 
ections  at  81  wing  locations.  Spanwise  sectional  plots  of  these  deflections 
showed  relatively  small  camber  deformations  In  the  modes  of  concern  and 
therefore  the  modal  amplitudes  were  transformed  Into  sectional  pitch  and 
plunge  degrees  of  freedom  In  order  to  simplify  the  subsequent  aerodynamic 
weighting  processes. 

A  sectional  least  squares  fit  was  used  to  determine  pitch  and  plunge 
amplitudes  for  a  sixteen  spanwise  section  mathematical  model.  The  measured 
modal  data  and  transformed  modes  are  shown  In  Section  3.1. 

As  a  check  against  the  transformation  process,  vibration  modes  for  the 
same  (sectional)  model  were  calculated  with  distributed  stiffness  and 
inertial  data  which  become  available  later  In  the  study  program.  The  result 
ing  mode  shapes,  frequencies  and  generalized  mass  agreed  well  with  the  trans 
'^formed  test  data  except  for  the  torsion  mode,  in  which  the  discrepancy  is 
due  to  a  probable  misinterpretation  of  torsional  stiffness  magnitude. 
Factoring  the  torsional  stiffness  levels  led  to  good  agreement  In  modal 
deflections  and  generalized  mass. 


2.3  STATIC  AEROELASTIC  ANALYSIS  PROCEDURES 

The  analysis  procedures  used  for  all  static  aeroelastic  analyses  proced¬ 
ures  used  herein  consist  of  a  subset  of  the  generalized  equations  used  to 
solve  for  the  linear  elastic  trim  conditions  and  deformation  of  a  free 
flying  aircraft.  It  is  assumed  that  total  balance  force  of  the  model  is 
specified  and  the  wind  tunnel  model  is  then  pitched  until  that  specified 
total  lift  is  achieved.  Drag  forces  are  not  accounted  for  herein  although 
the  procedure  is  simply  extended  for  these  forces  if  required.  The  number 
of  equations  to  be  solved  for  is  n  +  1,  where  n  is  the  number  of  normal 
elastic  modes  of  the  system  used.  Lift  and  weight  are  taken  as  positive 
down  and  moments  are  positive  nose  up  in  the  following.  The  generalized 
equations  are  formed  by  a  similarity  transform  with  the  normal  modes  upon 
the  mass  and  aerodynamic  equations  for  the  model .  The  total  force  on  the 
model  support  is  given  by 

^  {q})  +  q  Lj  +  W  =  (2.3-1) 

The  elastic  equations  are: 

[K^^]  {q}  -  q  ({D^^}a  +  [D^^]  {q})  -  gim^^}  -  q{Lj)  =  D  (2.3-2) 
where 

q  =  dynamic  pressure,  psi 

a  =  model  reference  angle  of  attack 

q  =  generalized  modal  amplitudes  (  an  n  x  1  vector) 

[.  =  total  jig  (reference  integrated  basic  lift),  lift/q 

J 

W  =  total  weight 

Fj  =  total  specified  force  on  the  modal  balance 

2 

=  generalized  stiffness  (m^^  t  u)  ^^) 
g.mEh  =  generalized  weight  forces 

L I  =  generalized  jig  (reference  basic  lift  and  moment),  forces/q 
J 

D.  =  total  lift  due  to  a/q 

ha 
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D,.  =  elastic  aerodynamic  force/a*q 

Ea 

=  total  lift/q*q 
=  elastic  aerodynamic  force/q*q 

The  model  reference  angle  of  attack  is  split  into  two  parts  for  convenience: 
a  = 

where 

a  =  model  insertion  angle 
0 

oj  =  model  trim  angle 

Equation  (2.3-1)  may  be  solved  for  in  terms  of  q,  jig  forces  and 
balance  force  and  substituted  into  Equation  (2.3-2)  leaving  n  equations  with 
n  unknowns  (q).  Post  multiplication  of  the  vector  q  into  the  modal  set  used 
to  obtain  the  equations  leads  to  the  deformed  surface 

{h}  =  [<!']  {q}  (2.3-3) 

where  h  are  the  deflections  (and  rotations)  of  the  surface  and  $  is  the  set  of 
normal  mode  shapes  of  the  surface.  The  spanwise  local  angle  of  attack 
distribution  then  is  given  by  the  sum  of  the  model  Insertion  angle,  the  model 
trim  angle,  the  local  elastic  twist  angles  and  the  effective  jig  or  built  in 
twist  distribution. 


2.4  AERODYNAMIC  WEIGHTING  FACTORS 

The  concept  of  modifying  the  matrix  of  theoretical  Aerodynamic  Influence 
Coefficients  [AIC's]  with  a  correction  factor  or  weight  factor  matrix  [WT] 
has  been  investigated  in  many  reports  one  of  which  is  Reference  3.  One  of 
the  methods  of  Reference  3  is  to  premultiply  the  theoretical  AIC  matrix 
with  the  weight  factor  matrix  as  follows: 

[AIC]^,.^=  [WT]  [AIC] 

This  technique  is  the  one  used  in  the  present  report. 

The  type  of  AIC  matrix  used  is  that  for  an  elastic  axis  representation. 

Thus,  the  elements  of  the  AIC  are  proportional  to  the  sectional  values  of 
lift  and  moment  coefficients.  In  what  is  to  follow  the  sectional  lift  co¬ 
efficient  is  discussed  primarily;  however,  the  concept  applies  equally  to 
the  section  moment  coefficient. 

The  weight  factor  for  a  particular  section  lift  coefficient  (at  a  parti¬ 
cular  wing  spanwise  location)  is  given  in  two  parts:  a  steady  part,  WTS 
and  an  unsteady  part,  WTUS.  The  total  weighted  value  for  lift  is: 

c.  (k  )  =  WTSCL  WTUSCL  c^  (k^) 

a  ^  ct 

where  WTSCL  and  WTUSCL  are  WTS  and  WT'JS  for  the  c^^  term. 

a 

The  steady  weight  factor,  WTS,  is  usually  obtained  from  wind  tunnel 
section  lift  coefficient  data.  The  idea  is  to  correct  the  theoretical  lift 
coefficient  so  that  it  gives  the  experimental  values  when  the  frequency  is 
zero.  This  term  is  obtained  from  the  ratio  of  experimental  to  theoretical 
lift  curve  slope  at  k^  =  0.  Thus: 

WTSCL  =  c^  (k^  =  0)/cj^  (k^  =  0) 

a  a 

Where  c  is  the  theoretical  result  and  c.  is  that  obtained  experimentally. 

a  “ 
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Similarly  for  the  moment 


WTSCM  =  (k  =  0)/c^  (k^  =  0) 

o  Cl 

Or,  in  terms  of  aerodynamic  center,  and  lift  coefficient. 


(exp) 

a 


WTSCM  >  c,  (Xg, 

a 


where  x  is  the  axis  about  which  the  moment  is  taken, 
ea 


The  aerodynamic  theory  Ubod  to  obtain  the  AIC  and,  consequently,  the 
c  (k  )  is  a  subsonic  lifting  surface  theory  containing  no  transonic  effects. 
Th0s,  the  effects  of  shock  motion,  etc.,  are  not  contained  in  this  theory. 
Transonic  effects  can  be  included  into  the  analysis  by  using  an  unsteady 
weight  factor  based  on  other  theories  which  contain  transonic  effects.  For 
instance,  if  a  two-dimensional  transonic  theory  is  available,  then  it  can  be 
used  as  a  basis  for  the  weight  factor,  WTUS,  on  a  strip  theory  basis  as 
follows: 


WTUSCL  (k^)  » 


(2DTRANS) 

“(MrANsT 


0) 


where  c  (2DTRANS)  result  obtained  from  the  two-dimensional  transonic 

theory  anS  c  is  the  result  obtained  from  the  two-dimensional  subsonic 

theory  (the  tfto-dimensional  equivalent  of  the  AIC).  The  reason  for  the 
division  by  the  values  of  c^  at  k^  =  0  (steady  values)  is  so  that  the  weight 
factor  WTUS  will  reduce  to  uPlity  for  the  steady  case.  This  is  required 
since  WTS  has  already  taken  care  of  the  steady  weighting. 


The  formula  above  for  WTUSCL  holds  equally  well  for  the  moment;  simply 

replace  WTUSCL  with  WTUSCM  and  c^^  with  c^  . 

a  a 

The  application  of  two-dimensional  theory  to  the  three-dimensional  swept 
wing  case  requires  further  analysis.  Consider  Figure  2.4-1. 
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Figure  2.4  -  1 


Strip  Theory  Approximation 
For  a  High  Aspect  Ratio  Wing 


The  constant  chord  two-dimensional  swept  wing,  that  approximates  the 
strip  on  the  wing  whose  chord  length  is  c,  is  shown  as  dashed  lines.  The 
flow  is  two-dimensional,  in  a  plane  normal  to  the  axis  of  sweep.  If  we 
designate  all  values  given  in  this  reference  system  with  a  superscript  tilda, 
then  the  section  lift,  moment  and  frequency  are: 


4Cpdx  =  Cj 


C  c 

=  -2  f  ''Cp(x-xjdx  -  J 


1C „(x-x„)dx  =  c„ 
p  m  m 


U)C  (jC  cos  A 

‘>•20  ' 


where  is  the  frequency  to  use  in  the  two-dimensional  theory.  The  two- 
dimensioSSl  frequency  is  related  to  the  three-dimensional  theory  as  follows: 


(oC  C  _ 

^  ?  " 

“  c 


£ 

C 


The  two-dimensional  transonic  theory  used  as  a  basis  for  this  strip  theory  is 
presented  in  subsections  2.5  and  3.3.  The  Mach  Number  used  in  the  two- 
dimensional  theories  is  related  to  the  free  stream  Mach  Number  as  follows: 


^20  “  ^ 

In  the  strip  theory  approach  the  chord  line  of  the  airfoil  section  considered 
is  forshortened  as  Figure  2.4-1  shows.  The  thickness  and  camber  distributions 
however  are  not  changed  thus  the  effective  thickness  and  camber  of  the  airfoil 
section,  as  a  percent  of  local  chord,  increases  as  1/cosA.  Also  the  angle  of 
attack  increases  in  the  same  way. 


The  steady  weight  factors,  WTS,  are  usually  determined  from  wind  tunnel 
data  taken  on  semi  rigid  models.  It  has  been  found  that  the  steady  incidence 
angle  of  each  wing  section  has  an  important  bearing  on  both  the  steady  and 
unsteady  lift  curve  slopes  especially  for  supercritical  wings.  Since  the 
shock  moves  as  the  section  angle  of  attack,  a,  is  changed,  there  appears  a 
large  increase  in  perturbation  pressure,  Cp^,  at  the  steady  shock  location 
which  is  proportional  to  the  jump  ifi  steady  pressure.  Figure  2.4-2  which  was 
taken  from  Reference  4,  illustrates  this  effect.  The  effect  of  the  peak  in 
pressure  is  very  noticeable  in  the  moment  coefficient  especially  for  super¬ 
critical  type  wing  sections. 


Figure  2.4-3  presents  the  lift  and  moment  coefficient  for  the  TF-8A  wing 
(y/(b/2)  =  .804)  operating  at  a  free  stream  Mach  Number  of  0.99.  Notice 
especially  the  nonlinearity  in  the  moment  coefficient.  Figure  2.4-3  also 
shows  nonlinearities  caused  by  viscous  or  separation  effects.  These  are  easily 
seen  in  the  plot  of  c^^  versus  a. 

The  slopes  of  these  curves,  ,  c^^^  ,  are  the  parameters  required  for 
the  weight  factor  calculation.  It“is  oBvious  that  the  slopes  vary  along  the 
curves.  Sometimes  this  variation  is  very  large.  Thus,  the  values  of  c^^  and  c^ 
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Figure  2.4-2  Steady  and  Quasi-Steady  Pressure  Distributions 
on  the  NLR  7301  Airfoil  in  Transonic  Flow 
(Taken  from  Reference  4  ) 
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used  depend  to  a  large  degree  on  the  section  angle  of  attack.  The  section 
angle  of  attack,  a^,  is  made  up  of  the  following  parts: 


which  consist  of  the  aircraft  angle  of  attack  a  (which  is  a  function  of 

weight,  etc.),  the  built-in  twist,  e.,  and  the  elastic  twist,  c  ,  (which  is 

J  ^ 

proportioned  to  loading,  dynamic  pressure,  etc.).  Once  the  section  twist  is 

known,  the  values  of  c„  and  c  can  be  determined  from  the  curves  obtained 

SL  m 

from  the  wind  tunnel.  ^  “ 

This  procedure  is  slightly  different  from  the  usual  procedure  of  reading 

c„  and  c  on  the  linear  part  of  the  curve.  The  usual  procedure  is  consider- 
£  m 

ed“conservStive  because  the  linear  value  of  c  is  most  often  the  maximum 
value,  resulting,  usually,  in  the  lowest  flutter  speed. 

Some  caution  must  be  exercised  when  using  to  find  the  lift  and  moment 
curve  slope.  The  reason  for  this  is  that  the  experimental  curves  of  and 
C|^,  versus  a  may  be  inaccurate  or  inappropriate  for  several  reasons.  First, 
the  pressure  model  is  not  entirely  rigid  and  model  bending  and  twisting  induce 
elastic  twist  in  the  model.  For  instance,  the  elastic  twist,  e^,  induced 
by  the  bending  of  the  swept  TF-8A  wing  can  be  calculated  from  bending  data 
given  in  Reference  2.  Figure  2.4-4  presents  a  plot  of  normalized 
elastic  twist  Cg/qa  versus  span  for  the  one  condition  ^t  which 

the  bending  was  determined.  This  curve  was  used  to  correct  the  section  angle- 
of-attack  to  an  effective  angle  of  attack,  for  all  cases  as  follows: 

eff  e 


where 


e  =  -T  a  q 
e  a 


or 


off  =  “  -  T„  ,) 
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where 


and  where  the  reference  condition  is  given  by  =  3.24'^  =  932  PSF, 

and  the  free  stream  Mach  Number,  M  =0.99. 

00 

Just  as  small  changes  in  steady  section  mean  angle  of  attack,  a^,  change 
the  loading  on  the  flutter  model,  small  changes  in  on  the  pressure  model 
have  the  same  effect.  Thus,  elastic  twist,  although  usually  assumed  small 
for  the  pressure  model,  wall  induced  angle  of  attack,  etc.  are  important 
considerations  when  using  the  pressure  model  data. 

Another  important  consideration,  especially  for  supercritical  wing 
sections,  is  the  camber  change  due  to  boundary  layer  thickness.  Shock 
location  and  the  loading  on  the  aft  part  of  the  supercritical  wing  section 
are  sometimes  very  critically  affected  by  the  boundary  layer  as  Figure 
2.4-5  shows.  This  indicates  that  the  boundary  layer  must  be  modeled  accurate¬ 
ly  in  the  wind  tunnel.  Therefore,  among  other  things,  the  test  Reynolds 
Number  must  be  sufficiently  large. 

Another  reason  for  requiring  that  the  test  Reynolds  Numbers  be  large  is 
that  the  proper  stall  characteristics  may  need  to  be  simulated.  Flutter 
models  are  usually  highly  flexible  so  that  flutter  speeds  can  be  reached  in 
the  wind  tunnel.  These  models  are  dynamically  modeled  but  are  not  statically 
modeled.  This  flexibility  can  cause  large  negative  twist  values  for  swept 
wings.  Figure  2.4-3  shows  that,  at  y/(b/2)  =  0.804  for  the  TF-8A  wing 
operating  at  M_^  =  0.99,  any  twist  of  the  flutter  model  (beyond  that  built-in 
to  the  pressure  model)  of  greater  than  a  -2  degrees  will  put  the  flutter 
model  into  an  apparent  stall  condition.  At  other  sections,  which  do  not 
stall  quite  so  easily,  the  slope  of  the  lift  curve  at  a  negative  twist  angle 
may  change  drastically  from  its  untwisted  position.  Reynolds  Number  plays  a 
part  here  since  the  nonlinear  shape  of  the  curve  is  a  function  of  it.  Thus, 
for  insufficient  Reynolds  Number,  the  lift  and  moment  curve  slopes  are  net 
to  be  trusted  except  in  the  linear  range.  Figure  2.4-6  presents  the  effect 


O  Experimental  Data  (“gj^p  = 
A  Revised  Cp  Distribution 
□  Cp  Using  Original  Airfoil 


5urf*c* 


Original  Airfoil 
Revised  Airfoil 


Fiqure  2.4-5  The  Effects  of  Boundary  Layer  Displacement  Thickness  on  a 
Section  of  the  TF-8A  Wing  {y/(b/2)  =  0.933,  M^q  =  0.90, 
Mon  =  0.72,  (X,  =  -3') 


of  Reynolds  Number  on  the  lift  curve  slope  for  positive  angles  of  attack. 
There  are  noticeable  changes  In  c^^  for  small  changes  In  Reynolds  Number. 
Thus,  If  the  Reynolds  Number  for  tWi’^pressure  model  test  1s  Insufficient 
(It  should  be  approximately  10^)  then  the  results  for  lift  and  moment  curve 
slopes  may  be  Inaccurate  for  angles  of  attack  In  the  nonlinear  range. 

For  the  case  of  the  TF-8A  wing,  the  pressure  test  was  run  at  a  Reynolds 
Number  of  1.96  x  10^  (based  on  the  chord  at  y/(b/2)  =  0.75)  where  as  the 
flutter  model  was  run  at  a  Reynolds  Number  ranging  from  3.0  to  9.0  x  10®. 

Thus,  It  must  be  concluded  that  the  nonlinear  portions  of  the  lift  and 
moment  curves,  as  obtained  from  the  pressure  model,  probably  do  not  apply 
to  the  f 1  utter  model . 

The  shock  location  also  has  a  large  effect  on  the  unsteady  and  steady 
weight  factors.  Thus,  the  twist  and  angle  of  attack  of  the  flutter  model, 
or  aircraft,  must  be  known  at  the  flutter  speed  If  these  weight  factors  are 
to  be  calculated  properly. 
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2.5  TRANSONIC  TWO-DIMENSIONAL  LIFTING  SURFACE  THEORY 


Three  approaches  to  the  unsteady  transonic  flow  problem  are:  (1)  direct 
numerical  solution  of  the  time  dependent  differential  equations  step-by-step 
in  time  using  finite  difference  techniques  (see  References  5  and  6);  (2) 
direct  numerical  solution  of  the  oscillatory  transonic  differential  equations 
using  finite  difference  or  finite  element  techniques  (see  References  7,  8, 
and  9);  and,  (3)  solution  of  the  transonic  problem  using  lifting  surface 
theory  (see  References  3  and  10). 

The  first  approach,  which  is  currently  operational  in  two-dimensional 
analyses,  accounts  for  transient  shock  motion.  The  cost  of  this  type  of 
approach  is  prohibitive  in  production  flutter  and  PSD  gust  analysis,  at  least 
with  current  generation  computers. 

The  second  approach  shows  promise  even  though  shock  motions  are  currently 
not  accounted  for  and  even  though  it  is  restricted  to  low  reduced  frequencies. 
The  dependence  of  the  oscillatory  transonic  solution  on  the  steady  state 
solution  multiplies  the  number  of , cases  that  must  be  done  in  a  routine  flutter 
investigation,  thus,  also  multiplying  the  cost. 

Classical  subsonic  lifting  surface  theory  has  gone  a  long  way  in  pre¬ 
dicting  flutter  speeds  even  in  the  transonic  region,  especially  when  steady 
weight  factors  are  applied.  It  seems  logical  that  modifications  of  lifting 
surface  theory,  to  account  for  transonic  effects,  is  a  profitable  avenue  to 
investigate.  The  third  approach  deals  with  this  area.  The  approach  of 
Cunningham  (Reference  10)  is  encouraging.  This  approach  uses  a  local  appli¬ 
cation  of  classic  subsonic  and  supersonic  lifting  surface  theories  coupled 
with  an  appropriate  shock  motion  analysis. 

The  method  of  Reference  3  addresses  the  problem  of  modifying  the 
classic  lifting  surface  theories  themselves  to  account  for  transonic  effects. 
The  lifting  surface  theory  approach  to  the  problem  is  economically  the  most 
desirable  if  it  produces  sufficient  accuracy. 

In  this  section  of  the  report,  an  investigation  into  a  two-dimensional 
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transonic  lifting  surface  theory  is  presented  which  is  an  extension  of  the 
initial  work  of  Reference  3.  The  theory  is  split  up  into  a  part  which  is 
independent  of  shock  wave  notion  and  one  part  that  deals  strictly  with  whe 
effects  of  the  shock  wave  motion.  The  method  is  correlated  using  the  NLR 
7301  airfoil. 

2.5.1  Basic  Theory  Without  Shock  Wave  Motion  Effects 

In  this  section,  the  basic  differential  equation  for  the  potential, 
is  solved  using  Green's  theory  which  results  in  a  surface  singularity 
approach.  The  nonuniformity  of  the  flow  field  is  accounted  for  in  the  express¬ 
ion  for  the  point  singularities  used.  The  airfoil  boundary  conditions  are 
satisfied  by  adjusting  the  singularity  strength.  The  basic  singularity  used 
in  lifting  surface  theory  is  the  doublet;  however,  the  doublet  can  be  obtained 
from  an  even  more  basic  singularity,  i.e.,  the  source.  The  derivation  of 
the  expression  for  a  transonic  source  is  given  in  the  following  and  in 
Appendix  A. 

The  shock  boundary  conditions  are  dealt  with  in  Section  2.5.2. 

2. 5. 1.1  The  Transonic  Source 

Lifting  surface  theory  is  based  on  the  assumption  th?t  the  disturbances 
to  the  fluid  caused  by  a  surface  moving  through  it  are  small.  These  small 
disturbances  travel  through  the  fluid  at  the  speed  of  sound  as  acoustic 
waves.  Classic  lifting  surface  theory  makes  the  further  assumption  that  the 
medium  through  which  the  waves  travel  is  uniform.  It  is  felt  that  one  of  the 
major  transonic  effects  arises  due  to  the  fact  that  acoustic  waves  encounter 
nonuniform  flow  fields.  This  is  especially  important  for  forward  traveling 
waves  in  a  near  sonic  stream.  Small  changes  in  fluid  speed  can  change  the 
direction  of  motion  of  these  waves  thereby  changing  the  character  of  the 
flow.  A  derivation  of  an  expression  for  the  potential,  (}>,  due  to  a  source 
in  a  nonuniform  stream  is  presented  in  Appendix  A.  The  result  is: 

N  ^ 

i)(x,y,z,t)  =  ^  ° 

i=l  ^ 
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The  term  o  is  the  source  strength  which  is  a  function  of  time.  To  understand 
this  expression,  it  must  be  noted  that  a  source  in  a  fluid  stream  is  actually 
an  infinite  string  of  acoustic  pulses;  one  born  at  each  instant,  traveling 
downstream.  Figure  2.5-1  shows  a  source  in  a  unifotnn  flow  located  at 
(^Q.  Hq.  Cq)  with  an  illustration  of  the  trailing  acoustic  pulses.  Notice 
the  two  pulses  that  reach  the  receiving  point  {x,y,z)  at  the  particular  time 
in  question.  These  pulses  were  generated  previously  one  at  t  =  and  one 
at  t  =  ^2'  sending  or  source  point  (Cq. 

Once  an  acoustic  pulse  is  generated,  its  total  strength  remains  unchanged. 
Of  course,  the  local  wave  strength  on  the  pulse  perimeter  varies  with  time 
as  the  v/ave  expands  and  distributes  its  strength  along  an  ever-lengthening 
perimeter.  Thus,  the  total  strengths  of  the  waves  that  have  reached  the 
receiving  point  (x,y,z)  are  a  function  of  the  time  at  which  they  were 
generated  at  the  point  (c^,  5q).  The  strengths  of  the  waves  reaching 

(x,  y,  z)  are  the  strengths  of  the  source  at  (C^.  n^.  at  the  previous 
times  t  =  T.|  and  T2.  This  idea  is  written  as  o(t^)  and  oCt^). 

For  the  uniform  flow  supersonic  case,  these  times  t^,  or  phase  lags, 
are  expressed  as; 

'1  -  -7  (X  -  t„)  *  *  t  i  -  A 

a  3^  I 

'  Uniform  super- 
1  sonic  case 

'2  =■  -7  (N.  (X  -  *  t  <  =  2  ) 

For  the  uniform  flow  subsonic  case,  only  one  wave  passes  the  receiving  point, 
thus 

T,  =  (M  (x  -  c  )  -  R)  +  t  i  =  1  Uniform  Subsonic  Case 

CO 

where 

=  (x  -  +  6^r^ 

r^  =  (y  -  +  (z  - 
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a  =  Speed  of  sound 

In  general,  there  are  N  such  waves  passing  the  receiving  point  for  the 
nonuniform  flow  case.  For  this  case,  a  special  acoustic  pulse  generator 
is  required  to  determine  the  time  histories  of  the  shapes  and  locations  of 
the  generated  acoustic  pulses  in  a  nonuniform  flow.  Such  an  acoustic  pulse 
generator  is  described  in  subsection  2. 5. 1.4. 

In  the  expression  for  the  source  potential,  ij),  the  term  R^.  is  expressed 
in  terms  of  acoustic  pulse  parameters  as  follows: 


a 


where  Vj  is  the  total  velocity  of  a  point  on  the  acoustic  pulse  surface, 

*n  its  surface  normal  vector  in  the  direction  of  travel  (see  Figure  2.5-2),  a, 
an  equivalent  speed  of  sound  and  p  an  equivalent  acoustic  pulse  radius.  The 
acoustic  pulse  generator  is  used  to  determine  Vy  and  n.  The  equivalent  speed 


Figure  2.5-2  Pulse  Surface  Velocity  Vectors 
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of  sound  a  is  dotcrmined  iinpllcltly  from  the  followinQ  formula 
Vj  =  V  +  a  n 

A 

This  formula  simply  says  that  the  total  velocity,  Vj,  of  the  wave  front  is 
a  vector  combination  of  the  fluit  velocity  V  and  the  acoustic  velocity  of 
the  pulse,  a  n.  This  is  observed  graphically  in  Figure  2.5-2.  The  nonuniform 
flow  field  iT  must  be  known  from  a  suitable  steady  transonic  flow  solution; 
in  this  case  the  Jameson  Method  (see  section  2. 5. 1.4).  Since  the  total,  V.^., 
is  known  from  the  acoustic  pulse  generator  the  equivalent  speed  of  sound,  a, 
may  be  solved  for.  When  no  approximations  forfi  are  used  a  =  a  and  the  formula 
for  a  is  not  needed.  The  expression  for  a  is  kept,  however,  for  later  refer¬ 
ence  as  approximations  for  n  are  introduced, 
i  =  (V^  -  V)  •  n 

The  effective  local  Mach  Number,  M,  is  then 

i .  4l=  ^-I-, 

a  Vj  •  n  -  V  •  n 

The  inverse  of  the  term  p  is  a  measure  of  the  local  wave  strength  along 
the  wave  parameter.  For  the  uniform  flow  case  this  term  reduces  to  the 
acoustic  pulse  radius,  p,  (see  Figure  2.5-1)  where: 

p(t)  =  a(t  -  T^) 

and  where  a  is  the  speed  of  sound  in  a  uniform  flow.  A  factor  f  is  introduced 
to  account  for  the  fact  that  the  wave  strength,  in  a  nonuniform  flow,  is  not 
exactly  proportional  to  the  inverse  of  the  acoustic  pulse  radius. 

p  =  fp  -  fa  (t  -  T^) 

For  uniform,  flow.  Appendix  A  shows  that  simply  becomes  the  hyperbolic 

radius  for  both  the  subsonic  and  supersonic  cases.  That  is, 

R,  ^  R  =  (x  -  C  r^  i  =  1,2  Uniform  Flow  Case 
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If,  further,  the  oscillatory  case  is  considered  for  uniform  flow,  i.e., 
o(t)  =  a  exp(iuT),  then  the  expression  for  the  potential  becomes: 


4 


1  £ 
471  R 


For  the  supersonic  case  N  -  2,  and  using  the  expressions  for  t2  gives 

"  2^  cos|^jexp[1uM^(x  -  f,p)  /  ae^^]  e^*^^ 


For  the  subsonic  case  N  =  1 ,  and  using  the  expression  for  gives 
<^(t)  =  ~  exp  [iu)M^  (x  -  /  aej]  exp  [-ii.R/aB^^]  e’*"^ 

These  are  the  classic  results  for  the  potential  of  the  supersonic  and  subsonic 
source. 


2. 5. 1.2  Equivalent  Sending  Point  Location 


The  acoustic  pulse  generator,  to  be  described  in  subsection  2. 5. 1.4, 
provides  a  means  of  calculating  the  time  history  of  an  acoustic  pulse  emitted 
from  a  sending  point  C^)  moving  in  a  nonuniform  steady  transonic 

flow  field.  When  the  pulse  has  arrived  at  the  receiving  point,  it  has  dis- 
tored  from  the  classic  circular  shape  and  its  wave  strength  has  deviated 
from  the  classic  value  of  a/p  where  p  is  the  pulse  radius.  However,  locally 
at  the  receiving  point  (x,  y,  z)  an  equivalent  circular  pulse  can  be  matched 
to  the  actual  pulse  (see  Figure  2.5-3). 


The  quantities  to  be  matched  are  the  terms  t-  and  R.  appearing  in  the 
expression  for  the  source  in  a  nonuniform  flow  potential: 


N 


o(t^) 

% 

47tR^ 


that  is,  the  terms  t.  and  R^.,  for  the  equivalent  circular  pulse,  are  to  be 
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matched  to  the  actual  pulse  at  the  receiving  point  (x,  y,  z).  Once  the 
match  is  made  an  equivalent  sending  point  location  (5^.  n^.  flow  field 
speed  V,  speed  of  sound  a  and  pulse  radius  p  can  be  found.  The  bar  over  the 
terms  indicates  equivalent  terms  and  not  actual  ones.  The  derivation  of 
these  equivalent  terms  is  given  as  follows. 


-T  } - ^ 

(Acoustic 
Pulse 
Center) 


Act'ial  Pulse 


Matched  Circular  Pulse 

Figure  2.5-3  Matching  a  Circular  Pulse  to  the  Actual  Pulse  at  the 
Receiving  Point  (x.y.z) 


From  Figure  2.5-3,  the  following  vector  identity  is  easily  observed: 
p  =  (x  -  1  +  r  =  V  (t  -  t)  +  pn 

or  by  components: 

(x  -  Tq)  =  V  (t  -  x)  +  pn^^ 
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and 


r  =  pn^ 

where  =  n  •  r/r 

Using  these  quantities  in  the  expression  for  the  equivalent  hyperbolic 
radius,  R,  gives: 

=  (V  (t  -  t)  +  pn^)^  +  n^^ 

~2  (9  -  -9  9  ) 

=  I  t  2n^  V5  +  V^n/  +  L  j 

where 


L  =  n^^  (V^  -  M^a^) 


Now  if  we  require  that 


M  = 


then  L  =  0  and  R  becomes 


R  =  p 


Vn. 


1  + 


It  can  be  shown  that  this  expression  reduces  to  the  expression  for  R/f  as 
follows: 


R 


•  n 


fR 


when  i^t  is  assumed  that  the  equivalent  velocity  is  parallel  to  the  x-axis, 
i.e.,V  =  y  i  and  when  we  take  note  of  the  expressions 


=  V  +  a  n 


•v 

P  = 


fp 
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which  were  given  in  subsection  2.5.1.  Thus,  the  equivalent  hyperbolic  radius 
is  fR  and  the  expression  for  the  potential  becomes: 


e 


N 


E 

i=l 


a(t^) 


The  quantities  V.  p.  f  and  a  are  yet  to  be  determined.  First,  assume 
that  V  is  the  local  flow  field  velocity,  at  the  receiving  point  (x,  y,  z). 
Secondly,  the  pulse  radius  p  can  be  selected  by  fitting  the  best  radius  to 
the  actual  pulse  or  simply  by  letting  it  be  a(t  -  t)  where  a  is  the  equiva¬ 
lent  speed  of  sound.  For  simplicity,  the  latter  is  assumed.  The  quantity  f 
is  then  obtained  by  noting  the  actual  strength  of  the  wave  and  comparing 
it  to  a(t  -  t),  i .e. , 

f  =  — ^ 
a(t  -  t) 


I _ _ _ _ _ ^ 

Figure  2.5-4  Pulse  Surface  at  Time  *  t  and  t+At  and  Associated 
Vector  Quantities 


Pulse  Surface 
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Thirdly,  the  equivalent  local  speed  of  sound  can  be  obtained  by  refering 
to  Figure  2.5-4.  This  figure  shows  the  position  of  a  wave  at  two  time  inter¬ 
vals  At  apart.  Inspection  of  the  figure  shows  that: 

Vt  =  n  a  +  V  1  =  ^  n 

The  term  An/At  can  be  obtained  from  the  acoustic  pulse  generator  in  a 
numerical  manner  from  the  pulse  shape  at  two  successive  times.  Lastly, 
knowing  the  values  of  V  and  n  allows  a  to  be  solved  for  as  follows: 

a-  = 

All  terms  in  the  expression  for  4>g  now  known  and  it  can  be  used  in 

the  form  presented.  However,  the  classic  potential  is  usually  not  used  in 

this  form.  Usually  the  more  classical  one  given  at  the  end  of  subsection 

2. 5. 1.1  is  used.  In  these  expressions,  the  terms  t.j  and  ^2  are  given  in 

terms  of  (x  -  M  and  R  as  given  in  Appendix  A.  According  to  Figure 

o 

2.5.3, 


+  [(x  -  5q)  +  V  (t  -  t)]^  =  S^t  - 


solving  this  expression  in  terms  of  (t  -  i)  gives 

(x  -  L)  V/i  i  R 


+  (t  -  -c)  = 


i  (1  -  M^) 


Since  M  =  V/a  then 


and 


R^  =  (x  -  f 


+  u(t  -  t) 


2M 

eg 


L(x  -  ?g)  M  IR  ] 
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where  =  <*>c/2U^ 


Thus,  the  classic  expresslor  for  4t  -  t)  Is  nidified  by  simply  replacing 
M  with  ii.  and  (X  -  5.)  and  r  with  (x  -  t„)  and  r  and  dividing  by  WU  . 
The  classic  oscillatory  subsonic  case  has  one  root,  i.e.,  1  =  1  an  e 
result,  taken  from  the  end  of  subsection  2. 5. 1.1,  Is. 


The 


({i  =  e 


^  exp  [  (M„  (x  -  -  R)  ^  3 


If  we  introduce  the  equivalent  qualities*  i.e., 

(X  -  Cq)  (x  -  Cq)  *  V(t  -  t)  +  P 

r  r  = 

M  R  =  V/a 


where 

a  =  (i^-V1  ••!?) 

p  =  (t  -  t)! 

and  where  5  Is  the  local  flow  field  velocity  at  the  receiving  point  and 
(t  -  t)  Is  the  time  required  for  the  acoustic  pulse  to  move  from  the 
sending  point  to  the  receiving  point;  then  the  equivalent  potential 
for  a  source  in  a  nonuni  form  flow  is 

=  eia)t  _L_  exp  [  (M  (X  -  ?_)  -  R)  /  5  V)  ] 

e  4TrRf 


where 

R^  =  (x-  f 

Here  the  quantities  (t  -  x).  An/At  andli  are  obtained  from  the  acoustic 
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pulse  generator  and,  V,  from  the  steady  transonic  flow  generator. 
2. 5. 1.3  The  Two  Dimensional  Transonic  Source 


The  three-dimensional  source  can  be  written  as: 


II 

-e- 

7 

♦e  ' 

^e  = 

5  1  M  I{x  -  i^) 

47rf  ® 

-iAR 

e 

R 

where  A 

=  A/V/U„ 

where  A 

2 

c 

The  two-dimensional  source  can  be  obtained  from  the  three-dimensional  one 
by  integr 
infinity. 

♦2D  =  / 


by  integrating  the  sending  point  in  the  ng-direction  from  plus  to  minus 


♦e  '‘”0 


If  we  assume  f  is  independent  of  n  then 

-  iMA(x-i  ) 


7D  4Trf 


°  f  e'^ 

J  D 


dn. 


but 


/ 


-i  AR  . 
e _  dn. 


=  2 


/ 


e^^*^  dR 

Kjd  0 


=  r  ft 

-a  Jr- 


-iAR 


dR 


'2D 


{W- 


2D 


where 


■  (X  -  1/  -  6^  (  .  -  )' 


< 


I 
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Introducing  0  =  R/^2D 


(See  Appendix  C  for  expression  for  supersonic  flow) 

The  only  difference  between  this  source  and  the  classic  uniform  flow  source 
are  the  quantities  f,  V/U..  fi  and  x  -  and  z  - 

The  assumption  Is  now  made  that  the  source  potential, 
erentlated  and  Integrated  to  produce  the  velocity  field  due  to  a  pressure 
doublet.  This  Is  necessary  because  the  pressure  doublet  Is  the  basic  singul¬ 
arity  used  In  lifting  surface  theory. 

This  assumption  Is  probably  accurate  for  d1fferent1at<on;  however.  Inte¬ 
gration  requires  consideration  of  a  sending  point  that  varies  from  the  origin 
to  downstream  Infinity.  The  first  step  In  obtaining  the  downwash  due  to  a 
pressure  doublet  Is  to  obtain  the  downwash  due  to  a  velocity  doublet  by  a 
double  differentiation.  I.e..  The  second  step  Is  to  Integrate  the 

effects  of  a  distribution  of  velocity  doublets  along  a  line  originating  at 
(C  .  c  )  and  passing  to  downstream  Infinity  parallel  to  the  x-axIs.  The  most 
accurate  way  to  handle  this  Is  to  obtain  a  new  (t  -  t).  dn/at.  etc.  for  points 
along  the  Integration  path  and  perform  the  Integration  numerically.  For  now 
however,  the  values  of  (t  -  x).  on/dt.  etc.  will  be  obtained  for  the  origin  only 

U  c  ).  In  essence,  this  assumption  states  that  the  flow  is  nonuniform  from 

the'seSdIng  point  (L.  I  j  to  the  receiving  point  (x.  z)  but  that  the  flow  is 
uniform  from  a  point  In  the  wake  of  the  sending  point  up  to  the  sending  point 
itself.  Figure  2.5-5  gives  a  graphical  illustration. 

With  these  assumptions  one  simply  introduces  the  equivalent  geometry. 

(x  -  Eq)  and  (z  -  Eq),  along  with  the  equivalent  Mach  Number,  M,  into  the 
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Thus  the  two-dimensional  equivalent  source  is 

«(x  - -g 
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Non-uniform  flow  region 
Figure  2.5-5 


Uniform  flow  region 


Result  of  Assuming  that  the  Equivalent  Geometry  (^  ,  f  ) 
Applies  Directly  to  a  Point  Pressure  Doublet.  ^ 


expression  for  the  point  pressure  doublet  (used  in  lifting  surface  theory) 
and  then  proceeds  to  solve  the  resulting  equations  in  tfie  usual  manner. 

2. 5. 1.4  Acoustic  Pulse  Generator 

In  the  previous  subsections  expressions  were  derived  for  equivalent 
distances  and  Mach  Number's  based  on  the  properties  of  acoustic  pulses  traveling 
in  nonuniform  flow  fields.  It  was  assumed  that  the  following  data  were  obtain¬ 
able  from  these  acoustic  pulses:  (1)  (t  -  t),  the  time  taken  for  a  pulse  to 
travel  from  its  origin  to  a  receiving  point  (x,  z);  (2)  n,  the  normal 

vector  of  the  pulse  surface  at  the  receiving  point;  and  (3)  f,  the  dialation 
of  the  pulse  over  and  above  the  usual  time  dialation.  These  quantities  are 
known  if  the  time  history  of  the  pulse  is  known.  In  this  subsection,  a  method 
of  computing  the  pulse  time  history  in  a  nonuniform  flow  field  will  be 
presented. 

Basically,  a  pulse  surface  is  made  up  of  a  series  of  points  each  of  which 
has  its  own  time  history  or  trajectory.  The  velocity  of  these  surface  points 
depends  on  the  flow  field  velocity,  the  speed  of  sound,  and  the  normal  direction 
of  the  surface  in  the  direction  of  motion. 
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Figure  2.5-6  Quantities  Associated  With  the  Trajectories  of 
Points  Lying  on  a  Pulse  Surface 


Specifically,  Figure  2.5-6  shows  that  the  total  velocity  Vj  of  a  pc^nt 
in  the  acoustic  pulse  surface  is  the  vector  addition  of  the  flow  field,  V, 
and  the  pulse  front  velocity,  which  is  the  speed  of  sound  in  the  direction  of 
the  surface  normal,  i.e.,  a"?.  Thus  the  total  is: 

Vj  =  V  +  na 
or  in  the  normalized  form. 


where  M  is  the  free  stream  Mach  Number. 
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The  normal  vector,  n,  at  a  pulse  surface  point  is  obtained  by  averaging 
the  normal  vectors  of  the  elements  on  either  side  of  the  [joint  as  Figure 
2.5-6  illustrates.  Thus, 

n  =  (n^  +1^2)  !  2 

More  sophisticated  curve  fitting  procedures  were  tried  but  they  were  found  to 
be  time  consuming  and  sometimes  unstable.  This  simple  procedure  has  worked 
best. 


Figure  2.5-7  Grid  Used  for  The  Interpolation  of  The  Transonic 
Flow  Velocity  Field 


The  velocity  V(  x,  z)  was  obtained  by  an  interpolation/extrapolation 
procedure  based  on  a  set  of  velocities  obtained  from  a  modified  Jameson 
steady  two-dimensional  transonic  computer  program  (Reference  11). 
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The  onset  velocity  field  V  is  broken  up  into  a  uniform  component  U^i  and 
a  nonuniform  increment.  aV.  The  interpolation/extrapolation  is  performed  on 
the  nonuniform  increment.  aV.  The  vector  notation  will  be  dropped  for  conven 
ience.  Also,  with  this  program,  it  was  assumed  that  AV^k  was  small  enough  to 
ignore.  The  interpolation  in  the  x-direction  is  linear  as  follows: 

aV(x.  z)  =  [AVj^.i(z)  -  AVj{z)]  mx  +  aV.(z) 

where 


and  where  XR.  are  the  x-locations  of  AVj(z).  i.e..  the  values  to  be  inter¬ 
polated  from.  Figure  2.5-7  Illustrates  the  grid  of  points.  The  tilda  over 
the  z  indicates  a  transformed  plane  where  the  airfoil  surface  transforms  into 
the  z  =  0  plane.  The  upper  and  lower  surfaces  are  calculated  separately  and 
treated  independently.  The  transformation  required  to  render  z  =  0  on  the 
surface  is  assumed  small  enough  to  ignore  subsequent  to  the  transformation. 


The  determination  of  AVj(z)  is  either  found  by  linear  interpolation,  for 
z  <  ZR(max-l)  fitting  for  z  >  ZR(max-l)*  expression  for 

interpolation  is: 

AVj(z)  =  [AV.lzR^^.^)  -  AVj(zR^)]  mz.  +  aV^IzR.)  z  lZR(n,ax-l) 


where 


mz, 


z  -  zR^ 


For  the  case  z  >  ZR 


(max-1) 


the  formula  to  use  is; 


iVjz)  =  iV,(ZR  ) 


b 
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This  formula  has  the  property  of  approaching  zero  like  1/z  as  z  approaches 
infinity.  This  formula  is  fitted  through  the  last  two  points;  and 

ZRmax_i •  Figure  2.5-7  also  illustrates  this  curve. 

The  values  of  aVj.(ZR^)  are  obtained  from  a  modified  version  of  the  steady 
flow  two-dimensional  transonic  method  of  Jameson  (Reference  11).  The 
method  had  to  be  modified  since  only  surface  values  of  velocity,  pressure  and 
Mach  Number  were  computed.  Appendix  B  presents  the  method  used  in  modifying 
the  Jameson  procedure  to  produce  the  flow  field  velocities  off  of  the  body 
surface.  Basically,  what  was  available  was  the  potential  at  various  grid 
points  in  a  plane  transformed  so  that  the  space  around  the  airfoil  was  mapped 
into  a  rectangle,  one  unit  high  and  two  units  wide.  The  procedure  involves 
multiple  mappings  and  Appendix  B  presents  the  properties  of  these  and  their 
use  in  obtaining  the  flow  field  velocities.  Figure  2.5-8  presents  an  example 
of  the  flow  field  incremental  velocity,  V  -  U^,  as  calculated  by  the  modified 
Jameson  procedure  for  the  NLR  7301  airfoil  operating  at  =  0.745  and 
a  =  -0.20. 


The  acoustic  pulse  generator  is  initiated  by  prescribing  an  initial  set 
of  points,  given  on  the  arc  of  a  very  small  circle.  As  time  passes,  the  pulse 
radius,  and  thus,  the  distance  between  adjacent  points,  becomes  large.  An 
automatic  procedure  for  filling  points  in  between  these  initial  points  has 
been  implemented.  Thus,  the  character  and  accuracy  of  the  pulse  surface  is 
maintained  even  if  the  surface  is  dialated  excessively  at  various  points. 
Points  are  filled  in  at  a  rate  such  that  the  distance,  D,  between  adjacent 
points  does  not  exceed: 

0  ^  A(t  -  T)a 

where  A  is  a  prescribed  parameter. 

Example  calculations  of  the  acoustic  pulse  generator  are  shown  in 
Figure  2.5-9  through  2.5-11.  Figure  2.5-9  presents  the  time  history  of  the 
surface  of  an  acoustic  pulse  in  a  nonuniform  flow  whose  origin  is  at  the  55% 
point  on  the  NLR  7301  airfoil  operating  at  =  0.7  and  a  =  1.75°.  Below 
the  x-axis  is  presented,  for  reference,  the  time  history  of  the  same  pulse  in 
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Figure  2.5-8  Non-uniforw  Flow  Field  Increment  Over  the  NLR  7301  Airfoil 
As  Calculated  by  the  Jameson  Transonic  Flow  Method 
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Figure  2.5-9  Time  History  of  An  Acoustic  Pulse  h'hose  Origin  Is  At 
x/c  ■  0.55  For  The  Nonunlfortn  Flow  (above  axis)  and 
Uniform  Flow  (below  axis)  Cases.  Nonuniform  Flow  Is 
That  for  The  NLR  7301  Airfoil  (Mco-  0.7,  o-  1.75°)  oi 
the  Upper  Surface 


Figure  2.5-10(a) Time  History  of  an  Acoustic  Pulse  Hhose  Origin  is  At 
^  x/c  ■  0.70  For  The  Non-Uniform  Flow  (above  axis)  And 

Uniform  Flow  (below  axis).  Cases.  Non-Uniform  Flow  Is 
That  For  the  NLR  7301  Airfoil  (Mco“  0.745,  o  «  -0.2) 
On  The  Upper  Surface. 


s 


Uniform  Alrfol 
Flow 


(b)  Time  History  of  an  Acoustic  Pulse  Whose  Origin  is 
at  x/c  =  0.4  for  the  Non-Uniform  Flow  (Above  Axis) 
and  Uniform  Flow  (below  axis)  Case.  Nonuniform 
Flow  is  that  for  the  NLR  7301  Airfoil  (Ma,=  0.745, 
a  =  -O.Z”)  on  the  Upper  Surface. 
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Time  History  of  An  Acoustic  Pulse  Whose  Origin  Is  At  x/c 
For  The  f!on-Un1forr  Flow  (above  axis)  and  Uniform  Flow  (I 
axis)  Cases.  Non-Uniform  Flow  Is  That  for  The  TF-8A  Win 
Section  at  y/(b/2)  -  0.657  (M20  -  0.711  «2d  “  o"  ^ 
Upper  Surface. 


Figure  2.5-11 


in  a  uniform  flow.  The  time  history  depicted  below  the  x-axis  is  the  type 
assumed  for  the  classic  lifting  surface  theory,  where  as,  that  shown  above  the 
x-axis  is  the  actual  one. 

The  nonuniform  flow  field  found  on  the  upper  surface  of  the  NLR  7301 
airfoil  operating  as  indicated  possesses  a  supersonic  zone  ahead  of  the  40% 
point  which  is  terminated  with  a  shock  wave  that  extends  to  a  height  of  0.2c. 

A  small  secondary  supersonic  zone  occurs  aft  of  the  shock  also.  The  acoustic 
pulse  surface  moves  forward  very  slowly  at  first  slowing  down  even  further  as 
it  approaches  the  shock  wave  location  (40%  point).  The  upper  part  of  the 
acoustic  pulse  moves  forward  of  the  lower  part  causing  a  rotation  in  the 
inclination  of  the  wave  surface.  As  the  wave  rotates  counter  clockwise  it 
starts  moving  into  the  supersonic  zone;  the  greater  the  rotation  the  faster 
it  moves. 

Figures  2.5-10{a)  and  (b)  presents  a  similar  calculation  for  the  NLR 
7301  airfoil  operating  at  the  design  "shock  free"  point  of  =  0.745  and 

a  ■  -0.20.  Even  through  there  is  no  shock  the  supersonic  zone  still  acts 
like  a  barrier  to  the  advancing  acoustic  wave.  The  wave  shapes  are  quali¬ 
tatively  similar  to  those  for  the  shock  wavs  condition  except  that  the  pulse 
passes  into  the  supersonic  zone  a  little  more  smoothly. 

Figure  2.5-11  presents  a  third  example  calculation.  This  involves  the 
TF-8A  wing. 

2. 5. 1.5  Application  to  Lifting  Surface  Theory 

The  transonic  approach  taken  in  this  report  consists  of  finding  equivalent 
geometry  and  Mach  Number  that  exist  between  a  sending  anc!  receiving  points  and 
then  placing  these  equivalent  quantities  into  the  classic  theory  for  the 
desired  transonic  result.  Thus  a  classic  method  is  required  that  will  cover 
both  the  subsonic  and  supersonic  cases.  The  Doublet  Lattice  Method  is  subsonic 
only  and  thus,  a  supersonic  extension  is  needed.  Such  an  extension  for  the 
two-dimensional  case  is  derived  in  Appendix  C.  As  a  test  case,  the  supersonic 
Doublet  Lattice  Method  was  compared  to  the  analytic  solution  for  an  airfoil 
pitching  about  its  mid-chord  at  M  =  1.2  and  k^  =  0.99.  The  result  shown  in 
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Figure  2.5-12  is  very  satisfactory. 


Figure  2.5-12  Comparison  of  the  Present  Method  With  the  Analytic 
Solution  for  The  Oscillatory  Supersonic  Case. 
Airfoil  Pitching  About  Mid-Chord  at  k  •  0.99  and 
M«,-  1.2 


Cunningham  (Reference  10)  has  shown  that  the  downwash  boundary 
conditions  should  be  modified  for  the  transonic  case. 


Or,  the  modified  downwash,  w,  is  the  unusual  downwash,  w,  times  the  ratio  of 
local  surface  Mach  Number,  M,  to  the  free  stream  Mach  Number,  M^.  Cunningham 
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goes  on  to  modify  the  reduced  frequency  in  a  similar  manner. 

In  the  method  developed  for  finding  an  equivalent  geometry  (subsection 
2. 5. 1.2)  the  possibility  exists  that  a  planar  case  (airfoil  in  z  =  0  plane) 
may  become  nonplanar  because  of  the  acoustic  pulse  matching  procedure.  That 
is,  the  matched  circular  pulse  may  have  its  origin  (5^.  1^)  off  the  airfoil 
surface  even  though  the  actual  origin  (t^i  c^)  is  on  the  surface.  This  is 
especially  true  for  acoustic  pulse  waves  traveling  upstream  over  the  top  of 
supersonic  zones.  As  Figures  2.5-9  through  2.5-11  show,  these  pulse  waves 
become  inclined  to  the  airfoil  surface.  Such  an  inclination  for  matched 
circular  pulse  waves  indicates  a  pulse  origin  that  lies  off  of  the  airfoil 
surface. 

At  this  stage  of  development  of  the  method,  such  a  complication  of  the 
calculation  is  unjustified.  Thus,  an  assumption  is  introduced  at  this  point. 
The  fitted  circular  pulse  origin  is  assumed  to  lie  on  the  airfoil  surface. 
This  means  that  the  wave  normal  vector  is  not  matched  but  the  time  of  arrival 
(t  -  t)  is  matched.  The  equivalent  values  of  geometry  and  flow  parameters 
then  become: 

(*  -  Ij)  Mt  -  i)  + 

where 

=  sign  (*  -  5„) 
r  -►  0 
M  =  V/a 

where 

a 

p  =  (t  -  t)  a 

Where  ax  is  the  incremental  distance  along  the  x-axis  that  the  pulse  has 
traveled  during  time  At.  With  these  quantities  (x  -  l^)  and  M,  can  be 
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written  as: 


(x  -  l^)  -  (t  -  t)  (1  +  6^/M)V/U„ 

M  =  1  /  ^  -  1 
AtV 

The  values  of  (t  -  t)  and  M  are  introduced  into  the  subsonic-supersonic 
two-dimensional  Doublet  Lattice  Method  to  account  for  transonic  effects.  To 
save  computation  and  input  time  and  expense,  the  values  of  t  -  t  and  M  are  not 
input 'for  every  sending/receiving  point  combination  (which  for  10  elements 
would  be  100  values)  but  are  input  at  a  selected  set  of  points.  The 
values  needed  are  then  obtained  by  a  double  linear  interpolation  scheme. 

Figure  2.5-13  presents  an  example  of  a  typical  time-of-arrival  plot.  Speci¬ 
fically,  the  time  of  arrival  of  an  acoustic  pulse  (whose  origin  is  located  at 
the  70%  point)  at  various  locations  on  the  airfoil  surface  is  presented. 

Notice  that  the  time  of  arrival  of  the  pulse  is  large  upstream  of  the  pulse 
origin  and  small  downstream  of  it.  This  is  because  upstream  moving  acoustic 
pulses  are  moving  against  the  flow  while  downstream  moving  pulses  are  moving 
with  the  flow.  The  result  for  the  uniform  flow  case  is  also  shown.  Even 
though  the  flow  at  the  airfoil  surface  is  supersonic  on  the  forward  portion 
of  the  airfoil,  the  acoustic  pulse  still  moves  upstream  (see  Figures  2.5-9 
through  11),  even  if  only  very  slowly. 

Notice  that  the  curve  of  (t  -  t)  vs  x/c  has  a  discontinuous  slope  at  the 
sending  point.  A  special  interpolation  procedure  is  used  because  of  this 
change  in  slope.  This  special  procedure  requires  that  the  curves  of  (t  -  t) 
vs  x/c  be  known  not  only  over  the  airfoil  surface  but  a  short  distance  up  and 
down  stream  of  it.  For  instance,  in  Figure  2.5-14,  if  the  desired  sending 
point,  XSI^,  lies  between  two  available  sending  points,  XSl  and  XS2,  then  an 
interpolated  curve  (shown  dashed)  must  be  constructed.  The  right  half  of  the 
interpolated  curve  is  obtained  from  the  right  halves  of  the  known  curves 
conversely  the  left  half  of  the  interpolated  curve  is  obtained  from  the  left 
halves  of  the  known  curves.  In  this  example  it  is  easily  seen  that  the  left 
half  of  the  known  curve  eminating  from  XSl  must  extend  forward  of  the  leading 
edge  so  that  an  interpolated  left  half  can  be  constructed.  The  converse  is 
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true  for  the  trailing  edge  when  dealing  with  the  right  half  of  an  interpolated 
curve  whose  origin  lies  between  XS4  and  XS5. 

The  local  effective  Mach  Number,  M,  is  also  a  discontinuous  function 
and  requires  special  treatment.  It  has  been  found  that  it  is  a  function 
mostly  of  the  receiving  point  location  and  whether  it  is  up  or  down  stream 
of  the  sending  point.  Figures  2.15(a)  and  (b)  present  the  local  effective 
Mach  Number  plotted  versus  the  receiving  point  location  for  various  sending 
point  locations  for  two  conditions.  The  following  conclusions  can  be  drawn 
from  these  figures.  First,  for  receiving  points  downstream  of  the  sending 
point,  the  local  effective  Mach  Number  is  close  to  the  local  surface  Mach 
Number.  Second,  for  receiving  points  downstream  of  the  supersonic  zone 
the  local  effective  Mach  Number  is  also  close  to  the  local  surface  Mach  Number 
independent  of  the  location  of  the  sending  point.  Third,  for  receiving 
points  upstream  of  sending  points  (in  a  supersonic  zone)  the  local  effective 
Mach  Number  falls  on  one  curve  independent  of  sending  point  location  and  the 
Mach  Numbers  are  subsonic.  The  reason  why  the  effective  Mach  Numbers  are 
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subsonic  for  this  condition  is  that  the  acoustic  pulse  wave  fronts  manage  to 
propagate  up  into  this  supersonic  zone  just  as  if  the  Mach  Numbers  were  sub¬ 
sonic.  Figures  2.5-9  through  11  illustrate  this  mechanism. 

Other  assumptions  have  been  made  to  facilitate  the  incorporation  of  the 
effective  geometry  and  Mach  Number  into  the  classic  subsonic  and  supersonic 
methods.  For  instance,  the  application  of  this  technique  to  pressure  doublets, 
which  are  used  in  the  Doublet  Lattice  Method,  requires  just  such  an  assumption. 
A  discussion  of  this  assumption  is  presented  in  subsection  2. 5. 1.3  but  it 
bears  repeating  here. 

Pressure  doublets  are  made  up  of  a  distribution  of  velocity  point 
doublets  along  a  line  starting  at  the  pressure  point  and  passing  to 

downstream  infinity  in  the  x-direction.  The  effective  geometry  and  Mach 
Number  should,  strictly  speaking,  be  applied  separately  to  each  of  these 
points  in  the  ‘wake'  of  the  pressure  doublet.  However,  to  simplify  this  pro¬ 
cedure,  it  is  assumed  that  the  effective  geometry  and  Mach  Number  for  the 
pressure  doublet  (with  wake)  are  those  based  on  the  leading  edge  of  the 
pressure  doublet,  i.e.,  the  sending  point  (Cqi  t^).  This  means  that  pulses 
that  elinate  from  the  wake  of  the  pressure  doublet  (aft  of  the  po^nt  (5^,  c^)) 
see  a  uniform  flow  between  them  and  the  point  (5^,  and  a  nonuniform  flow 
from  there  to  the  receiving  point  (x,  z).  The  effects  of  this  assumption  are 
not  easily  observed;  however,  in  general  they  should  not  be  too  large. 

Future  refinements  may  involve  the  elimination  of  this  assumption. 

Another  assumption  involves  the  basic  nature  of  the  transonic  pulse. 

The  basic  theory  of  transonic  nonuniform  flow  (subsection  2. 5. 1.1)  shows  that 
the  effect  of  an  acoustic  source  at  a  receiving  point  (x,  z)  depends  on  how 
many  times  the  pulse  passes  over  that  point.  The  summation  over  j  from  1  to 
N  in  the  equation  for  ^  in  subsection  2.5.11  and  Appendix  A  indicates  this 
fact.  In  classic  subsonic  flow,  the  acoustic  pulse  passes  over  a  given  point 
only  once  (see  Figure  2.5-16a).  The  upstream  facing  part  of  the  pulse  moves 
slowly  upstream  while  the  downstream  part  moves  rapidly  in  the  downstream 
direction.  In  supersonic  flow,  such  is  not  the  case.  The  pulse  wave  either 
passes  over  the  point  twice  (for  a  point  lying  in  the  aft  Mach  cone  of  the 
source)  or  not  at  all,  for  points  lying  outside  of  the  Mach  cone.  The  down- 
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Figure  2.5-16  Pulse  Wave  Patterns  for  Various  Flow  Conditions 
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stream  facing  part  of  the  wave  passes  first  and  then  the  upstream  facing 
part,  which  is  being  swept  downstream  (see  Figure  2.5-16b). 


For  transonic  flow,  the  qualitative  characteristics  of  pulse  wave  passage 
can  be  more  complicated  than  just  described  and  can  also  act  subsonically  in 
a  supersonic  zone  and  supersonically  in  a  subsonic  zone.  For  instance,  it 
it  has  been  shown  that  pulses  can  move  upstream  up  and  over  a  shock  wave  and 
down  into  the  supersonic  zone  (see  Figures  2.5-9  through  11  and  Figure  2.5-16c). 
The  basic  nature  of  this  pulse  is  subsonic;  acting  in  accordance  with  the 
local  effective  Mach  Number  M  and  not  the  local  supersonic  Mach  Number  M. 

For  this  case,  the  pulse  passes  over  the  receiving  point  once.  Thus,  the 
subsonic  Doublet  Lattice  Method  is  modified  for  this  case  using  (x  -  l^)  and 
M.  Refering  to  Figure  2.5-16d,  it  is  seen  that  pulses  that  start  out  in  a 
supersonic  zone  have  a  dual  character:  a  supersonic  character,  where  the 
pulse  passes  over  those  receiving  points  twice  that  are  within  the  Mach 
cone  (x^,  z^)  (not  quite  a  cone  for  this  case)  and,  a  subsonic  character,  as 
the  pulse  breaks  out  of  the  Mach  cone  and  moves  upstream  to  pass  over  points 
that  lie  outside  of  the  Mach  cone  (X2,  Zg).  The  receiving  point  (x^,  z^) 
then  has  the  pulse  passing  over  it  three  times.  An  assumption  currently 
adopted,  which  could  be  eliminated,  is  that  the  third  passing  of  the  wave  is 
neglected  for  points  lying  within  the  Mach  cone.  The  single  passing  of  the 
wave  at  (X2,  Z2),  however,  is  retained.  Figure  2.5-17  illustrates  this  case 
by  presenting  the  time-of-arrival  versus  airfoil  location  for  an  acoustic 
source  located  at  x/c  =  0.4  for  the  NLR  7301  airfoil  at  M  =  0.745  and  a  =  -0.20J 
In  the  area  lying  between  0.53  >  |l0.4,  there  are  three  values  for  the  time- 
of-arrival  of  the  acoustic  pulse.  Currently, in  the  present  method  the  highest 
of  the  values  of  (t  -  t)  is  ignored  in  this  area.  There  is  some  justification 
for  this  since  the  pulse  wave  strength  weakens  as  time  increases.  However, 
other  factors,  like  the  speed  with  which  the  wave  passes  over  the  receiving 
point,  have  a  large  effect  on  the  magnitude  of  the  pulse  contribution.  The 
second  highest  value  for  the  time  of  arrival  is  accounted  for  implicity  by 
using  the  first  time  of  arrival  to  obtain  (x  -  l^)  and  M  and  assuming 
supersonic  flow.  This  means  that  the  second  time  of  arrival  is  legislated 
using  (x  -  l^)  and  M  and  is  not  quite  equal  to  the  actual  value.  The  relation¬ 
ship  between°the  two  times-of-arrival  (assuming  uniform  flow  quantities  of 
V,  a,  etc.)  is 
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M  +  1 
M  -  1 


r 


(t  -  1)2  =  (t  -  t)^ 

\  / 

The  result  Is  shown  In  Figure  2.5-17  as  the  dashed  line.  This  assumption  can 
also  be  eliminated  by  breaking  up  the  supersonic  evaluation  into  its  two 
component  parts  and  applying  a  separate  time-of-arrival  and  local  effective 
Mach  Number  for  each  part.  Currently,  this  refinement  does  not  seem  warranted 
but  may  be  in  the  future. 

In  summary,  the  equivalent  (x  -  Tq)  and  M  are  used  in  the  classical 
formulas  for  either  the  supersonic  (M  >  1)  or  subsonic  (M  <  1)  downwash 
influence  coefficients  for  every  sending  receiving  (x)  pair  of  points. 

The  classical  method  used  is  the  Doublet  Lattice  Method.  This  simple  way  of 
using  the  equivalent  data,  i.e.,  (x  -  l^)  and  M  is  approximate  in  some  ways 
but  at  this  stage  of  development  this  seems  appropriate.  There  is  no  use 
refining  a  method  before  it  has  been  tried. 

This  method  was  applied  to  the  NLR  7301  airfoil  at  its  shockfree  design 
point.  Figure  2.5-15(a)  presents  the  values  of  M  used,  while  Figure  2.5-18 
presents  the  time-of-arrival  data  used.  Figure  2.5-19  presents  a  comparison 
for  the  steady  case  of  the  present  method  and  the  Jameson  Method.  The  results 
for  the  Jameson  method  are  obtained  by  taking  a  numerical  derivative  of  the 
results  at  a  as  follows: 

\  "  j  Cp  (a  =  ^  +  )  .  Cp  (a  =  5  -  |2-  )j/  La 

Figure  2.5-19  presents  data  for  the  upper  and  lower  surfaces.  On  the 
upper  surface  the  present  method  is  compared  to  the  numerical  derivative  of 
the  Jameson  method  (ao  =  0.10®)  for  two  average  angles  of  attack,  a,  that  are 
only  slightly  different.  The  difference  in  the  Jameson  results  for  a  of 
-0.2°  and  0.175*15  large.  The  large  response  of  to  small  a  changes  is 
not  currently  understood.  Because  of  this,  the  derivative  itself  may  be  in 
doubt.  The  present  method  is  in  better  agreement  with  the  Jameson  method 
than  is  the  classical  result  (which  is  the  negative  of  the  classic  result 
for  the  lower  surface).  It  is  the  opinion  of  the  author  that  the  large 
response  of  the  Jameson  method  to  small  changes  in  a  can  not  be  duplicated  in 
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Figure  2.5-19  Comparison  of  the  Present  Method  and  the  Classical 
Theory  with  the  Jameson  Transonic  Method  for  Two 
Average  Angles  for  the  Case  of  Shock  Free  Flow 


practice  by  the  present  method  with  Its  current  set  of  assumptions.  The 
present  method  depends  on  the  Jameson  methci  at  a  but  It  does  not  seem  possible 
that  the  flow  field  would  change  enough  (by  changlio  ^  from  -0.2  to  .17)  to 
change  the  entire  character  of  the  acoustic  pulse  solutions.  If  such  changes 
do  occur,  It  seems  that  a  much  more  accurate  representation  of  the  flow  field 
must  be  used.  In  the  acoustic  pulse  generator,  than  the  one  In  use  (double 
linear  Interpolation).  It  seems  that  a  fruitful  line  of  Investigation  would 

be  to  try  and  uncover  the  mechanism  of  this  rapid  change  of  character  with  a 
using  a  finite  difference  method  like  the  Jameson  method  or  better  yet,  the 
Traci -Farr -Albano  method  (Reference  8). 

On  the  lower  surface  of  the  NLR  7301  airfoil,  the  classic  theory  and 
present  method  are  compared  with  results  from  the  Jameson  Method.  As  expected, 
this  surface  exhibits  less  transonic  effects  than  the  upper  surface.  Also, 
changing  a  from  -0.2"  to  0.17®  does  not  effect  the  result  to  any  noticeable 
degree.  The  results  for  the  present  method  agree  fairly  well  with  the  Jameson 
method. 

Application  of  the  present  method  and  Jameson  method  to  the  NLR  7301 
airfoil  for  an  off  design  case  (M^  *  0.7,  a  ‘  1.75®)  Is  presented  In  Figure 
2.5-20.  In  this  case,  a  shock  wave  exists  at  x/c  =  0.4  (see  Figure  2.5-22). 
Also,  shown  In  this  figure.  Is  the  classic  uniform  flow  subsonic  result  which 
bears  no  resemblance  to  the  Jameson  result.  The  results  given  by  the  present 
method  do  not  correlate  with  the  Jameson  results  especially  in  the  region  of 
the  shock  wave.  Notice  especially  that  the  large  peak  In  pressure  at  the 
shock  location  Is  missing.  In  addition,  the  large  negative  peak  in  pressure 
Is  missing.  What  does  exist  is  a  pressure  that  peaks  at  a  value  that  is 
approximately  the  average  of  the  two  peaks  however  which  lies  midway  between 
them.  This  Indicates  that  shock  wave  motion  effects  are  missing  from  the 
present  method.  This  could  easily  be  the  case  since  no  boundary  or  compati¬ 
bility  conditions  were  Introduced  at  the  shock  location.  The  subject  of  the 
next  subsection  Is  the  Inclusion  and  satisfaction  of  shock  wave  motion 
boundary  and  compatibility  condUlons  In  the  present  method. 

Before  proceeding  to  these  subsections.  It  will  be  instructive  to  compare 
the  Jameson  method  with  experimental  data  gathered  for  the  NLR  7301  airfoil 
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Figure  2.5-20  Comparison  of  the  Present  Method  (Without  Shock 
Wave  Effects)  and  the  Classic  Theory  With  the 
Jameson  Transonic  Method 
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Comparison  of  The  Jameson  Transonic  Flow  Method  With  Data  for  The  NLR  7401  Airfoil 
(M  ■  0.745,  a»-0.2O).  Angle  of  Attack  for  Theory  Adjusted  for  Best  Correlation. 


—  DC^  Wind  Tunnel  Data 
(Ref.  4  )  a  =  0.85° 


o 
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o 

Figure  2.5-23  Comparison  of  the  Jameson  Transonic  Method  With  Wind 
Tunnel  Data  for  The  NLR  7301  Airfoil.  Jameson  Method 
Used  at  Two  Angles  of  Attack 
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and  presented  in  Reference  4.  Figures  2.5-21  and  -22  present  comparisons  of 
wind  tunnel  data  with  the  Jameson  Method  for  the  two  conditions  previously 
studied.  The  angles  of  attack  have  been  adjusted  to  reproduce  the  best 
comparison.  A  better  method  of  matching  however,  is  to  add  the  boundary 
layer  thickness  to  the  airfoil  and  rerun  the  transonic  calculations.  A 
second  method  is  to  perform  an  inverse  calculation  whereby  the  required 
pressure  distribution  is  input  and  the  modified  airfoil  shape  it  output.  This 
latter  technique  was  applied  in  the  cases  of  the  TF-8A  and  YC-15  wing  sections 
(to  be  discussed  later)  but  was  not  applied  in  this  case.  Figures  2.5-23 
and  -24  compare  data  with  the  Jameson  method  for  LC^/^a.  These  figures  show 
that  the  pressure  levels  output  by  the  Jameson  method  are  approximately 
correct  and  can  be  used  as  a  basis  for  comparison  with  the  present  method. 
These  figures,  especially  2.5-23,  show  that  the  transonic  perturbation  press¬ 
ures  are  very  volatile,  varying  rapidly  with  the  flow  parameters.  In  Figure 
2.5-24  note  especially  the  disagreement  between  the  Jameson  method  and  data 
aft  of  the  shock. 

The  question  may  arise:  if  the  Jameson  program  predicts  the  perturbation 
loads  adequately,  why  not  use  the  Jameson  program  alone  and  not  worry  about 
the  lifting  surface  method.  First,  the  Jameson  program  is  applicable  only  to 
steady  flow.  Second,  the  extension  of  the  Jameson  method  to  three-dimensional 
flows  will  increase  the  cost  significantly, whereas  the  lifting  surface  theory 
costs  are  very  reasonable. 

2.5.2  Effects  of  Shock  Wave  Motions 

Figure  2.5-20  illustrates  that  the  present  method,  as  developed  up  to 
this  point,  does  not  account  for  shock  wave  motion.  Inspection  of  this  figure 
shows  that  the  shock  motion  effect  plays  a  dominate  role  in  the  prediction  of 
transonic  loads.  Thus  in  the  next  few  subsections  the  theory  for  the  effects 
of  shockwave  oscillation  will  be  developed.  Along  these  lines  References  10, 
12  and  13  are  helpful  in  the  following  developments. 

2. 5. 2.1  Interference  Flow  Due  to  Shock  Motion 

The  following  jump  conditions  across  a  shock  wave  can  be  derived  from 
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Reference  14  (page  59  Eq.  2.47). 


where  u  and  w  are  the  x  and  z  components  of  the  velocity,  the  tilda  stands  for 
instantaneous  position  and  the  subscripts  stand  for  locations  just  up  and 
downstream  of  the  shock.  Therefore  these  shock  jump  conditions  hold  just  up 
and  downstream  of  the  instantaneous  shock  position.  In  order  to  translate 
these  conditions  to  the  steady  shock  surface  an  analytic  continuation  of  the 
steady  conditions  is  necessary.  Consider  the  velocity  potential  for  instance: 


A  +  e^“^  X  H  I  3x  +  2 

^1,2  ^1.2' 


where  <|i  is  the  steady  potential  just  up  {i>  )  or  downstream  (4>g  )  of  the 

^12  ^  ^  I, 

steady  shock  position,  x  is  the  amplitude  of  oscillatory  motion  of  the  shock, 
and  2  oscillatory  perturbation  potential  just  up  or  downstream  of 

the  steady  shock  location  (see  Figure  2.5-25).  Here  we  have  made  the  assump¬ 
tion  that  the  shock  moves  in  a  simple  harmonic  manner.  Since  w  =  94>/9z  the 
tangential  flow  boundary  condition  ^2  -  =0  can  be  written  as  follows 

using  the  above  expression  for  ?i^2' 


But  since  w  =  9^/9z  then 


steady 

Oscillatory 


9W  /9X  =  924,  /9Z9X  =  9U  /9Z 

S  5  ^ 

Thus  the  oscillatory  shock  jump  condition  can  be  rewritten  as: 

«2  -  “1  '  -’'h 
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Figure  2.5-25  Analytic  Continuation  of  Flow  Variables  Due  to 
Shock  Wave  Motion 


This  equation  shows  that  vorticity  exists  on  the  shock  surface  because  there 
Is  a  jump  In  ti^e  tangential  velocity,  w,  across  the  shock.  Let  y  be  this 
vorticity,  then 


=  -(Wg  -  w^)  =  A 


G  -  D  ) 

az  \  S2  / 


If  we  use  the  linear  expression  for  the  steady  pressure,  i.e.,  Cp  =  -2  u^/U^ 


then  the  vorticity  y  can  be  expressed  in  terms  of  pressure  jump  across  the 
steady  shock  wave. 


u  /  2 


(2. 5. 2-1) 


The  quantities  C  (z)  and  C  (z)  are  known  from  the  steady  solution,  there- 
fore  y/a  is  known. 


This  vorticity  distribution  Induces  a  flow  field  (see  Figure  2.5-26)  in 
the  vicinity  of  the  airfoil  which  must  be  negated  by  the  doublet  lattice 
elements  that  are  on  the  lifting  surface.  Cunningham  (Reference  12)  came  to 
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Figure  2.5-26  Shock  Vorticity  Interference  Flow 


the  conclusion  that  a  concentrated  vortex  was  to  be  placed  at  the  foot  of  the 
shock  wave  to  simulate  the  shock  motion  interference.  However  from  the  above 
derivation  it  is  seen  that  such  a  concentrated  vortex  is  an  approximation  of 
the  distributed  vorticity  that  lies  on  the  shock  surface.  The  strength,  r, 
of  this  equivalent  concentrated  vortex  is: 

^sonic  /  \ 

f  =  ^  9dz  =  X  (  Cp  (0)  -  Cp  (0)|u^  /  2  (2. 5. 2-2) 

6  ^  ^2  ^1 


since  C„  (Zcnnir^ 


some' 


across  the  shock  at  its  upper  end, 


This  is  because  the  pressure  jump 

is  zero.  The  quantities  Cp  (0), 

^2 


C  (0)  are  the  steady  pressure  coefficients  on  the  surface  (z  =  0)  in  front 
^^1 

C  ,  and  in  back  of,  ,  of  the  steady  shock. 

Pc.  Ps. 


For  simplicity  in  the  present  analysis  the  point  vortex  model  will  be 
adopted.  The  onset  downwash  flow  field  due  to  this  point  vortex  can  be 
computed  with  the  aid  of  the  equivalent  geometry,  (x  -  ).  and  Mach  Number 

R  where  the  subscript  "sv"  stands  for  "shock  vortex''.  SoS^  question  arises, 
however,  as  to  what  values  of  time  delay  (t  -  t)  should  be  used  in  this  down- 
wash  calculation  and  if  the  point  vortex  is  to  be  placed  on  or  off  of  the  air¬ 
foil  surface.  At  present  it  will  be  assumed  that  the  point  vortex  lies  on  the 
airfoil  surface  and  that  the  time  delay  to  be  used  is  that  which  is  associated 
with  this  surface  point.  Figure  2.5.27  presents  the  downwash  due  to  a  unit 


f 

/ 


Figure  2.5-27  Shock  Vortex  Interference  Flow  Downwash  For 
The  NLR  7301  Airfoil  Operating  at  Mq,  ■  0.7 
and  a  ■  1.75° 
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point  vortex  on  the  surface  of  the  NLR  7301  airfoil  at  the  shock  location 
x/c  =  0.4  for  the  off-design  case  of  =  0.7.  Shown  in  this  figure  are  two 
curves;  one  representing  the  downwash  due  to  a  transonic  vortex  and  one 
representing  the  downwash  due  to  a  classic  subsonic  vortex  in  a  uniform 
stream  (of  M^=  0.7)  for  reference.  As  a  check  on  t  -  t  for  the  shock  wave 
vortex  the  acoustic  pulse  generator  was  rerun  for  a  series  of  points  above 
the  airfoil  along  the  shock  surface.  The  results  are  that  the  time  delay, 
t  -  T,  did  not  change  significantly  for  points  along  the  length  of  the  shock 

wave. 


For  the  NLR  7301  airfoil  the  shock  motion  x  is  known  because  of  experi¬ 
mental  data  and  data  from  the  Jameson  method.  Thus  the  vortex  strength,  r, 

is: 

where  3X/3o  is  the  rate  of  change  of  shock  motion  from  its  steady  position 
with  respect  to  a.  In  general  this  term  is  time  complex  but  in  the  steady 
case  to  follow  it  is  simply  a  real  number.  For  the  NLR  airfoil  at  0.7, 
a  =  1.75°  the  Jameson  method  gives  3X/3a|Cp  (0)  -  Cp  (0))=  6.02.  Thus 

^2  ^1 

f  =  3.01  for  this  airfoil  and  condition.  The  resulting  interference  pressure 
is  presented  in  Figure  2.5-28.  Large  variations  are  observed  between  the 
classic  subsonic  vortex  and  the  transonic  one.  This  volatibility  indicates 
that  the  details  of  the  vorticity  on  the  shock  may  be  important.  Future 
refinements  should  then  include  refinements  in  the  modeling  of  the  shock 

vorticity. 

2. 5. 2. 2  Pressures  Due  to  Shock  Wave  Motion  Interaction  with  Steady  Flow 

The  pressures  C  caused  by  the  shock  vortex  interference  flow  (shown  in 

P  Qt 

Figure  2.5.28)  is  to  be  added  to  the  pressure  calculated  in  the  absence  of 
the  shock  (shown  in  Figure  2.5-20).  Before  this  is  done  however  there  is  one 
more  term  that  must  be  considered.  This  last  term  in  the  pressure  is  caused 
by  the  direct  action  of  the  shock  wave  on  the  surface.  As  the  shock  wave 


77 


IHiS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELIABE 

UNDER  DOD  DIRECTIVE  S200.20  AND 

Nu  RESTRICTIONS  ARE  IMPOSED  UPON 
!TS  USE  AND  DISCLOSURE. 

DISTRIBUTION  STATEMENT  A 

# 

APPROVED  FOR  PUBLIC  RELBASE; 
DISTRIBUTION  UNLIMITED. 


passes  over  a  pafticular  point  the  perturbation  pressure  rises  sharply.  The 
region  on  the  airfoil  where  this  occurs  is  confined  to  that  area  over  which 
the  shock  wave  moves. 


The  pressure  coefficient  Cp  on  the  airfoil  due  to  the  motion  of  the  shock 


wave  IS 


Cp(a.  X)  •  Cp  -  Cp 

\  $2  s 


where  H  is  the  step  function  at  the  instantaneous  position  of  the  shock  x^{a) 

and  the  zeros  have  been  dropped  from  C  (0)  and  C  (0)  for  convenience. 

Pc  Pc. 


The  derivative  is: 


(2. 5. 2-3) 


where  6  is  a  dirac  delta  function.  In  the  limit  as  the  amplitude  of  shock 
wave  travel  is  reduced  to  zero  then 

X  (a)  -►  X  =  steady  position  of  the  shock 

9  9 


The  force,  generated  by  this  pressure  is  then 


For  the  case  of  oscillatory  motion  ^  (-iwT^)  and  a  =  aexp(i(i)t). 

If  F/q  =  F/q  exp(iiot)  then: 

/  V 

F/q  =  F^/qS  =  SIC  -  Cp  expl-itp)  (2. 5.2-5) 

\  S2  / 
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where  4>g  and  are  functions  of  the  reduced  frequency  k^.  Notice  that 
equation  2. 5. 2-4  Is  exactly  equal  to  the  force  generated  by  the  shock  vortex 
of  strength  r^.  That  Is: 


3: 


r 


The  expression  for  the  pressure,  Cp  ,  given  in  equation  (2. 5. 2-3),  is 

not  suitable  for  plotting  and  applies  fo?  infinitesimley  small  motions.  For 

finite  motions,  consider  the  finite  form  of  C  : 

*^0 


(2. 5. 2-6) 


This  is  a  rectangular  shape  of  height  (C  -  C  )/La  and  width 

^2  h 

An  =  3X/3a  Aa.  This  equation  applies  for  shock  waves  that  are  not  defused. 

However  Figure  2.5-22  shows  that  the  shock  is  not  sharp  for  the  NLR  7301 

airfoil  at  =  0.7  and  a  =  1.75®,  but  is  spread  out  from  x/c  =  0.39  to 

x/c  =  0.45;  a  total  length  of  ax^  =  0.06.  An  approximate  way  of  accounting 

for  this  effect  would  be  to  plot  ACp/Aa  as  a  trapezoid  whose  base  is 

AX  +  AX  ,  whose  top  width  is  ax  and  whose  area  is  F  /q.  For  the  NLR  7301 
5  o  ® 

airfoil  at  =  0.7,  a  =  1.75  the  value  of  ax  is  0.065  for  a  corresponding 

Aa  of  0.0087  rad.  The  trapezoid  base  is  then  0.125,  the  top  width  is  0.065. 

The  total  force,  9X/3a  (C  '  height  of  the  trapezoid 

Ps2 

is  6.02(.065  +  .175)/2  =  63.4. 


All  of  the  contributing  terms  for  the  final  pressure  distribution  are 
now  known,  that  is,  1)  the  shock  free  pressure  distribution  derived  in  sub¬ 
section  2.5.1,  2)  the  pressure  distribution  caused  by  the  interference  flow 
of  the  shock  vortex,  and,  3)  the  pressure  distribution  just  discussed  which 
is  a  direct  result  of  the  moving  of  the  shock  wave.  Specific  application  to 
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the  NLR  7301  airfoil  at  =  0.7  and  o  =  1.75®  is  presented  in  Figure  2.5-29, 
where  the  results  of  Figure  2.5-20  (shock  free  flow)  are  combined  with  the 
results  of  Figure  2.5-28  (transonic  shock  vortex)  and  with  the  trapezoidal 
distribution  discussed  above.  The  agreement  is  good  except  in  the  region 
forward  of  the  shock  wave  in  the  supersonic  zone.  The  good  agreement  aft 
of  the  shock  may  not  be  entirely  trustworthy  since  it  appears  that  the  shock 
vortex  interference  pressures  are  presently  somewhat  in  doubt  and  this  is  the 
flow  that  produces  this  correlation.  The  correlation  with  the  shock  peak  ,is 
good  because  the  present  method  has  been  constrained  to  reproduce  the  proper 
total  force  F  /q  due  to  steady  shock  motion.  This  total  force  is  proporational 
to  shockwave  motion,  3A/8a,  and  steady  pressure  jump  (C  -  C  ).  Currently 

Ps2  Ps^ 

the  present  method  uses  either  experimental  data  or  results  from  the  Jameson 
steady  transonic  flow  method  to  obtain  the  quantities  needed,  i.e.,  ax/8a 
and  C  ■  subsections  describe  a  possible  method  for  deter- 

Ps2  Pst 


mining  the  shock  wave  amplitude  3X/ao  for  the  general  oscillatory  case  (which 
includes  the  steady  case).  The  Jameson  method  can  not  be  used  to  obtain  shock 
wave  motion  data  for  the  oscillatory  case  although  it  can  be  used  to  aid  in 
its  determination. 


The  object  of  the  present  method  is  to  furnish  an  inexpensive  lifting 
surface  theory  for  transonic  flow  that  can  be  run  for  various  modes  of  motion 
at  various  frequencies  and  which  is  based  on  steady  flow  data  which  is  used 
once  and  for  all.  This  steady  flow  data  can  be  from  experimental  sources 
or  from  transonic  steady  flow  methods. 


So  far  in  the  correlation  of  the  present  method  only  steady  flow  conditions 
have  been  considered.  However  at  this  stage  of  development  this  seems 
appropriate  since  the  steady  perturbation  results  must  agree  before  the  un¬ 
steady  results  are  attempted.  In  addition  the  unsteady  shock  wave  amplitude 
calculation  has  not  yet  been  tried  in  the  present  method  and  therefore  a 
complete  correlation  awaits  this  development. 

The  correlation  illustrated  in  Figure  2.5-29  shows  that  the  shock  wave 
motion  contribution  to  the  pressure  is  a  dominant  transonic  effect.  A 
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possible  method  of  including  transonic  effects  in  a  simple  way  would  be  to 
add  the  shock  wave  motion  effect  to  the  classic  theory.  This  idea  is  pursued 
in  the  next  subsection. 

2. 5. 2. 3  Investigation  of  a  Simple  Method  for  Accounting  fcr  Transonic  Effects 

In  this  subsection  the  results  of  adding  the  shock  motion  pressures  to 
those  determined  by  the  classical  subsonic  theory  is  investigated  for  the 
NLR  7301  airfoil.  The  pressure  distribution  is  presented  in  Figure  2.5-30 
for  reference. 

In  the  analysis  of  flutter  and  other  dynamic  aeroelastic  phenomena 
it  is  the  section  lift  and  moment  which  are  of  paramount  importance.  In 
Reference  4  the  experimental  lift  and  moment  are  plotted  versus  the  reduced 
freguency  for  various  Mach  Numbers.  The  results  of  the  classic  theory  are 
compared  to  these  experimental  values.  At  low  Mach  Numbers  the  data  and 
theory  should  agree  since  transonic  effects  are  small.  However  there  is  still 
a  difference  between  theory  and  data  and  this  difference  is  attributable  to 
viscous  boundary  layer  effects.  This  boundary  layer  effect  can  be  isolated 
by  subtracting  the  results  for  the  experiment  from  those  of  the  theory  for 
the  low  Mach  Number  case.  This  effect  can  then  be  taken  out  of  the  data  for 
transonic  Mach  Numbers  so  that  viscous  effects  do  not  confuse  the  issue. 


Transonic  effects  can  be  easily  isolated  now  that  the  viscous  effects  are 
eliminated.  This  is  accomplished  by  taking  the  difference  between  the  classi¬ 
cal  theory  and  the  corrected  experimental  data  in  the  transonic  speed  range. 
This  has  been  done  for  the  NLR  7301  airfoil  at  M^  =  0.7,  a  =  1.75  and  the 
results  plotted  in  Figure  2.5-31.  The  specific  formula  for  the  corrected 

transonic  increment  in  lift  ac^^^  is: 

a 

=  c.  (data)  -  c.  (classic  theory)  at  =  0.7 

where  the  corrected  data,  .  is  the  result  of  eliminating  the  viscous  loss 

of  lift,  ACj^  ,  from  the  data? 
a 

c.  (data)  =  (data)  +  Ac^^  at  M^  =  0.7 

*■0  a  a 


83 


•  Jameson 


Simplified  Theory 


Comparison  of  The  Jameson  Method  With  The 
Simplified  Transonic  Method 


•  BOO  Experimental  Transonic  Effects 
•  ■  Real  on  imaginary 


Effects  Due  to  Shock  Wave  Motion 

Steady  Shock  Strength  &  Location 
Based  on  Jameson  Method 
Steady  Shock 

Strength,  etc.,  Using  Data 


Tnifintry 


Comparison  of  The  Simple  Transonic  Method  With 
Experimental  Data.  Simple  Method  Uses  Steady 
Data  From  Two  Sources. 


Figure  2.5-31 


where  the  viscous  loss  of  lift  Ac  Is  determined  from  theory  and  data  at  M  =  0.5. 

X  00 

a 

AC  =  c.  (theory)  -  c,  (data)  at  M  =  0.5. 

XX  X  00 

a  a  01 

A  similar  procedure  applies  to  the  moment. 


This  data  is  compared  with  the  simplified  transonic  effect  consisting  of 
the  shock  motion  pressure  alone.  Thus 


a 


_  3^  -  C  ) 

^2  ^1 


/V/ 

AC 


m 


Xj/C 

a 


where  F  /qc  is  the  normalized  force  due  to  shock  wave  motion  and  x  is  the 
a  _  s 

shock  location  point  at  which  the  force  acts. 


The  term  X  is,  in  general,  a  complex  number  which  can  be  represented  as 
follows  (see  Equation  2. 5. 2-5): 


A  =  XQ(k^)exp(-i(<)j) 

where  is  the  phase  lag  between  the  shock  and  airfoil  motions.  Reference  4 
shows  that  this  phase  lag  is  proportional  to  the  frequency.  Based  on  an 
approximate  calculation  it  v/as  concluded  that  the  constant  of  proportionality 
was  the  time,  t,  it  takes  a  "Kutta  wave"  to  reach  the  shock  wave.  Reference 
4  defines  the  Kutta  wave  as  an  acoustic  wave  eminating  from  the  trailing 
edge.  Using  the  acoustic  wave  generator  described  in  subsection  2. 5. 1.4  the 
time  history  of  an  acoustic  wave  eminating  from  the  trailing  edge  was  gener¬ 
ated.  The  results  are  presented  in  Figure  2.5-32.  The  time  that  the  shock 
first  fee'is  this  acoustic  wave  is  approximately  t  =  tU^c  =  3.0.  The 
experimental  data  of  Reference  4  shows  that 


<1.5  =  4.32  k^ 
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If  it  is  assumed  that  (j)^  is  also  proportional  to  the  time  t  then  the  prop¬ 
ortionality  constant  is  1.44  or 

i|)g  =  1.44  T  (concluded  by  present  authors) 

The  authors  of  Reference  4  concluded  that 


because  they  had  only  an  approximate  value  for  x.  Irrespective  of  these  two 
formulas  the  phase  lag,  4)^,  is  taken  as  the  experimental  value  of  4.32. 
(Actually,  for  Figure  2.5-31,  the  value  4)^  =  4.56  was  used.  The  difference, 

5%,  should  not  cause  any  significant  difference). 

Reference  4  also  4)resents  experimental  data  (for  the  NLR  7301  airfoil) 

for  the  amplitude  of  the  shock  wave  motion,  x^,  as  a  function  of  frequency. 

This  data  can  be  generalized  by  normalizing  it  with  respect  to  the  steady 

amplitude  X^  (k^,  =  0).  A  plot  of  ^Q(l<j,)/^Q(k|.  =  0)  =  G  shown  in  Figure 

2.5-33.  This  curve  will  be  used  later  for  the  TF-8A  and  YC-15  wings.  The 

value  of  X  (k^)  is  obtained  from  this  curve  as  follows: 

0  r 

^^(kr)  =  x^(k^  =  0)  G  (k^) 

For  Figure  2.5-31  two  values  of  ^^(k^  =  0)  were  used:  the  experimental  value 
and  the  value  obtained  from  the  Jameson  method. 

The  agreement  for  the  incremental  transonic  effects  is  fairly  good 
considering  the  simplicity  of  the  theory.  The  percent  error  due  to  the 
differences  in  the  lift  increments  will  be  reduced  when  added  to  the  basic 
lift.  For  the  moment  however  the  shock  movements  represent  the  largest  portion 
of  the  final  result  and  thus  the  percent  error  will  not  change  much.  In 
general  it  is  observed  that  the  theory  over  predicts  the  lift  and  under  pre¬ 
dicts  the  moment. 

2. 5. 2. 4  Shock  Wave  Compatibility  Relations 
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In  subsection  2.5.2. 1  the  tangential  velocity  shock  jump-condition  was 
used  to  obtain  the  interference  flow  due  to  the  shock  motion.  It  was  shown 
that  this  interference  flow  is  caused  by  a  vorticity  distribution  of  known 
shape  which  is  proportional  to  the  shock  wave  amplitude.  The  shape  of  the 
vorticity  was  simplified  to  a  point  vortex  on  the  surface.  The  total  vorticity 
is  proportional  to  the  steady  component  of  pressure  rise  across  the  shock 
wave,  (C  -  C  ),  and  the  amplitude  of  shock  wave  motion,  3X/8a.  The  steady 
^2  '^^1 

pressures  are  known  but  the  amplitude  of  unsteady  shock  motion  is  not  Known 
and  it  is  the  object  of  this  section  to  suggest  a  procedure  for  obtaining  it. 


The  first  equation  of  subsection  2. 5. 2.1  presents  the  jump  conditions  for 
the  x-component  of  velocity  at  the  instantaneous  shock  position: 


-2  'v- 


1 


where,  as  before,  the  tilda  indicates  instantaneous  shock  position  and  the 
subscripts  indicate  positions  just  upstream,  1,  and  just  downstream,  2,  of 
the  shock.  Landahl  (Reference  15,  Equation  1.17,  page  4)  gives  an  expression 
for 

— - —  =  — ^  +  (1-y)[^ 

-V/  2  y,  2  I  ^  '^X,z'-'  '  <» 

M,  M 

1  “ 


where  0  (  )  indicates  order  of  magnitude.  The  following  analytical  continu 

ation  of  ^  is  used: 


^  =  <{)^  +  (A3(^ax  +  1)  exp(ia)t) 

'v  2 

If  the  proper  derivatives  are  taken  of  ^  and  the  nonlinear  terms,  0((^  ) 

X  >z 

are  assumed  small  enough  to  ignore  then 


+  (1  - 

^1 
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where  (|>j  and  ^  are  the  steady  and  oscillatory  components  of  the  potential. 

Applying  the  analytic  continuation  procedure  to  u^  gives 
u,  =  u  +  (x  3u^  /3x  +  uJ  exp(lut) 

1  Si  ' 

Introducing  the  expressions  for  u^,  and  1/M^^  "into  the  expression  for 
U2  -  u^  and  remembering  that  ^/ix  ■  u  gives: 

^2  ■  “1  ^  j A  +  [  XB  +  Cu^  +  D(|.^  ]  exp{i(i)t)  +  0>xp(i2a)t))  j 


where 

A  =  -u*  +  (1  -  y)  u.  ) 

s^  »  « 

B  •  -u,  (1  -  y)  (u,  +  1«u  )  -  u  (B^  /M^  +  (1  -  y)  u  ) 

'‘l  ®lx  ^1x  1 

C  •  -u  (1  -  y)  -  0  -  y) 

D  *  -u.  (1  -  y)  I0) 

^1 


Expanding  the  left  hand  side  of  the  equation  U2  -  u^  and  using  the  expression 
for  analytic  continuation  gives: 


U2  -  u^  =  E  +  (XF  +  U2  -  u^)  exp(ia)t) 


where 


E  =  u  -  u 
S2 


F  =  u  -  u 

^2x  ^Ix 


Equating  the  two  expressions  for  U2  -  u^  and  separating  steady  and 
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oscillatory  parts  gives: 


u 


^2 


(  i  /m!  +  (1  - 


steady  flow 


■  ^2  '  ^  ^  (2. 5. 2. 4-1) 

oscillatory  flow 

These  are  the  final  expressions  for  the  u-component  jump  conditions  across 
the  shock  referred  to  the  steady  position  of  the  shock. 


The  jump  conditions  for  the  w-component  were  given  in  subsection  2. 5. 2.1  as 


w.  -  w,  =  -  X3(u^  -  u^  )/dz 


'2  "1 


S2 


If  this  expression  is  integrated  clockwise  around  the  shockwave  then 

isonic 

^2  (z)  -  ^i(z)  =  y  ”i  ^  s  +  J  W2  d  s 


'sonic 


:somc 


'some 


=  -  y^  (W2  -  Wi)ds  =  X  J  bCUj^  -  Us^)/9s  ds 


>2 


4*2(2)  -  =  -  ^(Us^(z)  -  ^s^^z)) 


(2. 5. 2. 4-2) 


since  u  -  u  =  0  at  the  sonic  point  z  =  z  .j..  Here  the  assumption  has 
^2  ^1 

been  made  that  the  amplitude  does  not  vary  with  height  above  the  airfoil. 

If  this  is  not  the  case  then  this  assumption  can  be  eliminated  by  considering 
X  to  be  an  assumed  function  of  z  of  unknown  amplitude.  Equations  (2. 5. 2. 4-1) 
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and  (2. 5. 2. 4-2)  can  be  combined  to  eliminate  the  unknown  shock  wave  motion, 
X  as  follows: 

G  -  Ug  H  =  J  -  ^2  (2. 5. 2. 4-3) 


where 

/ 

G  =  [1-2  c/(y+1)]  (u 

S,  -  \) 

J  =  [F  +  2B/(y+1)]  +  ( ‘J's2  ‘ 

K  =  F+2B/(y+1) 

In  terms  of  pressures  C-  =  )•  Equation  (2. 5. 2. 4-3)  can  be 

P  \u»  '  written  as: 

■^GC  +^HC  •  (0  +  fg-S)  -  *2  (K  +  tr”'  (2-5.2-«-'t) 

2  Pi  c  ^2  '  "  " 

Equations  (2. 5. 2. 4-3)  or  (2. 5, 2. 4-4)  furnishes  a  compatibility  condition 
that  is  to  be  satisfied  by  the  oscillatory  flow  across  the  steady  shock 
position.  These  equations  hold  at  all  positions  along  the  shock  from  the 
surface  (z  =  0)  to  the  sonic  point  (z  =  Consistent  with  the 

assumption  that  the  shock  amplitude  of  motion  does  not  change  along  the 
shock  is  the  assumption  that  the  compatibility  equation  is  not  a  function 
of  the  vertical  direction  z.  These  equations  can  thus  be  satisfied  on  the 
airfoil  surface  only.  In  this  case  the  second  version  of  the  compatibility 
equation,  i.e.  (2. 5. 2. 4-4)  is  appropriate.  This  expression  furnishes  one 
additional  equation  in  the  set  of  equations  for  the  lifting  surface  pressures. 
The  wave  amplitude,  x,  is  the  one  additional  unknown  in  the  problem.  Thus 


the  problem  is  determlnent. 


The  lifting  surface  theory  without  a  shock  wave  (but  with  transonic 
effects)  can  be  written  in  matrix  form  as: 


[D]  {  I 


where  the  subscript  "o"  indicates  shock  free  transonic  flow.  The  shock 
vortex  onset  flow  is  given  by 


i  P 


I  = 


r  I 


[D]|w^ 


where  {Wj,}  is  the  downwash  due  to  the  shock  vortex  or  vorticity.  The  strength 
of  the  vortex  is  that  associated  with  shock  motion  of  unit  amplitude.  The 
compatibility  equation  can  be  written  as: 

where  the  term  {L}^  is  the  matrix  form  of  Equation  (2. 5. 2. 4-4)  and  where: 
j  I  =  (  ^  !  Cn  ! 

I  Pt  )  I  Po  )  I  Pr  i 


Combining  the  first  two  equations  with  the  third  produces  an  equation  that 
can  be  solved  for  the  shock  motion  amplitude  x.  Placing  x  into  the  expression 

produces  the  final  total.  This  total  represents 


i  r.  I 


for  Und  adding  '  C„ 

I  Pri  I  Po  ) 

the  combination  of  the  shock  free  transonic  and  shock  vortex  flows.  To  this 
is  added  the  shock  motion  pressure.  It  is  anticipated  that  the  steady  value 
of  X  can  be  used  as  an  aid  in  the  evaluation  of  the  quantities  required  in 
the  compatibility  equation. 


In  sumnary,  a  procedure  has  been  outlined  for  determining  the  shock  wave 
motion  amplitude  for  oscillatory  flow.  This  procedure  has  not  been  implemented 
under  the  current  contract  or  McDonnell  Douglas  Company  IRAD  funds.  Future 
investigations  of  this  procedure  are  thus  recommended. 
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2.6  FLUTTER  ANALYSIS  PROCEDURES 


All  of  the  flutter  analyses  performed  in  this  study  were  carried  out 
using  the  standard  American,  or  V-g  method.  The  analyses  used  the  measured 
aerodynamic  mode  shapes,  generalized  mass,  and  frequencies  of  six  elastic 
modes  of  the  TF-8A  flutter  model.  A  few  cases  were  also  analyzed  using  modal 
deflection  data  as  calculated  from  basic  distributions  of  inertial  and  rigidity 
data.  Section  3  below  discusses  the  basic  data  used  in  the  flutter  analysis 
in  detail . 

Analyses  were  run  at  Mach  numbers  of  .6,  .8,  .9,  .95  and  .99  for  the 
TF-8A  flutter  model  and  .76  for  the  YC-15II.  In  each  case  the  unsteady 
aerodynamic  influence  coefficients  (AIC's)  were  generated  theoretically  by 
the  subsonic  Doublet  Lattice  Method  and  weighted  appropriately. 

Two  sorts  of  weighting  were  applied  to  each  case;  they  are  described 
in  detail  in  Sections  2.4  and  2.5.  Briefly,  the  first  method,  loosely  des¬ 
cribed  as  "steady  corrections"  assumes  that  the  ratio  of  the  correct  data  to 
the  theoretical  data  at  all  values  of  reduced  frequency,  k,  is  the  same  as  at 
k  =  0.  The  correct  data  at  k^  =  0  is  of  course  the  available  static  pressure 
model  data,  converted  to  spanwise  distribution  of  sectional  lift  and  moment 
coefficient.  This  method  does  account  for  the  static  aeroelastic  twist  and 
the  non-linear  character  of  the  aerodynamic  derivatives.  The  second  method, 
described  as  "unsteady  corrections"  theoretically  accounts  for  unsteady 
transonic  effects  and  results  in  a  muliplicative  correction  factor  which  is 
applied  to  the  "steady  corrected"  AIC's.  This  multiplicative  correction 
factor  is  a  complex  valued  function  of  k^  and  reduces  to  ''.0  at  k^  -  0. 

The  V-g  solution  was  obtained  for  several  values  of  density  for  each 
Mach  number  so  that  the  calculated  flutter  speed  at  each  density  could  be 
cross-plotted  to  obtain  the  flutter  speed  which  matched  the  Mach  number 
under  consideration.  Structural  damping  was  assumed  to  be  g  =  0.02  in  all 
elastic  modes. 

2.7  FLUTTER  MODEL  PROCEDURES 
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The  analytical  work  documented  in  this  report  shows  the  need  for  further 
wind  tunnel  flutter  model  test  confirmation  of  the  analysis  methods.  Airfoil 
shape,  static  twist,  reduced  frequency  at  flutter,  and  Reynolds  numbers  are 
some  of  the  parameters  that  the  flutter  model  should  evaluate. 

A  cantilevered  wing  (with  adjustable  root  pitch)  would  be  an  ideal 
candidate  for  the  study  since  it  would  avoid  problems  of  scaling  lift  to 
weight  ratio  inherent  in  a  free-flying  model.  The  basic  planform  area  should  be 
one  for  which  high  speed  rigid  model  steady  aerodynamic  data  already  exists, 
such  as  the  proposed  C-15A  wing. 

Instrumentation  should  be  included  to  measure  the  wing  spanwise  static 
bending  and  torsion  deformation.  Optical  methods  can  confirm  the  tip  twist. 

The  Langley  Research  Center  Transonic  Dynamics  Tunnel  provides  the  current 
optimum  facility  for  this  study.  Reynolds  number  effects  can  be  evaluated  by 
building  models  of  difference  size  and/or  air/freon  mixes.  All  of  the  other 
parameter  of  interest  can  be  incorporated  into  a  single  spar  design. 


Section  3.5  below  discusses  some  details  of  a  proposed  flutter  model. 


SECTION  III  results  OF  ANALYSIS 


3.1  PHASE  I  RESULTS 

3.1.1  TF-8A  Vibration  Modal  Data 

The  flutter  analysis  of  the  TF-8A  flutter  model  was  based  upon  measured 
vibration  deflection  data  supplied  by  NASA.  Figure  3.1-1  shews  the  modal 
deflection  measurement  grid.  Table  3.1-1  lists  the  coordinates  of  this  grid. 
Tables  3.1-2  to  3.1-7  show  the  modal  deflections  at  the  deflection  control 
points  for  the  six  measured  TF-8A  inodes.  The  frequency  and  generalized  mass 
of  each  mode  as  supplied  by  NASA  are  also  shown  in  Tables  3.1-2  to  3.1-7. 

In  order  to  more  easily  weight  the  AIC's  it  was  decided  to  transform  the 
idealization  into  plunge  and  pitch  deflections  alono  a  set  of  streamwise 
strips.  This  proved  to  be  justifiable  since  an  inspection  of  the  original 
modal  data  showed  the  wing  to  be  essentially  rigid  in  a  structural  chordwise 
direction.  The  method  used  to  make  this  transformation  is  discussed  above 
in  Section  2.2.  Figure  3.1-2  shows  the  transformed  idealization.  Tables 
3.1-8  to  3.1-13  show  the  transformed  modal  deflections. 

The  vibration  modes  of  the  TF-8A  flutter  model  were  also  calculated 
from  NASA  supplied  rigidity  and  mass  data.  Figures  3.1-3  to  3.1-S  show 
these  data.  The  resulting  mode  shapes  agreed  well  with  the  transformed  test 
values.  For  the  bending  modes,  the  calculated  values  of  frequency  were 
within  five  percent  and  the  generalized  masses  within  four  percent  rf  the 
test  value.  For  the  torsion  modes  a  discrepancy  between  measured  and 
calculated  frequency  and  aeneralized  mass  existed.  Factoring  the  torsional 
stiffness  eliminated  this  discrepancy,  and  when  the  calculated  aerodynamic 
mode  shapes  were  used  in  the  flutter  analysis  along  with  the  measured 
frequencies  and  generalized  m.asses,  the  results  were  within  two  percent  of  those 
obtained  using  the  transformed  measured  deflections. 


Since, at  the  time  the  calculated  modes  were  obtained,  the  study  was 
already  underway,  it  was  decided  to  continue  using  the  transformed  measured 
modal  data  as  the  basis  of  the  flutter  analysis  instead  of  switching  in  mid¬ 
stream  to  the  calculated  shapes.  This  was  justified  in  light  of  the  good 
flutter  analysis  correlation  mentioned  above. 


See  Table  3.1-1  for  tabulated 
coordinates  of  grid 


+down  at 
each  point 


Ficure  3.1-1  TF-8A  Modal  Deflection  Measurenent  Grid 
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TABLE  3.1-1 

-8A  MODE  SHAPE  MEASUREMENT  S' 


INCHES 


STATION 


88.4250 

90.1830 

92.3614 

93.4315 


77.0373 

78.1839 

80.5916 

83.3432 

84.8337 


69.8996 

71.3137 

74.1036 

77.1992 

78.9572 


62.9207 

64.5641 

67.7361 

71.1757 

73.2012 


59.3910 

61.0725 

64.5121 

68.1046 

70.3212 


62.5528 
62.4873 
61 .0042 
59.1665 
58.2638 
59.8087 
59.0027 
57.3262 
55.3273 
54.2311 
56.1493 
55.2466 
53.4089 
51 .2488 
50.2494 
52.6329 
51.6657 
49.6345 
47.3132 
46.0559 
~49.1980 
48.0696 
45.8128 
43.3496 
42.0084 
45.4901 
44.2972 
41.9437 
39.3322 
37.8492 
42.0218 
40.7322 
38.2174 
35.4448 
33.8005 


38.4580 

37.0716 

34.3957 

31.4941 

29.7854 

34.8941 

33.4756 

30.5740 

27.5434 

25.6735 


54.2399 

58.0235 

61.9981 

64.6160 


48.7943 

50.7434 

54.7563 

58.9602 

61.8265 


45.1149 

47.2933 

51.5355 

55.8540 

59.6376 


41.7422 

44.0353 

48.3920 

52.9017 

58.0993 


37.5245 

40.4902 

45.1680 

49.8688 

53.3466 


37.1668 

41.9440 

45.5747 


34.1000 
29.9000 

25.1000 


31.2348 

29.6872 

26.6567 

23.4327 

21.4338 


26.0741 

22.8823 

19.5293 

17.3209 


22.3946 

19.0094 

15.4629 

13.0449 


18.7701 

15.1914 

11.5483 


15.1406 

11.4652 


11.5897 


.7441 


4.1000 

5.0000 

2.5000 


|l 


100 


THIS  PAGE  IS  BEST  QUALITY  PEACIICAflii 
mm  COPY  FURNISHED  TO  DDC _ -- 


Table  3.1-2 


TF-8A  Measured  Modal  Deflections 


Nl)  .  ^ 

1ST  1^1  NO  RENOIf.i 

MUDi 

FRtC  =  5.29 

RAY 

M 

EAY 

H 

1 

0 . 1005  +  0 1 

54 

0.609E-01 

2 

0 .676E+ 00 

55 

0.5825-01 

3 

0  .a37F+00 

56 

0.4175-01 

4 

C  .8925  +  00 

5  7 

0.3735-01 

5 

0.904E+00 

58 

0.3535-01 

6 

0.790E+00 

59 

0.3295-01 

7 

0.774E+00 

60 

0.299E-01 

B 

0  •  766  5  +00 

61 

0. 1375-01 

9 

0.7  755+00 

62 

0.17  )E-01 

10 

0.773E+00 

63 

U . 1615-0  1 

1  1 

0.693E+00 

64 

0. 1345-01 

12 

0.657E+00 

65 

0.1035-01 

13 

0 .65flE  +00 

66 

0.6305- 02 

14 

0.ti£j7E+00 

67 

0.6505-02 

15 

0 . 661 E  +  J  0 

68 

0 . 3  50  5- 02 

16 

0.561E+00 

69 

0.5305-02 

17 

0.5515+00 

70 

u .5705-02 

18 

0.589E+00 

71 

0. 7805-02 

19 

0.5365+00 

72 

0. 7205-02 

20 

0 . 5  74F  +00 

73 

0 .6705-02 

2  1 

) .5055+00 

74 

0 .8505-02 

22 

0.486E+00 

75 

0.9505-02 

23 

0*4885+00 

76 

0 . 8  5  0  5  -  0  2 

24 

0.4335+00 

7  7 

0.9  105-02 

25 

0 .4325+00 

78 

0.9305-02 

26 

0.411E+00 

79 

0.9  735-02 

27 

0.393E+00 

80 

0 . 34  )E-02 

28 

0 , 396E+00 

81 

0. 1115-01 

29 

0.416E+00 

30 

0 . 393E+00 

3  1 

0 . 320E+00 

32 

0 . 3045+00 

33 

C  .2925+00 

34 

0 . 286E  +00 

35 

C  .2885  +  00 

36 

0.2285+00 

37 

0 .2205+00 

38 

0 .2205+00 

39 

0.2075+00 

40 

0 .2  15E+00 

41 

0. 170E  +  00 

42 

0. 160E  +  00 

43 

0. 155E+00 

44 

0. 151E  +  00 

45 

0 . 1475+00 

46 

0 .1095+00 

47 

0. 1045+00 

48 

0 . 1035+00 

49 

0  .  IGIE+OO 

50 

0.9625-01 

51 

0.7135-01 

52 

0 .6545-01 

5'' 

0  .65  15-01 

■it  FRA  LI  7  =  0  = 

.00484 
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Table  3.1-3 

TF-8A  Measured  Modal  Deflections 


MlJOt  INIO.  2 
WAY  H 

1  0.lJ0  5<-0i 

2  0.686E+-0U 

3  O*O33Et0i0 

5  g,T}4J?*00 

7  it  *iJU 

8  D  *  1 25^E  * I]  iJ 

9 


2MO  WING  LENDING 

F»iEO  =  18.08  itNE^^ALIZEO  '^ASS^ 

H 


10 

0.47fiE 

63 

1 1 

^3.  142; 

Oa 

64 

12 

0.  t  SaEvOi] 

65 

13 

0 . 180E  4-00 

66 

14 

) .2225 

i+OO 

67 

15 

0 .2295 

■  +0.0 

68 

16 

—0 .6945 

-01 

6  9 

17 

-0.575E 

i-01 

70 

Id 

-0.301- 

rUl 

71 

19 

-a 

72 

20 

^  Hi 

r 

73 

2  1 

-0,;  jT- 

'tQiTi 

74 

22 

19  ^ 

r+OO 

75 

23 

f>g .) 

76 

24 

25 

*  Lh  4  L  ji  Ll  5  ►'j  0 

-o*  n^ii*oo 

77 

78 

26 

-3.  iU|+3J 

79 

27 

-ii.2 

:  *C0 

80 

28 

-o.?70|*jn 

81 

29 

3o 

=  «0  J 

3  1 
32 

33 

.□a 

34 

*a 

+0-3 

35 

36 

37 

'4:m\ 

+  30 
+  00 

38 

+  33 

39 

-0. 5flii 

+  t!3 

40 

-0 .2  *J0 

■tl 

-0. 1  J5e*00 

42 

43 

-.i-Aiinii 

44 

5*30 

45 

-0.24  7i 

^•3(1 

4o 

-0.247: 

f +00 

47 

-0 . 236i 

:  +00 

43 

-0.2221 

E+OO 

49 

-0.207jE+00 

50 

-0 . 198E+00 

5  1 

-0 . 184E4-00 

52 

-0 . 179i 

E  +  00 

5’ 

-0 . 1  75i 

:+oo 

54  -0.l60E«-0O 

55  -0, 154t<-U0 

56  -0.138E<-00 

57  -0.132E+00 

58  -0.124E+00 

59  -0.112E+C0 
oO  -0. loot *00 

61  -().854E-01 

62  -0.d52E-‘0l 

63  -0,735e-01 

64  -0.fc72E-01 

65  -0.S43E-01 

66  -0.5A1E-01 

67  -0.497E-J1 

68  -0.4o7E~ai 

69  --0.  i  92  E-01 

70  -0.3425-01 

71  -0.305E-01 

72  -o.2ao*-oi 

73  -0.269t-01 

74  -0.2  19  E-0  1 

75  -0.1535-01 

76  -0.18QE-01 

77  -0.1635-01 

78  -0.1505-01 

79  -0.1475-01 

80  -0.140E-01 

81  -0.130E-01 


.00268 
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Table  3.1-4 

TF-8A  Measured  Modal  Deflections 


MODE 

BAY 

1 

2 

3 

4 

5 

6 
7 
U 
9 

10 
1  1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
2  1 
22 

23 

24 

25 

26 

27 

28 

29 

30 
3  I 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 
5'’ 


NO.  3 
H 

0.  lOOt^-J  I 
-0.1  75E  +  01 
-0.5766>’00 
0.771E4-00 
0. 145E<-0  1 
-0.223E+01 
-0  p  I  4  4£  ♦  J  I 
-0,40 5? +00 
y*lJU*0 
a.tiiE+ol 

-0  p  fh  +0 
0,1 IfcC*  Ji 
0. 179E+0  1 
-0. I  82 E +01 
-O.lilE+Ol 


-0 

a* 


46E+00 

25C+0 


'v.m 

0 


-0 . iUbZ 

IE*  Jl 

J  ,  i  *  0 1 

-0  » i 3oc  +0 
-0 .d45E+00 
0. 179E  +  00 
0. 131E+01 
0. 185E+01 
-0, 109E+01 
-0.642E+0U 
0 .248E+00 
0. 123E+01 
0. 176t+01 
-0  ..338E+00 
-0 ,475t+00 
0.276E+00 
0. 1C7E+01 
0. 152E+01 
—0 .62  7E+00 
-0 . 360E+00 
0.259E+00 
0 .906E+00 
0. 127E+01 
-0 .462E+U0 
-J . 261c  +00 

3.201=+J0 
0.709t+00 
0.102E+J1 
-0.3  7ir;  +  00 
-0 . 195E  +  U0 
0 . 146 E  +0  0 


1ST  WING  TORSION 
FPEO  =  37.10 

HAY  H 

54  0.521C+00 

55  0.736E+00 

56  -0.265E+U0 

57  -0.157E+00 

58  0.8575-01 

59  0.344E+C0 

60  0.508E+00 

61  -0.1795+00 

62  -0. 1175+00 

63  0.337E-01 

64  0.190E+00 

65  0.286E+00 

66  -0. 1  185+00 

6/  -0.88aE-01 

d8  0.1575-01 
09  0.784E-01 

70  0.5805-01 

71  -0.940E-01 

72  -0.8245-01 

73  -0.335E-01 

74  -0. 1775-01 

75  0.275E-01 

76  -0.8425-01 

77  -0.527E-01 

78  -0.320E-01 

79  -0.955E-01 

80  -0.118E+00 

81  -0.142E+00 


Gf  NERALI  ZED 

.02735 


MASS= 
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Table  3.1-5 


TF- 

■8A  Measured  Modal  Deflections 

3Pn  WING  BENDING 

MJOt 

NU.  4 

FPEO  =  42.  80  GENERALIZED  MASS 

84Y 

H 

BAY 

1 

0 . 1  joe+0 1 

54 

0.163E»00  .00165 

2 

0.621E4-00 

55 

0.140E+00 

3 

0.539E+00 

56 

0.241E+00 

‘t 

0.591E+  >0 

57 

0.208Ef00 

5 

0.5  72E4-  )0 

58 

0.1 73E+00 

h 

0.267E*00 

59 

0.138E+00 

7 

0 . 255E  fOO 

60 

0.115E*GO 

8 

0 . 2  2  3  E  ^0  0 

61 

0.168E4-00 

9 

0. 185E+U0 

62 

0. 141E4-00 

10 

0. 167F.4-00 

63 

0.1  15E4-0  0 

11 

-0.  lOOE-01 

64 

0.915E-01 

12 

-0.303F-01 

65 

0.688E-01 

13 

-0.  J27t-0l 

66 

0.931E-01 

14 

-0, 140E-t-U0 

67 

0.333E-01 

15 

-0 . 18  IE  <-00 

68 

0. 708E-01 

16 

-0 . 166E  +  00 

69 

0. 546E-01 

17 

-0. 139C+00 

70 

0. 5005-01 

18 

-0. 250E+J0 

71 

0.440E-0  1 

19 

-0. 335E4-00 

72 

0.432E-01 

20 

-0. 3  76E4-00 

73 

0.365E-01 

21 

-0, 136E<-J0 

74 

0.2525-01 

22 

-0.211E+00 

75 

0. 149E-01 

2  3 

-J.292E+0J 

76 

0. 171E-01 

2h 

-0 . 3a9E«-J0 

77 

0, 140E-01 

25 

-0  .440E<-00 

78 

0.950E-02 

26 

-0 .  lOOE^OO 

79 

0.5 JOE-02 

27 

-0 . 1  30E  ♦•OO 

80 

0.1 305-02 

28 

-0.200E<-J0 

81 

0.340E-02 

29 

-0.294t^C0 

30 

-0  .  34JE<-00 

31 

0. 309F-01 

32 

-0 . 3  lOE-02 

33 

-0 .69  7E-01 

34 

-0 . 149E4-00 

35 

-0 . 204E+J0 

36 

J  .  1  63  E  ♦•0  0 

37 

0. 131L+00 

38 

0.595E-01 

39 

-0.11  lE-01 

40 

-0.497E-01 

41 

0.259E+-JU 

42 

0.223E«-00 

43 

0.  159E<-00 

44 

0 . lOOE+00 

45 

0 .6  76E-0  1 

46 

0.298E+00 

47 

0. 269E^0 J 

40 

0.2  1  7E<-J0 

49 

0. 158E+00 

50 

0 . 133En)0 

5  1 

0.288E+0J 

52 

0.249E  +  .)0 

5’ 

0 .2  J2E+00 

Table  3.1-6  ‘ 


TF-8A  Measured  Modal  Deflections 


MOUc 

HAV 

1 

2 

3 

4 

5 

6 
7 
6 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 


NJ.  5 
H 

O.lOOEi-Ol 
0.525E<-00 
0.478E+00 
0.444E4-0U 
0.418E^OO 
0.164E+00 
0.123E4-00 
0.5616-01 
-0.506E-01 
-0.979E-01 
-0.642E-01 
-O.IOIE+JO 
-0.194E+00 
-0.321E+00 
-0.391E+00 
-0. 109E+00 
-0. 139E*00 
-0.2  18E4-00 
-0.329E+C0 
-0.388E+00 
-0.162E-01 
-0.402E-01 
-0.9086-01 
-0. 140C+-0U 
-0.1796+00 
0.9636-01 
0.8166-01 
0 .6026-0  1 
0.5116-01 
0.4896-01 
0. 146E  + JO 
0.139E«-J0 
0. 154E+0J 
0.131E+O0 
0.207E+00 
0.948E-01 
0.  lOOE+OU 
0.  IhIE  +  JO 
0. 183E+00 
0.232E+C0 
-0.240E-01 
0. 140E-02 
0.604E-01 
0.123E+00 
0. 167E+00 
-0. 150E+J0 
-0 . 1076+00 
-0.304C-01 
0.319E-01 
0.731E-01 
-0.221E+00 
-0.173E  +00 
-0.1J5E+00 


4TH  WING  BENDING 
EREO  =  73.40 

BAY  H 

54  -0.398E-01 

55  -O.lOOE-02 

56  -0.219E+00 

57  -0.194E+00 

58  -O.IJIE+OO 

59  -0.730E-01 

60  -0.3806-01 
61  -0.1886+00 

62  -0.164E+00 

63  -0.117E+00 

64  -0.727E-01 

65  -0.397E-01 

66  -0.128E+00 

67  -0.114E+00 

68  -0.8546-01 

69  -0.559E-01 

70  -0.500E-01 

71  -0.715E-01 

72  -0.649E-01 

73  -0.5086-01 

74  -0.307E-01 

75  -0.250E-01 

76  -0.3086-01 

77  -0.303E-01 

78  -0.240E-01 

79  -0.200E-01 

80  -0.170E-01 

81  -0.146E-01 


GENERALIZED  MASS= 

.0008 


THIS  PACE  IS  BEST  QUALITY  PRACnCA^ 
COPY  FURK1SHE2)  TO  DDC 
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Table  3.1-7 


Solf  quality  PRACriCA^HI 

«SPM  COPY  PUR;tJisHf.ii  -  •■ 


TF-8A  Measured  Modal  Deflections 


2NC  WING  TORSION 

MODE 

NO.  6 

FREU  =  74.30 

BAY 

H 

BAY 

H 

1 

J.  lOOFfdl 

54 

-0.499E+00 

2 

-0.649E+00 

55 

-0.704E+00 

3 

-0. 761E-01 

56 

0.996F-01 

4 

0.5806+00 

57 

0.  lcilE-01 

5 

0.942E+J0 

58 

-0. 179E+  JO 

6 

-0 .984E+00 

59 

-0.393E+00 

7 

-0.721E+00 

60 

-0.568E+00 

6 

-0.200E+JO 

61 

0.426E-01 

9 

0.446fc+00 

62 

-0.1206-01 

lU 

0.778E+00 

63 

-0. 1326+00 

11 

-0.865E+ JO 

64 

-0.249E+00 

12 

-0.640E+G0 

65 

-0.337E+00 

13 

-0.216E+00 

66 

0. 177E-01 

1‘* 

0.254E+00 

b7 

0 , 107E-01 

15 

0  .4  76E  +  JO 

68 

-0.732E-01 

16 

-0.451 c+JO 

69 

-0. 132E+00 

17 

-0 . 3566  +00 

70 

-0. 180E+J0 

18 

-0. 131E  +  .J0 

71 

0. 113E-01 

19 

0. 135E+00 

72 

-0.190E-02 

20 

0 .233E+00 

73 

-0.241E-01 

2  1 

-0.283E-01 

74 

-0. 347E-0 1 

22 

-0.1 66E-01 

75 

0.8106-02 

23 

J. 165E-01 

76 

0 . 590E- J2 

24 

0.6l5t;-Jl 

77 

;J.310E-02 

25 

0. 900E-0  1 

78 

0.910E-02 

26 

0.333F+00 

79 

0.910E-02 

27 

0. 255E+U0 

80 

0.9106-02 

28 

0 . 1  1  Ii^+OO 

61 

0.  J 

29 

-0. 257f--0l 

3J 

-0. J77E-01 

31 

0 . 56  7E+00 

32 

0.403‘=  +  J0 

33 

0 . 123E  OU 

34 

-0. 169E+00 

35 

-0 . 334E+JO 

36 

0.602E+00 

37 

0. 397E+  )0 

38 

0.437E-01 

39 

0. 319E+00 

40 

-0.585E+00 

4  1 

0 .408E  +00 

42 

0.295E+ JO 

43 

-0 . 72 lE-01 

44 

-0.460E+00 

45 

-0.725E+00 

46 

0.  J25E+00 

47 

0 . 1  72E+00 

48 

0 . I60t+ JO 

49 

-0 . 525E+CO 

50 

-0. 75  IE  +00 

51 

0 . 1  J4E  +  0a 

52 

0.6d9E-0l 

5’ 

-0. 199C+00 

bfcNFRALI  ZEO  MASS,= 
.0057 
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Table  3.1-8 


TF-8A  Transformed  Modal  Deflections 


MJDE 

ND.  1 

1ST  rtING 
FKEO  = 

BAY 

H 

ALPH 

2 

0.8850E400 

0.1537E-01 

3 

0.77606+00 

0.1402E-01 

4 

0.6690E«-00 

0.12866-01 

5 

0.5690E+00 

0. 1180E-01 

6 

0.4750E4-00 

0.1080E-01 

7 

0.3850E4-00 

0.1014E-01 

8 

0.3010E4-00 

0.90706-02 

9 

0.2280E^00 

0. 76206-02 

10 

0. 1650E^00 

0.66406-02 

11 

0.1070E+00 

0.578CE-02 

12 

0.6100E-0  1 

0.4200E-02 

13 

0.3000E-01 

0.2610E-02 

14 

0, 1200E-01 

0.1170E-02 

1*5 

0.5000E-02 

0.14006-02 

16 

0.4000E-02 

0.10006-03 

17 

0.30006-02 

0.0 

GENERALIZED  MASS 
.00484 
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Table  3.1-9 

TF-8A  Transformed  Modal  Deflections 


mods  no.  2 

BAY  H 

2  0.6500E^00 

3  0.4050E^OO 

4  0.  1700E*00 

5  -0.3600E-01 

6  -0.1900E+00 

7  -0.2700E+00 

8  -0.3100E+00 

9  -0.3070E+00 

10  -0.2780E+00 

11  -0.2300E+00 

12  -0.1750E+00 

13  -0.1180E+00 

14  -0.7500E-01 

15  -0.4500E-01 

16  -0.2700E-01 

17  -0.2000E-01 


2N0  WING 
FfiEO  = 

ALPH 

0.4124E-01 
0.3890E-01 
0.3514E-01 
0.2932E-01 
0.2004E-01 
0. 1199E-01 
0.5930E-02 
0.9900E-03 
-0.2040E-02 
-0.4310E-02 
-0.5500E-02 
-0.4800E-02 
-0.3470E-02 
-0.2220E-02 
-0. lOOCE-02 
-0.7500E-03 


BENDI  NG 
18.08 


CcNc'^ALlZEO  MASS 


.00268 
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Table  3.1-10 


TF-8A  Transformed  Modal  Deflection' 


MOOfc  NCi. 

BAY 


1ST 

FRfO  = 

ALPf^ 


TJr^SION 
37.  lU 


2 

-0.6000E+00 

0.3269E+00 

3 

-0.4500E+00 

0.3127E+00 

-0.2950E>00 

0.2893E400 

5 

-0.1220E*00 

O.2567E+0C 

6 

0.6000E-01 

0.2234E+00 

7 

0. 1850E4-00 

0.1938E*00 

8 

0.2420E+00 

0.1637E+00 

9 

0. 2650E+00 

0.  n41E  +  00 

10 

0.2600E+00 

0.1074E400 

1 1 

0.2230E+00 

O.8476E-01 

12 

0. 1540E>00 

0.6260E-01 

13 

0.8400E-01 

0.4311E-01 

14 

0.3000E-01 

0.2684E-C1 

15 

-0. llOOE-01 

0.1444E-01 

16 

-0.3900E-01 

0.6590E-02 

17 

-0.5400E-01 

0.4230E-02 

GfNEKALIZEO  MASS 


.02735 


Table  3.1-11 

TF-8A  Transformed  Modal  Deflections 


MJD£  NO. 

BAY 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


3RC  WING 
FREO  = 

ALPH 


BENDl  NG 
42.80 


GENERALIZED  "^AS 


0.6000E»00 
0.2250E+00 
-0.8000E-01 
-0.2620E+00 
-0.2900E+00 
-0.2630E+00 
-0. 1750E+00 
-0.2500E-01 
0. 1650E+00 
0.2170E*00 
0.2120E*-00 
0.1780E4-00 
0.  1190E4-00 
0.6800E-01 
0.3000E-01 
0.  1200E-01 


0.4582E-01 
0.3556E-01 
0.1706E-01 
-0.3190E-02 
-0.1713E-01 
-0.2351E-01 
-0.2982E-01 
-0.3423E-01 
-0.2581E-01 
-0. 1157E-01 
-0.4290E-02 
0.8400E-03 
0.3520E-02 
0.344CE-02 
0.2310E-02 
0.186CE-02 


.00165 


A 


Table  3,1-12 

TF-8A  Transformed  Modal  Deflections 


OOt 

NG.  5 

4TH  Wl 
FRF  0 

UY 

H 

ALPH 

2 

0.4850E«^00 

3.4644E-01 

3 

0.5000E-01 

0.2150E-CI 

4 

-0.2080E+00 

-0.q580E-02 

5 

-0.2240E+00 

-0.290SE-01 

6 

-0.9300E-01 

-0.2967E-01 

7 

0.7200E-0  1 

-0.1890E-01 

8 

0. 1630E^00 

-0.1270E-02 

9 

0.1500E+00 

0.1427E-01 

10 

0.6500E-01 

0.2306E-01 

11 

-0.5000E-01 

0.2236E-01 

12 

-O.lllOE+00 

0.1554E-01 

13 

-0. 1310E  +  00 

0.8980E-02 

14 

-0.12O0Ef00 

0.3020E-02 

15 

-0.8600E-01 

-O.lllOE-02 

16 

-0.4800E-01 

-0.192CE-C2 

17 

-0.2900E-01 

-0.2070E-02 

NG 


B  E  ND  I  NG 
73.40 


GENERALIZED  MASS* 


.0008 


Table  3.1-13 

TF-8A  Transformed  Modal  Deflections 


MODE  NO.  ^ 

BAY  H 

2  -0.9000E-0  1 

3  -0.P070E+00 

4  -0.2280E+00 

5  -0.1220E+00 

6  0.2200E-01 

7  0.1270E4-00 

8  0.1370E+00 

9  0.5000E-01 

10  -0.7500E-01 

11  -0.1780E4-00 

12  -0.2110E+00 

13  -0.  19206^00 

14  -0.1300E+00 

15  -0.5500E-01 

16  -O.IOOOE-Ol 

17  0.5000E-02 


2N0  WING 
FKEO  = 

ALPI- 

0. 1952E400 
0. 1618E^OO 
O.IOSAE+OO 
0.4111E-01 
-0.8310E-02 
-0.3421E-01 
-0,4fl4eE-01 
-0.5206E-01 
-0.5045E-01 
-0.4785E-01 
-0.4190E-01 
-0.3576E-01 
-0.2500P-01 
-0.  1470E-01 
-0.4630E-02 
-0.1930E-02 


TORSION  , 

74.30  GENERALIZED  MASS 

.0057 
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TF-BA  Fluttv  Mode)  Vlng  ftigidl 


I  Mass 

O  n 

LB-Sec  LB-$ec 


Inertia 


r 


uiJlr 


44 


«**» 


3.1.2  TF-8A  Static  Aeroelast'ic  Solutions 

Section  2.3  above  contains  the  discussion  of  the  method  used  to  generate 
the  TF-8A  static  aeroelastic  solutions.  Figures  3.1-6  and  3.1-7  show  the 
resultant  twist  distribution  for  the  five  Mach  numbers  under  consideration. 
NASA  supplied  an  optically  measured  static  shape  for  several  cases  and 
Figure  3.1-8  shows  a  comparison  of  these  measurements  with  a  calculated 
shape  for  the  same  lift  and  Mach  number.  Good  agreement  is  seen  to  exist. 

The  use  of  these  static  shapes  is  discussed  above  in  Sections  2.4  and  2.5, 
and  below  in  Section  3.1.3. 


The  Doublet  Lattice  Idealization  of  the  TF-8A  wing  is  shown  in  Subsection 
2.1.  The  theoretical  results  for  the  spanwise  variation  of  section  lift 
coefficient  slope  and  aerodynamic  center  are  presented  in  Figures  3.1-13 


through  3.1-17  for  five  freestream  Mach  Numbers.  The  experimental  data  are 


also  shown.  As  expected,  the  theory  underestimates  the  lift  because  several 


transonic  affects  are  not  accounted  for. 


The  data  were  taken  from  Reference  2  and  corrected  for  aeroelastic 
effects  as  described  in  Subsections  2.1  and  2.4.1.  Limitations  of  the  data 
‘  due  to  insufficient  Reynolds  Numbers  './ere  also  described  in  Subsection  2.4.1. 
It  was  concluded  there  that  the  stall  characteristics  of  the  Ci  and  Cp^ 
versus  a  were  not  reliable.  Thus,  the  values  of  cj,^  were  read  only  in  their 
'non-stalled'  ranae  and  not  necessarily  at  the  section  angles-of-attack  to 
which  the  flutter  model  had  deformed  statically  during  the  flutter  test. 


For  example.  Figures  3.1-9  through  -12  present  the  aeroelastical ly 
corrected  pressure  model  data  (c^,  c^^  vs  fof  the  =0.99  case.  It  is 

obvious  that  the  values  of  c^^  and  especially  c^^  depend  on  the  section 
angle-of-attack.  The  section  angle-of-attack  distribution  across  the  span 
of  the  flutter  model  in  turn  is  the  sum  of  the  rigid  body  angle-of-attack, 
the  built-in  twist  and  the  elastic  twist  as  described  in  Section  2.4.  The 
distribution  of  section  angle-of-attack  for  the  TF-8A  wino  is  given  in 
Subsection  3.1.2.  The  effective  angle-of-attack  at  which  the  ci^  and  c^^ 
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are  to  be  read  is  the  difference  between  the  section  angle-of-attack  of  the 
flutter  model,  and  the  built-in  twist  of  the  pressure  model 

That  is, 


a 


(f) 

eff 


s 


a  (0  . 


where  is  given  in  Table  3.3.2.  The  reason  why  the  elastic  twist  of  the 
pressure  model  is  not  mentioned  here  is  because  it  has  already  been  accounted 
for  in  the  term  i.e.,  the  pressure  data  has  already  been  corrected  back 
to  the  rigid  state. 


The  values  at  which  the  lift  and  moment  curve  slopes  are  to  be  read, 

*^eff’  plotted  in  Figures  3.1-9  through  -12.  In  Figure  3.1-9,  the 

effective  section  angles-of-attack  of  the  flutter  model  lie  in  the  non  viscous 

range  for  the  c^  curve  and  these  are  the  values  read.  In  Figure  3.1-10,  the 

effective  angles-of-attack  of  the  flutter  model  fall  in  the  viscous  range 

only  for  the  section  at  y/{b/2)  «  0.933.  Since  the  curve  is  unreliable 

there,  the  value  of  c,  at  a  ..  *  0  is  used, 

*  a  eff 

The  curves  of  moment  coefficient  versus  effective  angle-of-attack  are 
nonlinear  everywhere  due  to  transonic  and  viscous  effects.  For  this  case 
the  slopes  of  c^^^  selected  for  the  0.99  case  are  shown  in  3.1-11  and 
3.1-12  as  tangent  lines  to  the  curves.  Due  to  the  lack  of  data  points, 
the  actual  shape  of  the  c^^^  vs  a  curve  is  open  to  various  interpretations. 

This  is  unfortunate  because  c„  scorns  to  have  a  large  effect  on  flutter 
speed. 


In  Figure  3.1-12,  three  separate  interpretations  of  the  c^^  versus  a  , 
curve  have  been  drawn  for  Stations  y/(b/2)  =  0.653  and  0.804.  The  tangent 
lines  or  slopes  of  these  curves  are  also  indicated.  The  solid  line  indicates 
our  best  guess  to  the  slope;  however,  the  other  slopes  are  also  retained  and 
used  in  the  flutter  analysis  to  show  the  effect  of  flutter  speed  on  the  slope 

c  .  Two  separate  interpretations  have  been  drawn  for  the  curve  at  y/(b/2)  = 

■"a 

0.48.  The  final  three  curves  of  aerodynamic  center  resulting  from  these 
various  interpretations  are  given  in  Figure  3.1-17. 

Special  consideration  is  given  to  the  curve  of  c^  versus  a  for  the  wing 
tip  Station  y/(b/2)  »  0.933.  As  indicated  previously,  the  point  at  which 
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Cl  and  c  are  to  be  measured  lies  in  the  viscous  region  for  the  pressure 
'ct  m(^  . 

model.  Due  to  the  Reynolds  Number  difference  between  the  pressure  and^ 

flutter  models^  (1.96  x  10^  for  the  pressure  model  and  3.0  -  9.0  x  10  for 

l^he  flutter  model),  the  pressure  model  curves  can  not  neally  be  trusted  in 

the  viscous  region. 

,»The  handling  of  this  difficulty  posed  a  problem  since  the  c^^  versus  « 
curve  is  very  nonlinear,  even  in  the  non-viscous  regions.  The  motions  of 
shock  waves  tend  to  change  the  character  of  the  c^^  curves  drastically  as  t 
passes  from  negative  to  positive  values.  One  possible  solution  is  to  extrap¬ 
olate  the  c  versus  a  curve,  based  on  its  shape  close  to  a  of  zero,  to  the 
point  in  question  ignoring  the  break  in  the  curve  at  '<  *  -1.5  degrees.  When 
this  is  done,  it  is  found  that  the  aerodynamic  center  lay  very  far  aft  on 
the  airfoil,  specifically,  at  /r  ^0.70.  This  result  is  considered 
unreasonable  since  the  aerodynamic  center  usually  moves  forward  at  the  wing 
tip,  not  aft.  If  this  result  /c  =  0.7)  is  used  in  the  flutter  analysis, 

the  resulting  flutter  speed  is  also  unreasonable;  lying  well  above  the 
results  for  the  unweighted  theory.  A  two-dimensional  transonic  solution  was 
applied  to  the  airfoil  section  (plus  boundary  layer)  at  y/(b/2)  =  0.933. 

The  section  data  gave  a  moment  curve  with  a  minimum  point  at  a  =  0  similar 
to  the  data.  However,  the  slopes  were  much  less  resulting  in  an  aerodynamic 
center,  at  a  =  -1.82,  of  x^  /c  =  0.3966.  Three-dimensional  effects  will 
move  this  point  forward.  It  is  felt,  then,  that  the  aerodynamic  center 
selected  for  this  station,  shown  in  Figure  3.1-17,  is  correct. 

A  similar  problem  was  encountered  for  the  case  of  Moo=  0.95.  Figure 
3.1-16  presents  the  unreasonable  aerodynamic  center  distribution  as  a 
dashed  line.  This  curve  was  flattened  to  the  solid  line.  Again,  if  the 
unreasonable  curve  of  aerodynamic  center  is  used,  the  flutter  speed  is  well 
above  the  unweighted  theory  result. 

The  final  values  of  Cf  and  x^  ^  used  for  all  Mach  Numbers  considered 
are  presented  in  Figures  3.1-13  through  -17  where  they  are  compared  with  the 
theory.  The  ratios  of  experimental  to  theoretical  c, ^  and  c^^  are  the  steady 
weight  factors,  as  described  in  Subsection  2.4.  The  curves  presented  in 
Figures  3.1-13  through  -17  were  interpolated  for  c^  ^  and  x^^^  at  the  bay 
positions  for  the  theoretical  as  well  as  experimental  data.  The  weight 
factors  for  the  five  Mach  Numbers  considered  are  given  in  Table  3.1-14. 
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-10  TF-8A  C.  vs.  M_  =  .99;  y/{b/2)  =  .653 


y/(b/2)  \ 

0.653 

0.304  N. 

0.933  > 

Tangent  Slope  (best  estimate) 
Tangent  Slope  (alternate  est.) 

Flutter  Model  aitl 


Fiqure  3.1-12  TF-8A  c„  vs.  a 


99;  y/(b/2)  =  .653,  .804,  .933 


PERCEN 


TABLE  3.1-14 


TF-8A  Steady  Weight  Factors 


M  =  0.6 

00 

M  =  0.8 

00 

Y/(b/2) 

WTSCL 

WTSCM 

Y/(b/2) 

WTSCL 

WTSCM 

.095 

1.141 

1.475 

.095 

1.227 

1.437 

.167 

1.107 

1.440 

.167 

1.238 

1.614 

.224 

1.070 

1.277 

.224 

1.165 

1.462 

.280 

1.075 

1.224 

.280 

1.146 

1.379 

.337 

1.053 

1.223 

.337 

1.145 

1.388 

.394 

1.033 

1.162 

.394 

1.134 

1.298 

.452 

1.034 

1.136 

.452 

1.106 

1.219 

.509 

1.057 

1.133 

.509 

1.146 

1.179 

.567 

1.123 

1.173 

.567 

1.185 

1 .163 

.624 

1.119 

1.205 

.624 

1.215 

1.118 

.680 

1.256 

1.232 

.680 

1.229 

1.108 

.737 

1.285 

1.228 

.737 

1.213 

1.093 

.795 

1.295 

1.181 

.795 

1.208 

1.068 

.852 

1.266 

1.162 

.852 

1.211 

.986 

.909 

1.218 

1.141 

.909 

1.26? 

.894 

.969 

1.105 

1.118 

.969 

1.378 

.689 

M  =  0.90 

00 

M  -  0.95 

00 

Y/(b/2) 

WTSCL 

WTSCM 

Y/(b/2) 

WTSCL 

vrrscM 

.095 

1.216 

1.336 

.095 

.862 

.342 

.167 

1.213 

1.552 

.167 

1.175 

.754 

.224 

1.161 

1.518 

.224 

1.271 

1.051 

.280 

1.151 

1.399 

.280 

1.332 

1.312 

.337 

1.203 

1.456 

.337 

1.383 

1.389 

.394 

1.211 

1.359 

.394 

1.396 

1.333 

.452 

1.250 

1.331 

.452 

1.433 

1.373 

.509 

1.279 

1.285 

.504 

1.443 

1.422 

.567 

1.340 

1.278 

.567 

1.436 

1.442 

.624 

1.385 

1.254 

.624 

1.424 

1.387 

.680 

1.438 

1.261 

.68 

1.430 

1.397 

.737 

1.479 

1.226 

.737 

1.448 

1.395 

.795 

1.501 

1.187 

.795 

1.496 

1.413 

.852 

1.495 

1.044 

.052 

1.485 

1.337 

.909 

1.436 

.808 

.909 

1.401 

1.145 

.969 

1.185 

.450 

.969 

.953 

.658 
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TABLE 

3.1-14  (Cont’d) 

TF-8A  Steady  Weight  Factors 

Y/(b/2) 

WTSCL 

WTSCM 

.095 

.902 

.553 

.167 

1.104 

.902 

.224 

1.157 

1.035 

.280 

1.195 

1.071 

.387 

1.256 

1.045 

.394 

1.288 

.917 

.452 

1.357 

.913 

.509 

1.445 

.949 

.567 

1.590 

1.026 

.624 

1.729 

1.107 

.68 

1.796 

1.133 

.737 

1.758 

1.103 

.795 

1.662 

1.102 

.852 

1.576 

.949 

.909 

1.493 

.878 

.969 

1.279 

.643 

-15 


1 


ight  Factors 


3 


y/(b/2)  =  .452 

WTUSCMfa 

K 

REAL 

IMAGINARY 

*^30 

1 

0 

0 

.968 

-.027 

.1 

.920 

.2 

.847 

-.081 

.333 

.827 

+  .037 

1.0 

y/(b/2)  =  .509 

REAL 

IMAGINARY 

*^30 

1 

0 

0 

.970 

-.025 

.1 

.926 

-.064 

.2 

,854 

-.082 

.333 

.817 

+  .035 

1.0 

y/(b/2)  =  .567 

WTUSCMea 

REAL 

IMAGINARY 

*^30 

1 

0 

0 

.971 

-.022 

.1 

.933 

-.060 

.2 

.863 

-.082 

.333 

.808 

_ 

+  .031 

1.0 

L 

TABLE  .T.l-IB  (Cont'd) 
TF-8A  Unsteady  Weight  Factors 


BAY  NO.  5 


M  =  .90 


y/(b/2)  »  .795 


REAL 

IMAGINARY 

1 

0 

.983 

-.017 

.959 

O 

• 

1 

.907 

1 

o 

o 

.856 

+  .020 

0 

-.012 
-.045 
-.077 
+  .004 


0 

.1 

.2 

.333 

1.0 


BAY  NO.  4 


y/(b/2)  =  .852 


MAG I NARY 


0 

-.014 

-.065 

-.117 

-.001 


IMAGINARY 


0 

-.008 

-.072 

-.139 

-.048 


0 

.1 

.2 

.333 

1.0 


BAY  NO.  3 


y/(b/2)  =  .909 


IMAGINARY 


0 

-.008 

-.057 

-.108 

-.023 


REAL 


1 

1.049 

1.012 

.935 

.626 


IMAGINARY 


0 

0.00 

-.062 

-.129 

-.077 


'3D 


0 

.1 

.2 

.333 

1.0 
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TABLE  3.1-15  (Cont'd) 


TF-8A  Unsteady  Weight  Factors 
.90 


BAY  NO.  2 

y/(b/2)  =  .964 

K 

WTUSCL 

WTUSCMpA 

RFAI 

TMARTNARY 

RFAI  . 

3D 

0 

1 

0 

1 

0 

.052 

.974 

-.004 

1.051 

1 

.105 

.951 

-.050 

1.019 

.2 

.175 

CO 

cn 

o 

• 

1 

.954 

-.116 

.333 

.520 

.661 

1 

. 

o 

.639 

-.109 

1.0 

TABLE  3.1-16 
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TABLE  3.1-16  (Cont'd) 
TF-8A  Unsteady  Weight  Factors 


.95 


BAY 

NO.  8 

y/(b/2)  =  .624  I 

K 

1  WTUSCL 

WTUSCMfa 

REAL 

IMAGINARY 

■dlHi 

1  IMAGINARY 

*^30 

1 

0 

1 

0 

-.044 

.968 

.1 

-.105 

.920 

.2 

-.134 

.823 

-.123 

.333 

+.049 

.707 

+  .041 

1.0 

BAY  NO.  7 

y/(b/2)  =  .68 

|( 

WTUSCL 

¥rU3 

™ _ 

'*20 

REAL 

1  IMAGINARY 

REAL 

IMAGINARY 

*^30 

0 

1 

1 

0 

0 

.076 

.969 

.970 

-.021 

.1 

.152 

.922 

-.100 

.927 

-.074 

.2 

.250 

.792 

-.137 

.840 

-.122 

.333 

.760 

.766 

+  .052 

.699 

+  .033 

1.0 

BAY  NO.  6 

y/(b/2)  =  .737 

K 

WTUSCL 

^2U 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

3D 

0 

1 

0 

1 

0 

0 

.071 

.971 

-.036 

.971 

BH 

.1 

.930 

-.093 

.935 

.2 

■I 

.802 

-.137 

.848 

-.119 

.333 

.710 

.750 

+.050 

.692 

+  .020 

1.0 
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TABLE  3.1-16  (Cont'd) 


TF-8A  Unsteady  Weight  Factors 
M  =  .95 

oo 


BAY 

NO.  5 

y/(b/2)  =  .7 

95 

1  WTUSCL 

UBm 

REAL 

real 

IMAGINARY 

■IBI 

0 

1 

0 

1 

0 

0 

.067 

.973 

-.030 

.973 

-.013 

.1 

.133 

.936 

-.088 

.942 

-.063 

.2 

.220 

.824 

-.136 

.865 

-.112 

.333 

.670 

.734 

+  .044 

.690 

+  .008 

1.0 

BA 

Y  NO.  3 

y/(b/2)  =  . 

909 

WT 

IJSCL 

SCMea 

*^20 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

■Gal 

0 

1 

0 

1 

0 

0 

.057 

.970 

+  .020 

.982 

+  .022 

.1 

.114 

.946 

-.018 

.971 

+  .011 

.2 

.190 

.888 

-.076 

.957 

-.013 

.333 

.570 

.641 

-.092 

.844 

-.054 

1.0 
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TABLE  3.1-16  (Cont'd) 
TF-8A  Unsteady  Weight  Factors 

M  =  .95 

00 


BAY  NO.  2  y/(b/2)  =  .969 


WTI 

wtuscmea 

*^20 

REAL 

■■■PV  iV  M 

IRBHI 

0 

1 

0 

1 

0 

.052 

.971 

+  .022 

.983 

+  .022 

.105 

.951 

-.010 

.972 

+  .013 

.175 

.900 

-.064 

.961 

-.008 

.520 

.658 

00 

o 

• 

1 

.853 

-.054 

TABLE  3.1-17 

TF-8A  Unsteady  Weight  Factors 
.99 


BAY  NO.  11  y/(b/2)  =  .452 


k, 

•^20 

1  WTUSCL 

k^ 

•^30 

RFAI 

inii^^iiiimn 

iiiRTVirni 

0 

1 

0 

1 

0 

0 

.095 

.977 

0.00 

.992 

.1 

.190 

.920 

-.070 

.978 

.2 

.314 

.822 

-.122 

.959 

.333 

.950 

.675 

+  .025 

.892 

-.024 

1.0 

BAY 

NO.  10 

y/(b/2)  =  .509 

k^ 

L  WTUSCL 

WTUSCMea 

KM 

*^20 

REAL 

IMAGINARY 

REAL 

1  IMAGINARY 

0 

1 

0 

1 

■■ 

0 

.090 

.979 

+  .002 

.992 

.1 

.180 

.930 

-.062 

.980 

.2 

.300 

.831 

-.120 

.961 

.333 

.900 

.668 

+  .013 

.897 

1.0 

BAY 

NO.  9 

y/(b/2)  =  .567 

k„ 

1  WTUSCL 

WTUSCMea 

k^ 

*^20 

REAL 

IMAGINARY 

REAL 

IMAGINARY 

*^30 

0 

1 

0 

1 

0 

0 

.085 

.980 

+  .006 

.993 

+  .011 

.1 

.170 

.938 

-.054 

.982 

.2 

.280 

.846 

-.115 

.964 

.333 

.850 

.664 

0.00 

.903 

-.030 

1.0 

TABLE  3.1-17  (Cont'd) 
TF-8A  Unsteady  Weight  Factors 


M.=  .99 


BAY  NO.  2  y/(b/2)  =  .969 


k,_ 

WTUSCL 

*^20 

REAL 

IMAGINARY 

mmmm 

•^30 

0 

1 

0 

-1 

0 

0 

.052 

.976 

0.00 

.986 

.1 

.105 

.952 

-.064 

.974 

.2 

.175 

.876 

-.140 

.942 

.333 

.520 

.569 

-.058 

.793 

-.066 

1.0 
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The  unsteady  weight  factors  are  described  in  Subsection  3.3.  Figures 

3.3-5  through  -10  present  unsteady  weight  factors  for  three  spanwise  stations 

as  a  function  of  section  reduced  frequency,  where  =  -i-  k^  .  The 

>^20  *^20  c  ’"30 

results  at  these  three  stations  are  applied  to  all  wing  bays  as  described 
in  Subsection  3.3.  Tables  3.1-15  through  -17  present  the  unsteady  correction 
factors  for  the  TF-8A  wing.  There  are  no  tables  for  =  0.6,  0.8  since  the 
unsteady  weight  factors  for  these  Mach  Numbers  are  unity  (no  shock  wave 
motion).  In  the  tables  for  the  other  Mach  Numbers  some  inboard  bays  are 
missing.  Again  the  reason  for  this  is  the  fact  that  the  weight  factors  there 
are  unity  (no  shock  wave  motion).  These  data  were  applied  at  the  four 
reduced  frequencies  of  k^^^  =0,  .1,  .2,  .333,  and  1.0.  The  total  weight 
factor  for  any  bay  is  the  product  of  the  steady  and  unsteady  factors: 

UTCL  =  WTSCL*WTUSCL 
^^■CM  =  WTSCM*WTIJSCM 

Normally,  these  totals  are  not  formed  directly.  The  procedure  followed 
is  to:  First,  modify  the  theoretical  AIC  matrix  with  the  steady  weight 
factors,  store,  and  use  them  as  required;  Second,  modify  the  modified  AIC 
matrix  with  the  unsteady  weight  factors,  store  again,  and  use  as  required. 

3.1.4  TF-8A  Flutter  Results  -  Steady  Weighting 

Flutter  analyses  were  run  using  the  "steady  weighted"  AIC's.  The 
discussion  of  the  weinhting  procedures,  along  with  the  unexpected  difficulties 
experienced  in  connection  with  Reynolds  Number  discrepancy  between  the 
pressure  and  flutter  models,  appears  above  in  Subsections  2.1,  2.4,  and 
3.1.3.  In  those  sections  it  was  pointed  out  that  at  high  Mach  numbers  the 
static  aenoelastic  twist  distributions  and  steady  aerodynamic  lift  curve 
slope  data  seemed  to  indicate  that  the  last  few  wing  tip  stations  were 
stalled;  e.g.,  c^^  =0.  A  flutter  analysis  was  run  at  =  0.99  under  this 
assumption.  The  flutter  dynamic  pressure  was  not  reasonable  (qp=  195  psf). 

Since  the  measured  flutter  speed  trends  indicated  no  tip  stalling, 
flutter  analysis  was  rerun  using  AIC's  weighted  under  the  assumption 
that  the  actual  lift  slopes  at  the  static  twist  angles-of-attack  are 

f 

the  same  as  shown  in  the  pressure  model  data  at  smaller  pre-stall ed 
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values  of  angle-of-attack.  The  basis  for  this  assumption  is  the  Reynold's 
Number  discrepancy  between  the  two  models  (See  Section  2.4  above  for 
discussion). 

Figure  3.1-18  shows  a  summary  of  computed  flutter  speeds  under  this  assump¬ 
tion  of  no  tip  stalling  for  the  value  of  C|^  .  Unweighted  theoretical  results 
(using  subsonic  Doublett  Lattice  Method)  are  shown  for  comparison. 


As  discussed  above  in  Section  3.1.3.  the  available  aerodynamic  data  was  not 
adequate  to  unambiguously  define  the  section  moment  slopes,  c^^^.  At  M^ 

0.99,  Figure  3.1-18  shows  the  calculated  flutter  speed  for  three  different 
interpretations  of  c^  ,  and  illustrates  the  sensitivity  of  flutter  speed 
to  c  .  Section  3.1^3  contains  a  detailed  discussion  of  this  moment  data 
and  its  interpretations.  At  M„’  .95  and  0.99,  one  additional  possible 
interpretation  of  c^  (namely,  using  the  steepest  negative  slope  between 
a  .  0  and  a  -  -1.5  )'was  used  which  resulted  in  an  unreasonable  aerodynamic 
center  distribution  (70%  chord  at  the  tip)  and  an  unreasonable  flutter 
speed  (near  300  qsf.  for  both  Mach  numbers). 

Section  5  below  discusses  recommended  work  to  resolve  these  ambiguities. 

The  steady  weight  factors  (best  guess  for  c^  )  used  in  the  TF-8A  analysis 
are  shown  in  Table  3.1-14. 

Figure  3.1-18  shows  an  apparent  discrepancy  between  the  calculated  and 
measured  flutter  frequencies  at  high  Mach  number.  Actually,  the  discrepancy 
is  not  severe  because  the  analysis  shows  that  the  flutter  mode  decreases  in 
frequency  very  rapidly  just  before  flutter  onset.  For  example,  the  analysis 
shows  the  frequency  to  be  dropping  at  a  rate  of  2  Hz  per  3  kts  at  flutter. 
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.8 


.9 
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O 

A 


□  CALCULATED  RESULTS 
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"a 


MEASURED  RESULTS 
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STEADY  WEIGHTED  AIC'S  (OTHER  ESTIMATES  OF  c  ) 

ot 
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n 
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.9 


Figure  3.1-18  TF-8A  Flutter  Results 


MACH  NUMBER 


3.2  RESULTS  OF  ANALYSIS  -  PHASE  II 


Flutter  analyses  were  also  run  using  AIC's  reflecting  the  unsteady 
weighting  factors".  Section  2.5  discusses  the  theory  and  Section  3.3 
discusses  the  application  of  that  theory  to  the  TF-8A  wing. 

Tables  3.1-15  through  3.1-17  show  the  unsteady  weighting  factors  for 
reduced  frequencies  of  0,  .1,  .25  and  0.5.  The  "steady  weighted"  AIC's  for 
these  values  of  were  weighted  by  these  factors  and  then  the  AIC's  for  the 
total  list  of  reduced  frequencies  required  by  the  flutter  analysis  was  ob¬ 
tained  by  interpolation  between  these  basic  four. 

The  results  of  the  flutter  analyses  showed  that  for  the  TF-8A  the 
unsteady  weighting  factors  changed  the  flutter  speed  by  less  than  two 
percent  at  all  Mach  numbers. 
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3.3  PHASE  III  RESULTS 

The  theories  developed  under  Phase  III  of  the  contract  are  presented 
in  Subsection  2.5.  It  was  noted  there  that  the  major  transonic  effect  is 
the  perturbation  pressure  induced  by  shock  wave  motion.  The  effects  of 
shock  wave  motion  were  added  to  the  classical  lifting  surface  theory  and 
the  results  discussed  and  correlated  in  Subsection  2. 5. 1.3.  The  correlation 
of  this  procedure  using  the  NLR  7301  airfoil  for  the  range  of  frequencies 
considered  in  the  experiment  is  very  encouraging. 

The  full  transonic  theory  described  in  Subsection  2.5  is  not  entirely 

complete  since  the  procedure  for  determining  the  shock  wave  motion  for 

oscillatory  flow  has  not  been  tried.  The  nondimensional  shock  wave  motion 

function,  A(k^)/A(k^  =  0)  obtained  experimentally  for  the  NLR  airfoil  at 

M  =  0.7  and  a  =  1.75°,  can  temporarily  be  used  to  obtain  the  shock  wave 
00 

motions  for  the  TF-8A  and  YC-15II  wings.  It  was  felt,  however,  that  if 
this  approximation  is  to  be  used,  then  a  less  general  theory  could  be  used. 
Specifically,  the  method  described  in  Subsection  2. 5. 1.2,  as  mentioned 
above,  would  be  the  logical  choice  because  it  contains  the  major  transonic 
effect,  i.e.,  pressures  due  to  shock  wave  motion  and  because  of  its  good 
correlation  with  data.  A  second  reason  for  using  this  method  is  the  effort 
involved.  At  its  current  level  of  development  the  general  theory  would 
require  more  time  than  available  to  run  all  of  the  conditions  for  the  TF-8A 
and  YC-15II  wings. 

The  theory  of  Subsection  2. 5. 1.3  gives  the  incremental  lift  and  moment 
due  to  shock  wave  motion  as: 

S  ■  ''1*2  <‘^p  •  S  >  ^ 

o  S2 

6-1  ,  j  1  upper  surface 

I  -1  lower  surface 

62  ,  j  1  shock  moving  aft  with  increasing  a 

( -1  shock  moving  forward  with  increasing  a 
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Where  Cp 


are  the  steady  pressures  on 


either  side  of  the  shock 


wave,  Xq  is  the  amplitude  of  shock  motion,  and  t  is  the  non-dimensional 

time  (tU  /c)  for  an  acoustic  wave  to  reach  the  shock  wave  from  the  trailing 
60 

edge.  The  oscillatory  value  of  x^  is  approximated  using  the  curve  of 
Xo(kr)/Xo(kr  *  0)  for  the  NLR  7301  airfoil  and  the  value  of  x^(k^  =  0) 
for  the  airfoil  under  consideration  (see  Figure  2.5-33).  Thus,  the 
for  any  wing  is 


Xo(kr)  =  ‘ 


where 

G(kr)  *  5iQ(k^)/XQ(k^  »  0)  for  NLR  7301  airfoil  at  =  0.7,  a  '  1.75°. 


If  the  acoustic  time  7  is  computed  as  in  Reference  4,  then  the  constant 
A  can  be  set  to  1.0.  To  simplify  the  calculation,  the  method  of  Reference  4 
was  adopted  instead  of  using  the  acoustic  pulse  generator  (along  with  the 
estimated  value  for  A).  The  approximate  method  of  Reference  4  consists  of 
using  the  following  formula: 

X  ■  X 

1  f  dx 

c  J 

X  =  c 

where 

\  -  0.7  |m  -  M.]  +  H. 

Where  M  is  the  local  surface  Mach  Number 

The  weight  factor  for  the  lift  coefficient  is  given  in  Subsection  2.^ 
as : 

WTUSCL(k^)  -  ^  (k^) 

a  /  ® 
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Where  the  transonic  i.e.,  (2DTRANS)  given  as  the  total  of 

o  a 

the  classic  theory,  plus  the  transonic  increment,  lc^  . 


The  weight  factor  for  the  lift  then  becomes: 
WTUSCL(k^)  “  (1  +  BCL  G(k^)LCL(l(^)/  (1  +  BCL) 


(3.3-1) 


where 


BCL 


LCL(k^)  = 


(S  -‘^P  )  ’V'r  ■  “>/ ’“A 

S2  / 

^  ('“'(V  •  0)/c, 

a  /  ® 


,  •  0) 


exp  (-i2k^T) 


The  derivation  of  the  weight  factors  for  the  moment  coefficient  about 
the  elastic  axis  follows  and  is  similar. 


‘'ll,  ■  (V'  - 

o'  a 


Where  x  is  the  shock  location.  The  unsteady  weight  factors  for  c_  is  then; 
s 


WTUSCM^^  "  0  ^  l-CM(k^)y  (1  +  BCM) 


(3.3-2) 


where 


(x,/c  -  X  /c) 

Bcm  =  — - S-3-^ —  BCL 

(V4  -  «.Jc) 


LCM(k^)  . 


o  /  o 


exp  (-i2k^T) 


The  values  of  the  steady  flow  classic  theory  are 

^.(2D)(k^  =  0)  =  2ir/6 


i  '  r 
a 


c(20),kr  .  0)  -  (*,/c  -  V4) 
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The  weight  factors  just  described  are  calculated  for  three  sections 
along  the  wing  and  applied  to  all  of  the  aerodynamic  bays.  Figure  3.3-1 
shows  the  bays  that  are  associated  with  the  three  sections  for  the  TF-8A 
wing  while  Figure  3.3-2  shows  it  for  the  YC-15II  wing. 

The  elastic  axis  data,  ^  /c.  to  be  used  in  the  formulas  for  the 
unsteady  weight  factors  are  found  in  Table  3.1-1.  The  values  of  c^^  (k^) 

and  c  (k  )  are  obtained  from  the  two  dimensional  compressible“subsonic 

lifting  surface  theory. 


The  shock  wave  data  to  be  used  in  these  formulas  could  not  be  taken 
directly  from  steady  wind  tunnel  data  because  of  the  elasticity  of  both  the 
pressure  model  and  the  flutter  model  and  because  of  the  difference  in  built- 
in  twist  between  the  two  models.  Shock  wave  location  and  movement  are 
highly  dependent  on  the  local  section  angle-of-attack  and  boundary  layer 
build-up.  Since  both  the  built-in  twist  and  elastic  twist  of  the  pressure 
model  are  different  than  those  of  the  flutter  model,  the  following  procedure 
was  followed. 


The  first  step  is  to  find  the  total  section  angle-of-attack , a 


(f) 


of 


each  of  the  three  wing  sections  on  the  flutter  model  at  an  estimated  dynamic 
pressure  for  flutter.  The  section  angle-of-attack,  is  composed  of 


rigid  body  angle  a 

Jf) 


elastic  twist  and  built-in  or  'jig'  twist, 


where  'f  indicates  flutter  model. 


a(f)  _  a(f)  ^  €(f)  +  c(f) 
s  j  e 


From  this  value  of  section  angle-of-attack  for  the  flutter  model,  is 
subtracted  the  value  of  for  the  pressure  model. 

J 

La  (eff)  _  a(f)  _  e(p) 
s  s  j 
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Selected  Wing  Sections 


>n  of  Aerodynamic  Bays  with  Selected 
for  the  YC-15II  Wing. 


TABLE  3.3-1 


Elastic  Axis  Locations 


TF-8A  YC-15 


Y/(b/2) 

^ea/c 

Y/(b/2) 

^ea/c 

.309 

.5 

.16 

.38 

.657 

.5 

.42 

.43 

.933 

.5 

.85 

.5 

TABLE  3.3-2  -  TF-8A  Incremental  Angles-of-Attack  Between 

the  Elastic  Flutter  Model  and  Pressure  Model 


y/(b/2) 

Aa 

5 

M«»  0.8 

M«=  0.9 

M„  =  0.95 

M<^=  0.99 

.309 

0.2° 

0 

0.6° 

1.2° 

.657 

-0.4° 

-0.7° 

0 

.933 

-1.62° 

-0.8° 

-2.4° 

-1.8° 

TABLE  3.3-3  -  YC-15II  Incremental  Angles-of-Attack  Between 

the  Elastic  Aircraft  and  Pressure  Model 


y/(b/2) 

Aoj  M  »  0.76 

.163 

O 

o 

• 

.42 

-1.0° 

.85 

-2.0° 
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The  incremental  angle-of-attack,  ,  is  the  angle-of-attack  on  the 

pressuv'-e  model  required  to  match  the  flutter  model  conditions.  (This  angle 
is  an  effective  angle  and  must  be  reduced  to  an  actual  angle  by  using  the 
aeroelastic  correction  of  the  pressure  model  in  reverse.  This  process  was 
not  done  in  the  calculations;  however,  the  resulting  error  is  small.) 

The  aeroelastic  calculations  of  for  the  two  flutter  models  is 
given  in  Subsection.  3.1.2  for  the  TF-8A  wing  and  3.4.2  for  the  YC-15II 
wing.  The  value  of  for  the  pressure  model  is  given  in  Figure  3.1-6. 

The  values  of  Aa  for  the  TF-8A  wing  at  various  Mach  Numbers  are  presented 
in  Table  3.3-2.  The  values  of  Ao^  for  the  YC-15II  are  given  in  Table  3.3-3. 

The  next  step  in  the  process  is  to  find  experimental  pressure  model 
data  for  pressure  distributions  that  correspond  as  c^ose  as  possible  in  a 
to  the  Aa  tables  above  and  to  match  this  data  to  theoretical  calculations 
using  the^ Jameson  and/or  Garabedian  methods  (see  Reference  16  and  Appendix  D 
for  a  description  of  the  Garabedian  method).  The  inverse  process  produces 
a  modified  airfoil  shape  when  the  experimental  pressure  distribution  and 
Mach  Number  are  used.  The  modified  shape  exhibits  a  thickening  in  the 
vicinity  of  the  trailing  edge  which  is  typical  of  boundary  layer  growth 
(see  Figure  3.3-3  for  examples). 

Early  in  the  contract,  calculations  were  performed  using  the  inverse 
Garabedian  method  on  TF-8A  wing  sections  but  at  angles-of-attack  that  did 
not  match  the  values  of  Aa^  given  in  Table  3.3-2.  In  order  to  take 
advantage  of  this  work,  the  resulting  modified  airfoil  shapes  (shown  in 
Figure  3.3-3  for  M^  =  0.99)  were  used.  Instead  of  revising  the  airfoil 
shape  to  match  the  data,  the  angle-of-attack  and  Mach  Numbers  were  varied 
to  produce  a  best-fit  of  the  theory  (Jameson  Method)  to  the  data.  A  table 
of  these  increments  is  given  in  Table  3.3-4. 

Some  of  the  section  pressure  distributions  of  the  TF-8A  wing  could 
not  be  matched  using  this  procedure  and  for  these  the  inverse  Garabedian 
procedure  was  reapplied  near  the  correct  angle  of  attack.  For  the  case  of 


ji 
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Figure  3.3-3  Comparison  of  Original  Pressure  -  Model  Airfoil  Section  Shapes  With  Those 
Required  to  Match  Experimental  Data  For  The  TF-8A  Wing.  Required  Airfoil 
Generated  by  Garabedian  Inverse  Process. 


the  YC-15II  all  of  the  airfoils  were  generated  using  the  inverse  Garabedian 
at  angles -of-attack  which  were  close  to  the  desired  ones  given  in  Table 
3.3-3.  The  output  of  the  inverse  Garabedian  procedure  also  includes  a 
change  in  angle-of-attack,  «a,  from  that  which  corresponds  to  the  input 
experimental  pressure  distribution.  In  Tables  3.3-4  and  3.3-5,  these  are 
indicated  by  an  asterisk  (*). 

The  term  6a  is  now  added  to  the  values  of  Aa^  and  is  added  to  the 
two-dimensional  value  of  Mach  Number  (M,,p  =  M^  cos  A). 


AOg  *  AOg  +  6a 


^20  “  ^2D 


The  values  of  Aa^  and  Hgp  input  into  the  Jameson  steady 

transonic  method  and  the  necessary  shock  data  calculated. 

The  shock  data  required  for  the  calculation  of  the  weight  factors  are: 

(1)  The  jump  in  pressure  across  the  steady  shock,  C  -  C  ;  (2)  The  shock 

>^s^  *^s, 

£.  I 

location,  x^/c,  (3),  the  local  Mach  Number  distribution  (so  that  t  may  be 
calculated);  and  (4)  The  change  in  shock  location  with  respect  to  angle-of- 
attack,  3AQ(k^  =  0)/3a.  The  last  term  is  done  in  a  numerical  manner  by 
varying  the  angle-of-attack  slightly  about  the  mean  position.  Thus, 

_ 0  _  A^ 

9a  ■  Aa 

An  example  Is  presented  in  Figure  3.3-4.  The  values  of  interest  for 
the  TF-8A  and  YC-15II  wings  are  given  in  Tables  3.3-6  and  3.3-7. 

The  quantities  in  these  tables  are  used  to  produce  the  desired 
incremental  unsteady  weight  factors.  These  factors  are  plotted  for  the 
TF-8A  wing  in  Figures  3.3-5  through  3.1-10.  These  figures  are  for  free 
stream  Mach  Numbers,  M  *  0.9,  0.95  and  0.99.  The  results  for  M^  =  .8  are 
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TABLE  3.3-4  -  (TF-8A)  Corrections  to  Jameson 

Transonic  Theory  to  Give  Best 
Match  to  Data 


y/(b/2) 

M  *  C 

_ 

.B 

4M 

M  «  0 

m 

6a 

.9 

£M 

M,=  0 

6a 

.95 

6M 

R 

.99 

6M 

0.309 

0.657 

0.933 

■ 

■m 

1 

0.58® 

2.08® 

0 

0.009 

B 

0 

0 

HQ 

TABLE  3.3-5  -  (YC-15II)  Differences  Between  Data 

and  Garabedian  Theory  for  a  and  M 


y/(b/2) 

M.- 

6a 

.76 

6M 

.163 

-0.85®* 

0* 

.42 

-1.5®* 

0* 

.85 

-0.18®* 

0* 

*  Could  not  obtain  good  match  by  changing  a  and  M  in  the  Jameson  Method 
and  had  to  perform  Garabedian  inverse. 
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Figure  3.3-4  Steady  Pressure  Distribution  Near  the  Shock  Wave  for 

Angles  of  Attack  Near  the  Required  Value  (l.S")  for  the 
TF-8A  Wing  (y/(b/2)  =  0.657,  =  M3P  =  0.99,  =  .80). 

Illustration  of  the  Calculation  9X/9a  (C^r  =  Cpr  7- 


TABLE  3.3-6 


Shock  Location  &  Motion 
Variables  for  the  TF-8A  Wing  Section 


"3D 

(MJ 

^20 

y/(b/2) 

.80 

- 

.309 

.80 

- 

.657 

.80 

- 

.933 

.90 

- 

.309 

.90 

.711 

.657 

.90 

.72 

.933 

.95 

.309 

.95 

.75 

.657 

.95 

.75 

.933 

.99 

- 

.309 

.99 

.80 

.657 

.99 

.77 

.933 

^2 

Xs/c 

0 

Shocks 

• 

No  Shock 

-1 

0.1351 

-1 

0.2141 

- 

No  Shock 

-1 

0.2373 

-1 

0.3894 

- 

No  Shock 

-1 

0.4445 

-1 

0.3894 

$2  Sj 


3.154 

1.766 

1.912 

2.685 


1.2900 

3.1515 


2.9223 

4.7413 

3.8964 

4.8705 


Table  3.3-7 


Shock  Location  &  Motion 

Variables  for  the  YC-I5II  Wing  Sections 


S2  Sj 


*^20 

y/(b/2) 

62 

x^/c 

T 

Bl 

.163 

1 

1 

.3644 

3.29 

mm 

.42 

1 

1 

.4122 

3.35 

.759 

.42 

-1 

-1 

.4350 

2.44 

.761 

.85 

-1 

-1 

.4644 

2.07 

1.776  Upper  surface 
1.031  Lower  surface 
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unity  because  there  are  no  shock  waves  at  this  Mach  Number,  and  thus,  are 
not  plotted.  The  results  for  y/(b/2)  -  0.309  for  all  Mach  numbers  are  not 
shown  for  the  same  reasons;  i.e.,  no  distinct  shock  exists  for  this  inboard 
wing  section. 

Figures  3.1-11  through  3.1-14  present  the  weight  factors  for  the  three 
sections  of  the  YC-15II  wing  plotted  versus  reduced  frequency  for  one  Mach 
Number  (M  =  0.76).  The  section  at  y/(b/2)  =  0.42  requires  two  plots  because 
there  are  shock  waves  on  both  top  and  bottom  surfaces. 

Notice  that  the  incremental  unsteady  weight  factors  are  unity  at  k^  =  0. 
These  factors  are  used  to  produce  a  weight  factor  correction  to  the  steady 
weight  factors  previously  discussed.  As  stated  in  Subsection  2.4,  the  total 
weight  factor  is  the  product  of  the  steady  weight  factor  and  the  unsteady 
weight  factor. 

The  incremental  unsteady  weight  factors  are  plotted  versus  a  two- 
dimensional  sectional  reduced  frequency.  The  two-dimensional  sectional 
reduced  frequency  is  related  to  the  three-dimensional  reduced  frequency 
as  follows  (see  Subsection  2.4): 

kr  *  k 
•^20  *^30  c 

The  values  of  k^  for  the  TF-8A  wing  are: 

"^30 

k  (TF-8A)  =  0,  0.1  ,  0.333,  1.0 
^3D 


The  values  of  k  for  the  YC-15II  wing  are: 
^30 

k  (YC-15II)  =  0,  0.1 ,  0.25.  0.5 
''30 


Si 


163 


164 


O  WTUSCL 

B  WTUSCMea 

—REAL 

— IMAGINARY 
LOWER  SURFACE 


.92 

WTUSCL 

WTUSCMea 

REAL 


WTUSCL 
02  WTUSCMea 
IilAGINARY 


Figure  3.3-5  Unsteady  Weight  Factors  for  the  TF-8A  Wing 
(y/(b/2)  =  0.657,  =  0.90,  M.n  =  0.71 
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Figure  3.3-6  Unsteady  Weight  Factors  for  the  TF-8A  Winq 
(y/(b/2)  =  0.933,  Mjp  =  0.9,  '  °-72) 
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Figure  3.3-7  Unsteady  Weight  Factors  for  the  TF-8A  Wing 
(y/(b/2)  =  0.657,  =  0.95,  =  0.75) 
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Figure  3.3-8  Unsteady  Weight  Factors  for  the  TF-8A  Win 
(y/(b/2)  =  0.933,  Man  =  0.95,  Mon  =  0.75 


Figure  3.3-9  Unsteady  Weight  Factors  for  the  TF-8A  Wing 
(y/(b/2)  =  0.657,  M3D  =  0.99,  M2D  =  0.80). 
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Figure  3.3-10  Unsteady  Weight  Factors  for  the  TF-8A  Wing  (y/(b/2)  =  0.933 
M3P  =  0.99,  Mgjj  =  0.77) 
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Figure  3.3-11  Unsteady  Weight  Factors  for  the  YC-15II 
(y/(b/2)  =  0.16,  M20  =  0.759) 
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Figure  3.3-12  Unsteady  Weight  Factors  for  the  YC-15II 

(y/(b/2)  =  0.42,  M«n  =  0.759),  Upper  Surface 
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Figure  3.3-13  Unsteady  Weight  Factors  for  the  YC-15II 

(y/(b/2)  =  0.42,  =  0,759),  Lower  Surface 
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3.4  RESULTS  OF  ANALYSIS  -  PHASE  IV 

The  methods  of  Phases  I  -  III  were  applied  to  the  YC-15II  aircraft. 

3.4.1  YC-15II  Basic  Data 

Reference  17  documents  the  YC-15II  basic  data. 


3.4.2  YC-15II  Static  Aeroelastic  Solutions 

The  methods  of  Section  2.3  were  applied  to  generate  the  YC-15II  wing 

static  aeroelastic  shape  at  the  critical  wing  Mach  Number,  M  =  0.76.  Figure 

00 

3.4-1  shows  the  results  for  a  typical  payload/fuel  configuration.  As  can  be 
seen  the  unswept,  relatively  stiff  wing  undergoes  quite  small  deflections  even 
at  this  high  dynamic  pressure.  The  maximum  magnitude  of  the  wing  tip  angle- 
of-attack  is  about  2.5  degrees  (minus). 

3.4.3  Aerodynamic  Data  and  Weighting  Factors  for  the  YC-15II  Wing 

The  Doublet  Lattice  Method  idealization  is  presented  in  Subsection  2.1. 

The  resulting  theoretical  distribution  of  section  lift  coefficient  slope  and 
aerodynamic  center  are  presented  in  Figure  3.4-4.  The  theoretical  distributions 
are  compared  to  those  obtained  by  wind  tunnel  tests. 


Figures  3.4-2  and  3  show  that  the  lift  and  moment  curve  slopes  are 
functions  of  the  section  angle-of-attack.  Specifically,  the  angle-of-attack 
distribution  required  is  one  that  is  relative  to  the  pressure  model  at  zero 
rigid  body  angle-of-attack.  The  distribution  of  section  anqle-of-attack  for 
the  YC-15II  aircraft  consists  of  (1)  the  rigid  body  angle,  ,  (2)  the 


built-in  twist,  ,  and  (3)  the  elastic  twist  calculated  at  a  given 

dynamic  pressure  corresponding  to  a  flutter  speed  previously  estimated.  The 
superscript  (a)  refers  to  the  aircraft.  This  curve  is  shown  in  Figure  3.4-1 
of  Section  3.4.2.  The  angle-of-attack  distribution  of  the  pressure  model  at 
zero  rigid  body  angle-of-attack  is  just  the  pressure  model  built-in  twist,  j^is 
is  also  shown  in  Figure  3.4-1.  The  desired  section  angle-of-attack  is  then 
the  difference  between  the  aircraft  and  pressure  model  distributions,  i.e., 


.(a) 


+  e  +  e  since  e- 

j  e  j  c  j 


(a)  _  ^  (p) 
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This  difference,  Acx  which  is  also  plotted  in  Figure  3.4-1,  is  the 
final  result.  The  points  at  which  c^^^  and  cm^  are  to  be  calculated  are  now 
known  and  are  indicated  in  Figures  3.4-2  and  -3.  The  resulting  Cs  and  x 

01  O  •  w  • 

distributions  are  plotted  in  Figure  3.4-4.  The  theoretical  distribution, 
also  shown  in  the  figure,  was  used  as  an  aid  in  finding  the  experimental 
distribution  because  only  three  spanwise  points  were  available. 

The  ratio  of  experimental  to  theoretical  values  of  Cj^  and  c^  at  the 

various  bay  locations  are  the  steady  weight  factors.  Table  3.4-1  presents 

the  results  of  these  calculations  for  two  estimates  of  Cm  . 

'"a 

These  two  estimates  come  from  the  slopes  on  either  side  of  the 

points.  These  two  values  of  are  carried  through  the  flutter  calculation, 

"a 

The  unsteady  weight  factors  are  described  in  Subsection  3.3  and  plots 
of  them  versus  the  reduced  frequency  are  given  in  Figures  3.3-11  through 
-14  for  three  spanwise  stations.  The  data  for  these  three  spanwise  stations 
are  applied  to  all  of  the  bays  in  a  manner  also  shown  in  Subsection  3.3. 

The  values  of  the  two-dimensional  reduced  frequency,  l^r2Q’  related  to 
the  actual  values,  3S  shown  in  Subsection  2.4  as,  k^^p  "  *^^3D’ 

The  reduced  frequencies  considered  for  the  generation  of  the  AIC's  are 
kj,  =  0,  0.1,  0.25,  0.5.  The  tables  of  the  unsteady  weight  factors  used, 
WTUSCL  and  WTUSCMg;^  (about  elastic  axis)  are  given  in  Table  3.4-2.  There  is 
a  slight  disparity  between  the  tables  and  curves  for  the  inboard  half  of  the 
wing.  Also,  the  shock  motion  on  the  lower  surface  at  y/(b/2)  =  0.42  was 
omitted,  see  Figure  3.3-13.  The  effects  of  this  omission  and  disparity 
should  be  very  small  as  far  as  flutter  is  concerned  since  the  flutter  speed 
changes  due  to  the  unsteady  weight  factors  is  itself  small  and  due  to  the 
fact  that  only  the  inboard  portion  of  the  wing  was  affected. 
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3.4.4  YC-15II  Flutter  Analysis  Results 

Figure  3.4-2  shows  the  measured  (high  speed  wind  tunnel  model)  sectional 
lift  curve  coefficient  versus  angle-of-attack  for  three  span  stations.  As  can 
be  seen,  the  statically  deformed  wing  is  still  in  the  essentially  linear  range 
of  the  aerodynamics.  Figure  3.4-4  shows  a  plot  of  sectional  c^^  vs  percent 
semispan  using  the  approximately  linear  portion  of  the  slopes  from  the  previous 
figure  (from  between  angles-of-attack  of  zero  and  -2.5  degrees).  The  AIC's 
were  weighted  to  match  Figure  3.4-4  for  the  "steady  weighting"  exercise.  This 
is  the  same  weighting  that  was  used  in  the  production  YC-15II  flutter  analysis, 
which  resulted  in  the  previously  reported  minimum  calculated  flutter  speed 

3t  M  =  0.76.  Flutter  analysis  results  are  discussed  in  detail  in  Reference 

00 

17.  Thus,  the  static  aeroelastic  twist  effects,  which  were  ignored  during 
the  original  YC-15II  flutter  analysis,  are  indeed  negligible. 

Table  3.4-2  shows  the  "unsteady  correction"  multiplicative  weighting 
factors  for  the  YC-15II  wing.  These  were  included  in  the  YC-15II  flutter 
analysis  resulting  in  a  flutter  speed  reduction  of  some  2-3  KEAS.  Hence, 
the  YC-15II  results  are  similar  to  the  TF-8A  results  in  that  the  unsteady 
transonic  corrections  are  virtually  negligible  in  effect  compared  to  the  steady 
corrections.  The  steady  corrections  were  of  course  very  important  on  the  YC-15II 
as  they  resulted  in  a  calculated  flutter  speed  equal  to  89  percent  of  the  flutter 
speed  predicted  by  purely  theoretical  AIC's  at  •=  0.76.  The  two  values  or 

c  used  had  little  influence  since  they  are  close  to  each  other  on  the  out- 

m 

boSrd  wing.  The  final  value  of  flutter  speed  was  taken  as  the  average  of  that 
obtained  using  the  two  values  of  c^^  . 
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Figure  3.4-1  YC-15II  Static  Twist  Distribution, 


min 


1 

i 
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YC-15II  Steady  Weight  Factors 
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TABLt  3.4-2  (Cont'd) 
YC-15II  Unsteady  Weight  Factors 
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TABLE  3.4-2  (Cont'd) 

YC-15II  Unsteady  Weight  Factors 
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TABLE  3.4-2  (Cont'd) 
YC-15II  Unsteady  Weight  Factors 
M  =  .761 


TABLE  3.4-2  (Cont'd) 
YC-15II  Unsteady  Weight  Factors 
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TABLE  3.4-2  (Cont'd) 
YC-15II  Unsteady  Weight  Factors 
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3.5  FLUTTER  MODEL  TEST  PLAN 


A  flutter  model  test  is  presented  to  support  the  analyses  in  determining 
the  effects  of  some  of  the  critical  parameters  on  flutter  of  supercritical 
wings.  A  brief  flutter  analysis  of  the  proposed  flutter  model  was  done  to 
verify  that  the  desired  flutter  speed  range  could  be  achieved  with  a  reason¬ 
able  stiffness  distribution. 

3.5.1  Description  of  Model 

Design  and  build  a  cantilevered  wing  flutter  model  representing  a  typical 
STOL  transport  wing.  The  model  is  to  be  tested  in  the  NASA  Langley  Transonic 
Dynamics  Tunnel.  The  model  will  have  the  following  features: 

1.  Interchangeable  airfoil  segments. 

a.  Supercritical  shape 

b.  Conventional  shape 

2.  Variable  angle  of  attack  at  the  wing  root. 

3.  Variable  wing  twist  distribution. 

4.  Instrumented  to  provide  spanwise  internal  load  distribution. 

5.  Removable  engine  weights. 

The  model  will  be  a  single  spar  wing  with  twelve  balsa  segments  represent¬ 
ing  the  airfoil  shape(s).  Two  dummy  engines  representing  the  important  mass 
properties  win  be  mounted  to  the  wing  using  rigid  pylons  to  eliminate  tuning 
effects  of  the  engines.  These  dummy  engines  will  also  be  removed  so  as  to 
gather  data  for  a  higher  frequency  flutter  case.  The  spar  will  be  instrumented 
with  a  series  of  strain  gages  along  the  span  to  measure  bending  in  two 
directions  as  well  as  torsion.  Figure  3.5-1  shows  the  general  arrangement 
of  the  model . 

The  variations  in  angle  of  attack  of  the  wing  will  be  attainable  by  using 
either  the  tunnel  floor  turntable  or  the  wide  wall  mounting  fixture  in  the 
NASA  Langley  16  ft  TDT. 
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The  variations  of  wing  twist  distribution  will  be  made  by  shimming  between 
the  wing  segments  and  the  spar. 

3.5.2  Wind  Tunnel  Test 

The  wind  tunnel  facility  anticipated  for  testing  the  cantilevered  wing  is 
the  NASA  Langley  16  ft  TOT  which  uses  Freon  for  a  test  medium.  Normal  testing 
method  of  changing  Mach  No.  with  total  pressure  constant  will  be  used.  Tunnel 
speeds  will  be  increased  until  flutter  onset  is  obtained  or  some  predetermined 
maximum  is  reached  due  to  tunnel  operating  limit  or  model  applicability. 

3.5.3  Model  Configurations 

The  two  basic  configurations  will  be  with  dummy  engines  and  without  dummy 
engines.  The  three  major  parameters  for  each  configuration  will  be  the  airfoil 
shape,  wing  root  angle  of  attack,  and  wing  twist  distribution. 

3.5.4  Wind  Tunnel  Test  Data 

The  model  test  data  for  any  configuration  will  consist  of  the  following: 

1.  Zero  speed  vibration  modes  and  damping. 

2.  Wind  tunnel  parameters. 

3.  Flutter  speed  and  frequency 

4.  Wing  spanwise  internal  load  distribution. 

5.  Wing  aeroelastic  deformation  if  measureable  with  optical  devices. 
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CANTILEVERED  WING  FLUHER  MODEL 
FOR  NASA  LANGLEY  16  FT.  TDT 


REMOVABLE 
SOLID  BALSA 
SEGMENTS 


MOUNTING 


ENGINE  INSTALLATION 


Figure  3.5-1  Proposed  Cantilevered  Wing  Flutter  Model 


SECTION  IV  CONCLUSIONS 


4.1  CONCLUSIONS  RELATED  TO  THE  FLUHER  ANALYSIS 

1.  The  theoretical  subsonic  method  seriously  underestimates  the  relevant 
steady  aerodynamic  data  (lift  curve  slopes  and  aerodynamic  center  distribution) 
at  the  critical  Mach  number  for  a  wing  with  "supercritical"  airfoil  sections. 
This  is  in  contrast  to  recent  results  for  subsonic  transport  designs  with 
conventional  airfoils  where  the  theory  has  been  in  good  agreement  at  the 
critical  Mach  number  (and  overestimated  the  derivatives  at  low  Mach  numbers). 
For  example,  the  YC-15II  wing  steady  lift  curve  slope  at  the  critical  Mach 
number  (0.76)  was  17  percent  higher  than  predicted  by  the  Doublet  Lattice 
Method,  while  for  the  TF-8A  the  wing  Cj^  test  data  was  approximately  double 
the  Doublet  Lattice  prediction  (at  M  =  0.99).  Figure  4.1-1  below  summarizes 
the  general  trends  for  conventional  and  supercritical  airfoils  of  recent  ex¬ 
perience. 


2.  The  sectional  lift  curve  slopes  for  a  supercritical  airfoil  possess 
an  inherently  nonlinear  character  over  the  pre-stall  range.  However,  it  appears 
that  the  pre-stall  range  can  be  approximated  by  two  slopes;  a  small  slope  over  a 
small  angle  of  attack  around  zero,  and  a  steeper  slope  at  higher  angles  of  attack 
(but  still  before  viscous  effects  are  seen).  Figure  4.1-2  shows  typical  data. 
This  characteristic  "dual  slope"  phenomenon  has  shown  up  consistantly  on  three 
different  supercritical  designs  and  is  more  pronounced  the  more  "supercritical" 
the  design.  The  steeper  slope  is  a  transonic  effect,  and  cannot  be  predicted 
by  subsonic  lattice  methods. 

/: 
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3.  The  "steady  corrections",  wherein  the  AIC  for  =  0  was  weighted  to 
match  the  steady  test  aerodynamic  data,  and  then  these  weight  factors  applied 
to  all  reduced  frequencies,  resulted  In  flutter  speeds  which  were  In  reasonable 
agreement  with  measured  flutter  speeds.  Analysis  using  purely  theoretical 
AIC's  resulted  In  significantly  non-conservative  flutter  speed  predictions, 
for  both  the  TF-8A  flutter  model  and  the  YC-15II  prototype. 

4.  The  unsteady  transonic  weight  factors  had  a  negligible  effect  on  the 
flutter  speed  of  both  the  TF-8A  and  YC-15II.  This  may  be  because  the  reduced 
frequency  at  flutter  is  relatively  low  for  both  aircraft;  =  0.15  for  the  TF- 
8A  flutter  model  and  k^  *  0.125  for  the  YC-15II.  Since  the'^di  recti  on  of  the 
effect  was  detrimental  §n  both  aircraft  (even  though  negligible  in  magnitude), 
these  techniques  should  be  applied  to  a  configuration  with  high  flutter  freq¬ 
uencies  before  any  final  conclusions  are  made. 

5.  For  advance  design  configurations  which  use  supercritical  wings  and 
for  which  no  high  speed  aerodynamic  data  exists,  application  of  the  steady 
three-dimensional  Jameson  method  (reference  11)  or  similar  method  should  be 
made.  This  method  has  been  implemented  Into  a  production  status  at  Douglas 
and  much  correlation  work  done.  Steady  and  unsteady  corrections  could  then 
be  based  on  this  predicted  transonic  aerodynamic  data  In  lieu  of  test  data. 
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6.  The  effect  of  static  twist  on  flutter  is  obviously  significant  since 
the  final  section  angle  of  attack  may  end  up  in  the  smaller  slope  range,  the 
steeper  slope  range,  or  even  the  viscous  effect  range.  It  seems  clear,  however, 
that  the  certification  analysis  must  use  the  steepest  slopes  regardless  of  the 
Igcruise  twist  distribution,  so  as  to  cover  the  maneuvering  range  of  the  air¬ 
craft,  and  other  unknowns. 

Accounting  of  static  twist  effects  are  obviously  very  important  in 
a  correlation  analysis  of  a  flutter  model  with  a  relatively  flexible  wing. 

7.  Incorrect  Reynolds  number  simulation  in  pressure  models  and  flutter 
models  may  lead  to  errors: 

a)  If  the  high  speed  flutter  model  does  not  fly  at  a  high  enough 
Reynolds  number  the  wing  tips  may  stall  prematurely  resulting  in 
non-conservative  flutter  speeds  predictions  or  in  buffet  or  stall 
flutter  condition  which  are  not  realistic. 

b)  If  the  pressure  model  and  flutter  model  were  flown  at  different 
Reynolds  numbers  (as  in  the  TF-8A  case)  they  will  experience  tip 
stall  at  different  angles  of  attack,  confusing  the  correlation 
analysis. 

c)  Pressure  models  may  not  be  "rigid".  The  TF-8A  pressure  model  was 
significantly  flexible.  In  the  future,  the  flexibility  of  the 
pressure  modal  should  be  calculated,  and  if  significant,  should 
either  be  re-designed  out  of  the  model  or  else  the  pressure  data 
from  the  test  should  be  properly  modified.  See  Section  2.1  for  a 
discussion  of  this  problem  and  its  ramifications. 

8.  Sufficient  section  lift  and  moment  data  should  be  available  vs  a  to 
unambiguously  define  the  curves  so  that  accurate  c^^  and  c^  data  can  be 
obtained  at  the  proper  twist  angles.  This  is  especially  trSe  of  the  moment  data 
since  it  is  so  highly  nonlinear.  In  practice,  this  means  getting  the  data  request 
to  the  relevant  Aerodynamics  personnel  in  a  timely  manner,  since  left  to  their 
own,  they  may  not  get  the  required  pressure  data  at  a  fine  enough  a  grid  or 

at  large  enough  negative  angles  of  attack. 
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4.2  CONCLUSIONS  RELATED  TO  TRANSONIC  AERODYNAMICS 

Three  significant  transonic  effects  on  lifting  surface  aerodynamics 

% 

Nonuniform  flow  fields  modify  the  basic  aerodynamics  such  that 
liftiR?  surface  elements  can  propogate  their  effects  upstream 
thrx)ugh  a  supersonic  zone.  This  seems  to  be  very  important  for 
^supercritical  airfoils  operating  at  their  design  (shock  free) 
condition. 

b)  The  shock  wave  generates  a  flow  field  through  the  action  of  a 
vorticity  distribution  on  the  shock  wave  itself.  The  effects  on 
the  wing  surface  are  especially  large  just  aft  of  the  shock  wave. 

c)  The  movement  of  the  shock  v/ave  with  angle  of  attack  creates  a 
very  large  effect  on  the  perturbation  pressures  and  loads.  This 
load  is  proportional  to  the  shock  movement  times  the  pressure 
jump  across  the  steady  shock  wave.  For  supercritical  airfoils  the 
shock  movement  versus  angle  of  attack  can  be  very  large.  For  air¬ 
foils  with  shock  waves  this  effect  is  probably  the  largest  of  the 
three  effects. 

2.  Transonic  load  effects  are  highly  dependent  on  the  shock  wave 
location  because  much  of  the  loading  is  concentrated  in  the  immediate  vicinity 
of  the  shock  wave.  For  supercritical  wing  sections  the  location  of  the  shock 
wave  is  a  sensitive  function  of  angle  of  attack,  boundry  layer  thickness, 

u 

wind  tunnel  wall  corrections  etc.  Thus  the  static  aeroelastic  twist  dis¬ 
tribution  is  very  important  to  the  flutter  analysis. 


are: 


a) 
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SECTION  V  RECOMMENDATIONS 


5.1  RECOMMENDATIONS:  FLUTTER  ASPECTS 

This  study  confirms  that  supercritical  wings  exhibit  signicantly  lower 
flutter  speeds  than  a  conventional  wing  of  equal  size  and  rigidity.  It  was 
demonstrated,  furthermore,  that  these  lower  flutter  speeds  could  be  accurately 
predicted  if  the  Doublet  Lattice  unsteady  aerodynamics  were  properly  weighted 
by  the  relevant  test  data. 

The  major  problem  encountered  during  this  study  was  that  the  relevant  aero¬ 
dynamic  data  required  to  unambiguously  define  the  important  sectional  lift 
and  moment  slopes  were  not  available.  A  secondary  complication  was  that  the 
static  aeroelastic  twist  of  the  TF-8A  flutter  model  was  extreme.  Finally, 
both  examples  treated  in  this  study  exhibited  flutter  at  low  reduced  frequencies; 
thus,  the  importance  of  the  unsteady  weight  factors  in  general  was  not 
established. 

Future  studies  should  be  planned  so  as  to  address  these  problems.  In 
Section  3.5  a  proposed  cantilevered  wing  flutter  model  is  discussed  that 
should  prove  to  be  a  useful  tool  to  further  evaluate  transonic  effects  on 
flutter  speed.  The  static  twist  distribution  of  this  proposed  wing  will  be 
more  representative  of  current  transport  aircraft  than  was  the  TF-8A  flutter 
model.  The  wing  geometry  is  the  YC-15  prototype,  for  which  adequate  definition 
of  aerodynamic  data  exists.  Furthermore,  analysis  has  shown  that  removal  of  the 
engine  weights  will  change  the  flutter  frequency  from  a  low  value  to  a  high 
value  of  k^,  thus  allowing  a  better  evaluation  of  the  unsteady  weighting  factors 
as  a  function  of  reduced  frequency. 
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5.2  RECOMMENDATIONS  FOR  THE  DEVELOPMENT  OF  TRANSONIC  LIFTING  SURFACE  THEORY 

It  has  been  demonstrated  that,  for  the  configurations  considered  in  this 
report,  the  steady  flow  transonic  effects  are  the  most  important.  The  basic 
reason  for  this  is  the  fact  that  these  configurations  fluttered  at  relatively 
low  reduced  frequencies  where  unsteady  effects  are  small. 

However,  it  is  still  important  to  obtain  the  unsteady  effects  for  two 
reasons:  1)  increased  accuracy  of  flutter  results,  and  2)  application  to 
aircraft  configurations  that  flutter  at  moderate  to  high  reduced  frequencies. 

The  unsteady  effects  may  not  be  important  enough  to  develop  and  employ  very 
expensive  finite  element  methods  but  it  does  seem  worth  the  effort  to  develop 
less  expensive  methods.  An  example  of  such  a  method  is  the  transonic  lifting 
surface  method  discussed  in  Subsection  2.5  of  this  report.  Specific  re¬ 
commendations  to  fu'-ther  the  development  of  this  method  are  outlined  as 
follows: 

1.  Extend  the  method  so  that  the  amplitude,  and  phase,  of  the  unsteady 
shock  wave  motion  can  be  determined. 

2.  Refine  the  shock  motion  interference  flow  representation.  Replace 
the  point  vortex  with  the  proper  distribution  along  the  shock. 

3.  Investigate  the  reason  for  poor  correlation  of  the  method  in  the 
supersonic  zone  and  make  appropriate  corrections.  Pertinent  to  this 
investigation  is  the  evaluation  of  the  Impact  of  the  various  simplifying 
assumptions  made  during  the  analysis. 

4.  Extend  the  method  to  the  three-dimensional  case. 

Transonic  lifting  surface  theories,  such  as  the  one  described  in  this 
report,  and  the  one  in  References  10  and  12,  require  reliable  production  sub- 
and  supersonic  linear  methods.  Subsonic  methods  are  well  developed  and  very 
reliable.  The  same  can  not  be  said  of  supersonic  methods,  however.  Thus  it 
is  recommended  that  supersonic  methods  be  brought  up  to  the  same  status  as 
the  subsonic  ones. 
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APPENDIX  A 


SOUPXE  IN  A  NON-UNIFORM  FLOW  FIELD 


Lifting  surface  theory  is  based  on  the  assumption  that  the  disturbances  to  tie 
fluid  caused  by  a  surface  moving  through  it  are  small.  These  small  disturbances 
travel  through  the  fluid  at  the  speed  of  sound  as  acoustic  waves.  Classic 
lifting  surface  theory  makes  the  further  assumption  that  the  medium  through 
which  the  waves  travel  is  uniform.  It  is  felt  that  one  of  the  major  transonic 
effects  arises  due  to  the  fact  that  acoustic  waves  encounter  non-uniform  flow 
fields.  This  is  especially  Important  for  forward  traveling  waves  in  a  near 
sonic  stream.  Thus,  an  expression  is  derived  here  for  an  acoustic  source  in 

a  non-uniform  flow  field. 


To  start  the  derivation,  consider  a  moving  coordinate  system  xyz  such  that 
the  flow  at  infinity  is  unifonn  and  equal  to  U,  in  the  x  direction.  Consider 

an  acoustic  pulse  emitted  at  n„,  c„  at  time  t  i.  If  the  flow  is 

unifom  near  the  source,  then  the  pulse  will  be  circular.  If,  m  general,  the 
flow  is  slightly  perturbed  then  the  wave  will  be  deformed  (see  Figure  A-1). 


CTlCIPlliQi  PKB  NO^  fUiSD 


For  no  perturbations  in  the  flow  field,  the  potential  is: 


1  6  (  (t  -  t)  -  p/a) 


4ti  °0 


where  =  (x  -  (t  -  t)  +  r^ 


=  acoustic  radius 


r'  -  (y  -  t  ^z  -  z/ 


where  a  is  the  speed  of  sound  and  where  is  the  source  strength.  When  per¬ 
turbations  exist,  in  the  flow  the  potential  is 


-  ,  1 

^  4tt 

p 


where  g(x,  y,z,  t-i)  =  0  is  the  equation  of  the  surface  of  the  wave,  ^  is 
an  equivalent  acoustic  radius  such  that  o/'q  is  the  wave  front  stre-^rth,  and 
a  is  the  effective  wave  speed. 

Now  as  a  point  on  the  lifting  surface  moves  along,  it  creates  sources  of 
varying  strength  at  various  points  along  its  flight  path  at  various  time 
(see  Figure  2.5-1 ). 


To  find  the  total  potential  field  due  to  this  string  of  pulses  simply  sum  or 
integrate  them  along  the  path  as  a  function  of  time.  The  potential  at  time 

i 

t  is  a  sum  of  all  the  sources  generated  previously  along  the  flight  path. 


The  results  for  perturbed  flow  is: 


r\ 
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A  theorem  for  the  proper  integration  of  delta  functions,  6  ,  is; 

N 

f  (t)  6  |g(T)  I  di  =  ^ 

i=l 

where  t.  are  the  real  roots  of  g{T)  =0  lying  between  t  =  a  and  b  and  where 
N  is  the  number  of  these  roots. 

The  term  ag/Sx  is  the  negative  of  ag/at  which  can  be  shown  to  be 


where  ’  n  is  the  velocity  of  the  wave  normal  to  itself.  The  velocity  Vy 
is  composed  of  the  acoustic  propagation  velocity,  an,  and  the  fluid  velocity 

Thus, 


and  thus , 

N 

(A-1) 


where  g(T^ )  =  0  i  =  1  ,N  ' 

The  solution  for  g  (x^)  =  0  i  =  1,  ..,  N  amounts  to  finding  thex^  values  for 
all  of  the  waves  that  pass  over  point  (x,y,z).  Figure  2.5-16  (b)  presents  a 
uniform  flow  supersonic  case.  For  a  point  lying  within  the  Mach  cone  (Xi,z^) 
two  waves  are  felt.  For  points  lying  outside  the  Mach  cone  (X2,Z2)  no  waves 
are  felt.  In  the  uniform  flow  subsonic  case  only  one  wave  passes  a  point  (see 
Figure  2.5-1  6(a)).  In  mixed  non-uniform  transonic  flow,  the  problem  is  com¬ 
plicated  as  figure  2.5-16(c) ,(d)  shows.  In  particular  Figure  2.5-16 (d)  shows 
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that  point  (x^z^)  has  three  waves  passing  over  it  while  (x2,Z2)  has  one 
wave  passing  over  it  even  though  it  lies  outside  the  Mach  cone. 


For  uniform  flow: 


p(t^  )  =  p(t^  )  =  a(t  -  ) 

V  =  T  U 


n  ^  = 


■  ^0  ■  /  p(t^) 

g  (t^)  =  0  =  t  -  -  p(t^)  /  a 


then  p 


1  + 


^  ^ 


P  +  (x  -  Cq) 


U  (t  -  tJ  M 

00  '  T  '  00 


and 


are  found  solving  g  (x^)  =  0. 


a2(t-T^)^  =  p2  «  (x  -  (t  -  t))^  +  r^ 


or  solving  for  a  (t  -  x^)  gives: 
a(t  -  T^)  =  +  [  M^  (x 

where  =  (x  -  +  6^r^ 


So)  t  “ij 


Now,  acoustic  waves  generated  in  the  future,  (t  -  t^)  <0  have  no  effect  at 
time  t,  thus. 


(t  -  Tj)  >0 

For  supersonic  flow  (M  >  1)  both  roots  (+R)  are  admissible  since  M  (x  -  r 

00  —  00  ^  Q 

Also,  for  supersonic  flow  (x  -  1  ®  since  waves  do  not  travel  up  stream. 

2  ^ 

Finally,  since  6  is  negative  the  negative  sign  outside  of  the  braces  must 
be  used. 
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a  (t  -  T.) 


4-  [  M.  (X  -  5„)  *  R]  1=1 

5 

-4-  [H.  (X  -  5„)  -  R]  1=2 


Supersonic  flow 


For  the  subsonic  case  <  1 »  ^ 

solutions  which  fit  the  condition  (t  -  t^)  l  0  but  only  one  makes  sense 

physically,  1 .e. , 


a  (t  -  T,)  =  [-H,  (x  -  c„)  t  R]  1  =  1 


Subsonic  flow 


Knowing  the  roots  allows  us  to  move  on  to  find  the  expression  for  the 
potential  as  follows; 

-  I  V  '  n  I  . 

Using  the  equation,  p  =  a  (t  -  t^)  in  the  expression  for  p  1  1  +  — g—  |gives 


1  +  =  P  + 


(t  -  T^)  as^  +  (x  -  ^0^  I 

Using  the  expression  for  the  roots  gives: 


H  •  n 


(t  -=  T.)  a  (x  - 


Thus,  for  uniform  flow  the  expression  for  the  potential  becomes 


aB^  _ I _ 1  a£  _  ■ 

4ttR 

!  t  ^  ^  (Hjx  -t^)  -  R>( 


Supersonic  source 


Subsonic  source 
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These  are  the  classic  expressions  for  a  source  in  a  uniform  flow.  The 

evaluation  of  the  potential  4)  for  the  case  of  non-uniform  flow  is  dependent 

on  the  ability  to  numerically  calculate  acoustic  waves  propagating  in  a  non- 

uniform  stream.  Such  an  acoustic  wave  generator  can  be  used  to  determine  the 

local  normal  vector  n,  acoustic  radius  p,  and  intercept  time  of  a  wave 

passing  over  the  receiving  point  (x,y,z)  that  originated  at  the  sending  point 

(c  •  n  .  c  ).  Whether  the  flow  is  uniform  or  non-uniform,  the  acoustic  wave 
'0  0  0 

expands  at  the  speed  of  sound.  The  wave  strength  along  the  wave  front 
diminishes  with  increased  time  as  the  wave  expands.  Because  of  non-uniform 
flow  conditions,  the  wave  distorts  from  the  circular  form  and  changes  in 
strength  along  the  periphery  of  the  wave.  One  way  to  express  this  fact  is  to 
say  that  the  wave  dilation  p  is  expressed  in  terms  of  the  uniform  flow  acoustic 
wave  radius  p  and  a  correction  factor  f. 


p  =  pf  =  a(t  -  t)  f 


The  term  p  =  a (t  -  t)  gives  the  increase  in  radius  with  time  while  f  accounts 
for  the  wave  strength  distorsion  along  the  wave  front.  The  expression  for  4> 
then  becomes: 


4ti  p(t^)  f 


ST 

1  +  V  •  n 


(A-3) 


where 

g  (t.)  =  0  i  =  1 ,  N 


(A-4) 
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APPENDIX  B 


FLOW  FIELD  CALCULATIONS  USING  THE  TWO-DIMENSIONAL, 

STEADY,  TRANSONIC  JAMESON  PROCEDURE 

The  two-dimensional  Jameson  procedure  transforms  the  space  around  the  airfoil 
into  a  rectangle,  one  unit  high  and  two  units  wide.  This  rectangle  is  divided 
up  into  a  grid  of  points.  The  modified  differential  equation  is  solved  for 
by  finite-difference  procedures  in  this  rectangular  space.  The  solution  pro¬ 
duces  a  value  of  potential  at  each  point  in  the  rectangular  plane.  The  pro¬ 
cedure  in  finding  the  velocity  at  a  particular  point  in  physical  sapce  is  to: 

1.  Find  the  location  of  the  point  of  interest  in  the  transformed  space. 

2.  Find  the  derivative  of  i|>  in  this  space  at  the  point  of  interest. 

3.  Find  the  velocity  transformation  to  be  applied  to  the  differentials 
of  to  produce  the  velocity  vector  in  the  physical  plane. 

The  following  equations  and  sketches  describe  the  process  of  transforming 
from  the  physical  to  the  computational  plane. 
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A  ■  xp 

B  ■  C.6256  ♦  C 
C  -  39-0.625  6 

'^1  +  * 

D  «  (xp-0.6256  )  /0.375 
6  ■  sign{xp) 


Computational 

Grid 


xc 


zp  /  >/\  ♦  zp* 

(  A  If  IxpU  0.625 
(  B  If  |xp|>  0.625 


TE 


TE 


Figure  B-5  Computational  Plane 


Once  the  desired  points  are  transformed  into  the  computational  plane  the 
desired  derivatives  3({>/  zc  can  be  numerically  obtained  from  the  tabulated 
values  of  vs.  zc  and  xc.  The  following  procedure  Is  required  to  transform 
these  derivatives  in  the  computational  plane  into  velocity  components,  Vx 
In  the  physical  plane. 

=  96(1>/3X0  +  COSa 

=  3A({|/3Z0  +  sina 


where  3A(}i/3xo  and  3A4>/3zo  are: 


3A(J)/3XO  =  3A(}>/3xp  {-^)  -  3A4!/3zp  (-^^) 

3A({i/3ZO  =  3A<|i/3xp  (  )  +  3A((./3Zp  (-^) 

2  2 
where  p  =  xp  +  zp 


and  where  3A<|)/3xp  and  3A4>/3zp  are; 

3A4>/3xp  =  3A!t>/3X^  -  3A4>/3zp  3s(xp)/3Xp 

3A4i/3Zp  »  3A(f./3Z^ 

The  qu  ntlties  3A4)/3xp  and  3A(})/3zp  are  found  as  follows: 

3A(Ji/3)^  =  3A(ti/3XC  A 

~  1/1 
3A<;/3Zp  =  3A<t>/3ZC  ^  /  (j  +  Zp) 

where  A  =  ^  2/2 

1  /  [1  *  (^)^]  .  I  xc  I  >  .625 


A  finally  the  terms  3A(t)/3xc,  3A(j)/3zc  are  determined  numerically  from  the  grid 
of  A<t>  vs.  xc  and  xc. 
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APPENDIX  C 


SUPERSONIC  TWO-DIMENSIONAL  DOUBLET  LATTICE  METHOD 


For  supersonic  flow  the  following  integral  equation  holds: 


w 


3_  ^  I 
8Z  35 


where  I  = 


/ 


x-B(z-;) 
A<^(5)  ^2D 


i.e. 


dc 


and  where  is  the  potential  jump  across  the  airfoil  and  is  the 

potential  due  to  a  source  of  unit  strength,  i.e., 

♦20  =  {-ixM^(x-c))/2B 

where  =  (x-c)^”B^(z-c) 

2  2 

B  =  M  -1 
00 

^  =  u)Mco/UcoB^ 

where  w  is  the  downwash  boundary  condition  applied  to  the  lifting  surface 
and  Jo  (xR)  is  a  Bessel  Function  of  the  first  kind. 

The  first  step  in  the  derivation  of  a  Doublet  Lattice  Method  for 

supersonic  flow  is  the  evaluation  of  the  double  differential  of  the  integral, 
2 

i.e.,  3  I/3Z35. 

Performing  the  first  differential  gives: 

||- --  BA(^(c=x-B(z-5))  A(^2d  (5=x-B(z-5)) 

^X-B(z-5) 

I  A4>(^)  9<|)2o/H  dc 

{ 

i  .e. 
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Now 

A(J)2q  (5=x-B(z-c))  =  exp(-ixM^B(z-^)) 

since  Jo(0)  =  1 .0 


Taking  the  second  derivative  of  I  gives 

8^I/8zac  =  B  aA(j)(5=x-B(z-c))/Sz  ^'^2D  (^~x-B(z-c)) 

+  B  A4i(5=x-B(z-c))  9A<t'2[)(^''x-B(z-;))/3z 

-  B  A(|)(5=x-B(z-c))  9A<J.20/9^  I 

'C=x-R(z-C) 

x-B(z-c) 

+  r  2 

J  A4>(5)3 

C 

i.e. 

The  terms  of  this  expression  are  evaluated  as  follows:  (use 
A<t)  =  A(^(c=X-B(z-i;)),  ^2D  “ 

3A4i/3Z  =  ^  =  a/  (-B) 

3^20/ 3z  =  -i(xMy2)  exp(-iAM^B(z-c)) 


3A(f)20 


/3c  =  -Oi(xR)  X3^  exp  (-ixM^(X-c)) 
3C 


3A(|>20/3c 


2B 

=  -BX  U~0-  exp  (-ixM_^B(z-c)) 
5=x-B(z-c) 


since  lim  (xR)/xR  —  1/2  as  xR  —  0 
^  (-ixM^(X-c)))/2B 

where 

^  ,  2 
P  =  X  -Ji  X3  R 

az  ac  sZac 
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2, 

and  ^  _  B  (z-c)  _ 

8c  “  R  "  3Z 

9CSz  R  r3 

Before  these  terms  are  substituted  into  the  expression  for  3^1/8285  the 
limit  i -►  c  must  be  applied  since  the  receiving  points  (x,z)  must  range  over 
the  airfoil  surface  (which  is  z  =  ?).  In  the  limit  as  z  t; 

R  -  I  x-c  I 

3^2o/9^  -i^M^/2 

lc=x-B(z-0 

F  -*■  -  Ji(x  |x-c| )  xB^/|x  •?! 

A<^(C=X-B(z-c))  -*■  A^(x) 

A(J)"(5=X-B(z-?))  A(^l  (x) 

A<t2D(5  =  X  -  B(z  -d)  ^ 

The  final  expression  for  32l/3z3c,  line  z  -►  c  is: 

3^I/3z3C  -►  — ^  I  A(|i'(x)  +  iXM^  AiJ>(x) 

f  A(|>(5)F(x|x  -c|  )exp(-ixM_^(x-c))clf.  j  (C-1) 

£.e. 

=  J^(,\|x-c|)/x|x-c| 


lim 
z  -►  c 


where  F 
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■-  i- 


Notice  that  c  is  always  less  than  or  equal  to  x  thus  x  -  ^  is  always  positive 
and  the  absolute  value  sign  can  be  dropped.  For  simplicity  of  notation 
introduce 


P 

w 


x(x  -  i)  then 

^  j  A<>"(x)  +  ixM^A(t.{x)  +  X 


A(ji(5)  F(p)exp(-iM  p)  dp 


{C-2) 


where 

A 

F(p)  =  Ji(p)/P.  and  5  =  x  -  p/x  and  (x)  =  8A4./3X 

Equations  (C-l)  or  (C-2)  are  the  final  results  for  the  downwash  due  to  a  dis¬ 
tribution  of  potential  jump  a<^(x),  on  a  lifting  surface  (and  wake).  Lifting 
surface  theories,  such  as  the  Doublet  Lattice  Method  (DLM)  do  not  use  a  dis¬ 
tribution  of  potential  jump  but  use  a  distribution  of  pressure  jump, 

C  (x).  This  distribution  of  pressure  jump  is  simulated  in  the  DLM  by  a  series 
o/point  pressure  doublets.  The  point  pressure  doublet  (shown  in  Figure  C-l) 
possesses  the  following  properties:  (1)  the  pressure  jump  and  potential  jump 
are  zero  for  points  forward  of  the  doublet;  (2)  the  pressure  is  a  delta 
function  at  the  location  of  the  doublet;  and  (3)  the  pressure  jump  is  zero 
downstream  of  the  doublet  but  the  potential  jump  is  not  (wake).  Mathematically, 
this  is  expressed  as  follows: 

0  X  <Xq 

2a^P  Mx  -  Xq)  X  = 

-2(a<p  i  a<p)  =  0  X  >  X.Q 
00 
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z 


Figure  C-1  Potential  Jump  Due  to  a  Point  Pressure 
Doublet 


Solving  the  differential  equation  for  A<j>  downstream  of  and  taking  into 
account  the  properties  at  and  ahead  of  x.^  gives: 

Act)  “  A<t)p  H(x  -  Xq)  exp  (-  iw  (x  -  Xq)  /  U^)  (C-3) 

Placing  (C-3)  Into  the  first  two  terms  of  (C-2)  gives: 

A<t)'+ixM_^A(})  =  A(t>  {{(x-Xq)  +  JfW  ^^(x"XQ)exp(-iu)(x-Xp)/U^)  {C-4) 

r  00^ 

when  it  is  noted  that  3H/ax  =  6.  The  expression  for  Ac{i(5)  given  in  terms 
of  (C-3)  and  the  variable  p  is: 

A(^(C)  +  Act)p  HU  -  x^j)  exp  (-iu)  (c-Xp)  /  U^)  = 

Act)  H(5-x  )  exp  (-ic*)(x-x^)  /  U^)exp  (+  p) 

Combining  the  last  part  of  this  term  with  the  term  exp(-iM^p)  in  (C-2)  gives: 
exp(io)p/xU^)  exp(-iM^p)  =  exp(-  ip/M^) 
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Placing  this  term  and  (C-4)  into  (C-2)  gives  the  downwash,  Wp,  for  a  point 
pressure  doublet. 


w 


-B 


+  iWh(x-x.)  exp{-ia)(x-x„)  /  U„) 


p  -  T  "^p 


'o'  u  ”0' 

00 

*^p..e. 


0'  '  oo' 


+  xexp(-ia3(x-XQ)  /  U^)  J"  H(£.-XQ)e”^P^*^«>F(p)dp 


The  step  function  in  the  integral  eliminates  that  part  of  the  integral  lying 
in  front  of  x^ 


HU-%)  = 


0  5  <  Xq  or  p  >  x(x-Xq) 


1  5  ix  or  p  £  p 


Thus,  the  upper  limit  becomes  p^,  thus 


P^.e.  "  Po 


X(X  -  Xq) 


and  the  term  H  (S-x^)  can  be  eliminated. 


Let 

p 

^(Po’  =  y*  e"'*P/^»  F(p)  dp 

o 

where  F(p)  =  Ji(p)/p 

Then  the  expression  for  w  becomes 


I  ^-^p  px-Xg)  + 

X[i/M^  +  \  E(Po.M„)]exp  {*U{x-x^)  /U^) j 


(C-5) 


(C-6) 
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i 


I 


The  force  f  generated  by  the  point  doublet  in  the  uniform  flow  is 


iCpdc 


/ 


'  ■  / 

Introducing  equation  (C-3)  gives 

f  *  -2A(() 


2  +  i  ^  A<^)  dc 


/ 


6(C-Xq)  d5  =  -  2a^p 


In  the  Doublet  Lattice  Method  a  point  dotblet  is  placed  on  each  element.  The 
resulting  force  generated  by  this  doublet  is  averaged  over  the  clement  to 
produce  an  average  pressure,  ^p.  Thus, 

^p  =  f/AX  =  -  ^  A(^p 


Solving  this  for  A<tip  and  introducing  the  result  into  {C-6)  gives  the  result 

for  the  local  downwash  w  due  to  a  point  pressure  doublet  of  average  pressure 

^  .  If  we  also  average  the  downwash  over  the  element,  vL,  then  the  result 

p  P 

is 

6  +  x[i/M^  +  E/2]  exp  ('iw(x-x^)  /  U^)  |  (C-7) 


w. 


Since  the  average  of  the  delta  function  is 

AX 

1 

6 


•-fe  / 


6(x-Xq)  dx 


Ji_ 

AX 


on  the  element 
off  of  the  element 


Equation  (C-7)  then,  along  with  the  expressions  for  E  and  6  are  the  final 
results  for  the  downwash  due  to  a  pressure  element  in  supersonic  flow.  The 
expression  for  E,  i.e.. 
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P  =  X(x-x  ) 

0 

F(p)  =  Jt(p)/P 


is  evaluated  numerically. 

The  steady  casei  u  =  0,  reduces  to  the  classic  results 


APPENDIX  D 


THE  GARABEDIAN  TRANSONIC  AIRFOIL  INVERSE  PROCEDURE 

The  program  method  of  analysis  is  based  on  the  work  of  Korn,  Garabedian 
and  Jameson  and  utilizes  a  Poisson  fast  solver  differencing  scheme.  Because 
viscous  effects  are  not  considered  in  the  2-D  Jameson  flow  solution,  the 
analytical  pressure  distribution  obtained  was  often  an  unacceptable  match 
with  pressures  obtained  experimentally.  The  effect  of  the  viscous  boundary 
layer  can  be  seen  in  Figure  2.4-5  where  the  pressures  on  both  the  upper  and 
lower  airfoil  surfaces  are  seen  to  deviate  from  the  experimental  data.  In 
order  to  account  for  the  effects  of  viscosity,  it  was  decided  to  alter  the 
geometry  of  the  airfoil  used  in  the  2-D  inviscid  analysis.  To  do  this,  use 
was  made  of  a  transonic  airfoil  design  computer  program  (ABBE)  containing 
an  iverse  Garabedian  design  option.  Although  basically  inviscid,  the 
program  accounts  for  viscosity  by  adding  a  displacement  thickness  to  the 
original  airfoil  geometry. 

To  decrease  computational  time  and  lower  cost,  a  Poisson  fast  solver 
differencing  scheme  was  incorporated  into  the  program.  Tlie  advantage  of  the 
fast  solver  lies  in  the  fact  that  convergence  occurs  in  fewer  cycles  through 
the  flow  field  than  with  other  methods.  Because  the  governing  compressible 
flow  equation  is  unstable  fcr  transonic  flows  when  written  in  the  form  used 
with  the  fast  solver,  fast  solver  cycles  are  alternated  with  successive  line 
over-relaxation  cycles.  The  SLOR  cycles  handle  the  supersonic  regions,  and 
the  fast  solver  cycles  propogate  information  quickly  through  the  flow  field 
and  satisfy  the  kutta  condition.  A  restriction  to  using  the  fast  Poisson 
solver  is  that  the  equations  must  be  solved  in  a  closed  area  with  Neuman, 
Dirichlit  or  periodic  boundary  conditions.  This  restriction  is  satisfied  by 
using  a  circle  as  the  computational  plane.  The  circumference,  which  corres¬ 
ponds  to  the  airfoil  surface,  has  Neumann  conditions,  the  center  of  the 
circle,  infinity  in  the  physical  plane,  has  Dirichlit  conditions,  and  periodic 
conditions  are  established  along  a  line  from  the  tail  to  infinity. 

To  aid  in  the  analysis,  the  program  was  set  up  to  run  from  the  IBM  2250 
graphics  terminal.  The  original  airfoil  geometry,  2-D  Mach  number  and  desired 
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2-D  lift  coefficient  are  input  and  a  20  cycle  flow  solution  is  performed 
first  on  a  coarse  grid  and  then  on  a  grid  refined  by  halving  the  crude  grid 
mesh  spacing.  The  resulting  pressure  distribution  is  then  displayed  on  the 
tube  and  is  changed  using  a  light  pen  to  correspond  to  the  experimental  data. 
Twenty  cycle  inverse  flow  solutions  are  then  performed  on  both  the  coarse  and 
refined  grids  to  obtain  a  revised  airfoil  geometry.  The  updated  geometry 
can  then  be  run  forward  through  a  flow  solution  and  the  resulting  pressure 
distribution  again  compared  and  updated  if  necessary.  This  process  is  repeated 
until  the  revised  geometry  produces  the  desired  pressure  distribution. 

Examples  of  airfoil  geometries  revised  using  the  inverse  Garabedian 
program  are  shown  in  Figure  3.3-3.  The  revision  to  the  station  .309  geometry 
is  a  typical  boundary  layer  displacement  correction  for  attached  flows  in  that 
the  trailing  edge  thickness  and  camber  are  altered. 

The  station  .657  revision  is  untypical  due  to  a  thickening  of  the  entire 
airfoil  section.  Because  the  two-dimensional  airfoil  geometry  is  scaled  from 
the  3-D  coordinates  as  Zgp/c  =  Z^q/c  /cos  A^pp*  the  thickness  correction 
indicates  the  use  of  an  incorrect  (too  small)  effective  sweep  angle  in  the 
initial  scaling  operation.  Because  the  2-D  Mach  number  is  scaled  down  as 
Moo  =  M  cosArcc.  the  thickness  displacement  can  also  be  seen  to  indicate 

Cu  *  Lrr  , 

that  too  large  a  2^D  Mach  number  was  used  in  the  inverse  calculation.  This 
uniform  thickening  would  be  eliminated  if  the  2-D  Mach  number  were  lowered 
an  appropriate  amount  and  the  case  rerun  in  the  Garabedian  inverse. 
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